Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,267,395 A→C 19.2% S192R (AGC→CGC)  rlmG → Ribosomal RNA large subunit methyltransferase G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,267,3950AC19.2% 129.4 / 3.6 71S192R (AGC→CGC) rlmGRibosomal RNA large subunit methyltransferase G
Reads supporting (aligned to +/- strand):  ref base A (17/38);  new base C (13/3);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 4.70e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCTGCCGCGCGCGGCCGGCGACCTGATGGAGAAGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCGGG  >  ACB122/5267263‑5267524
                                                                                                                                    |                                                                                                                                 
cACCGGCCGCGCGCGGCCGGCGACCTGATGGAGAATTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAgc                                                                                                                                 <  1:451841/135‑1 (MQ=255)
 aCCTGCCGCGCGCGGCCGGCGACCTGATGGAGAAGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGCGCGGCCGCTCGCCCAGTCCCCCTGCCCCCgct                                                                                                                                >  3:341807/1‑135 (MQ=255)
       cgcgcgcgGCCGGCGACCTGATGGAGAAGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCACCGGTTGCCGCGCGGCCCCTCGCCCAGTGCCCCTACCCCCGCGGCTAcc                                                                                                                          >  7:366240/1‑135 (MQ=255)
           cgcgGCCGGCGACCTGATGGAGAAGTCCATCGGCCCCGTACAGGCTTCGCGGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCGCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCt                                                                                                                      <  8:309945/135‑1 (MQ=255)
           cgcgGCCGGCGACCTGATGGAGAAGTACATCGGCCCCGTACATGCGGCGCTGGCCCTGAAAAAGGCCCGTTGGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGccg                                                                                                                      <  6:217000/135‑2 (MQ=255)
           cgcgGCCGGCGACCTGATGGAGAAGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCt                                                                                                                      <  6:224244/135‑1 (MQ=255)
            gcgGCCGGCGCCCTGATGGAGAAGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCCACGGTTGCCGCGCGGCCGCGCGCCCAGGCCCCCGACCCCCCCGGCGACCGccgc                                                                                                                     >  7:445502/1‑133 (MQ=255)
            gcgGCCGGCGACCTGATGGAGAAGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCCGCTGCTACCGCCTc                                                                                                                     >  6:94769/1‑135 (MQ=255)
                ccggcGACCTGATGGAGAAGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAACGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGAcg                                                                                                                 <  7:242149/135‑1 (MQ=255)
                ccggcGACCTGATGGAGAAGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGAcg                                                                                                                 <  8:346467/135‑1 (MQ=255)
                     gACCTGATGTCGAAGTCCATCGGCCCCGTACAGGCTTCGCGGGCACTGAAAACGGCCCGTTTGCGGACAGCAACGGTTGCCGAGCGGCCGCGCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCg                                                                                                            <  3:217573/135‑1 (MQ=255)
                            tGGAGAAGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGCGCGGCCGCTCGCCCAGTCCCCCAACCCCCGCGGCTACCGCCTCGCCGCGCCGGCccttg                                                                                                     >  7:63578/1‑135 (MQ=255)
                                 aaGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGCCAGCAACGGTTGCCGCGCGGCCGCTCGCCAAGTCCCCCAACCCCACCTGCAACCCCCTCCCCCCCCCGGCccatgacttg                                                                                                >  2:218675/1‑135 (MQ=255)
                                   gTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGt                                                                                              <  8:68110/135‑1 (MQ=255)
                                    tACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCACCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGt                                                                                              >  4:432528/1‑134 (MQ=255)
                                     aCATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCa                                                                                            <  7:417115/135‑1 (MQ=255)
                                           gCCCCTTACAGGCTTGGGTGGCACTGAAAACGGCCGGTTTGCTGACTCCAAGGTTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACg                                                                                      <  3:17358/135‑1 (MQ=255)
                                              ccGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCACCGGTTGCCGCGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCAGCTACCGCCTCGACGCGCCGGCACGTGACGTGGTCAACCACGccc                                                                                   >  5:113606/1‑134 (MQ=255)
                                              ccGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGCCTTGGTCAACCCCGCCa                                                                                   >  4:112974/1‑135 (MQ=255)
                                              ccGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGACCCGCCCCTCGCCAACTCCCCCAACCCCAGCGGCTACCGCCGCGACCCGCCGGCCCTTTAATTTGTCAACCACGCCa                                                                                   >  8:23466/1‑135 (MQ=255)
                                                gTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAAc                                                                                 <  6:109653/135‑1 (MQ=255)
                                                gTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAAc                                                                                 <  2:23285/135‑1 (MQ=255)
                                                 tACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGGCGCTCCCCAAGTACCCCGACCCCCGCTGCTACCGCCTCGACGCGCCCGCACTTGACTTGGTCAACCACaacaacc                                                                                >  3:145432/1‑134 (MQ=255)
                                                     ggCTGCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAAGGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGtt                                                                            <  1:272608/135‑1 (MQ=255)
                                                         tCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGCGCGGCCGCTCGCCCAGTCCCCCGACCCCACCTGCGACCCCCGCGCCGCGCCGGCACGTGACTTGGTCAACCACCCCCACGTGTTATGc                                                                        >  7:430696/1‑135 (MQ=255)
                                                         tCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGCGCGGCCGCTAGCCAAGTCCCCAAACCCCAGCTGCTACCGCCGCGACGCGCCGGCACATGACGTTGTACACCACGCCAAAGTGTTCTGc                                                                        >  4:23394/1‑135 (MQ=255)
                                                         tCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCCAGTCCCCCTACCCCCGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGc                                                                        >  6:365377/1‑135 (MQ=255)
                                                         tCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCCGCTGCTACCGCCTCGACGCGCCGGCCCTTGACTTGGTCAACCCCGCCAACGGGTTCTGc                                                                        >  4:228381/1‑135 (MQ=255)
                                                           ggtGGCAATCAAACAGGCCCGTTTCATGACAGAAAGGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGcc                                                                       <  7:78087/132‑1 (MQ=255)
                                                           gCTGGCACTGAAAAAGGCCTGTTTGCTGACAGCAAGGGTTGCCGAGCGGCCGCTCGCCAAGTCCGCATACCCCAGCTGCTACCGCATCGACGCGCCGGCACTTGTCTTGGTGAACCACGCCAACGTGTTCTgac                                                                       <  8:459872/134‑3 (MQ=255)
                                                             tGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGt                                                                     >  4:337731/1‑134 (MQ=255)
                                                             tGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGt                                                                     >  2:377309/1‑134 (MQ=255)
                                                             tGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCGCGCCCAGTCCCCCTACCCCCCCAGCTACCGCCGCGCCCCCCCGGCCCTGGAGTGGGTGTACAACGCCCACGGGTTCTGGCGt                                                                     >  1:310048/1‑134 (MQ=255)
                                                              ggCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGCGCGGCCGCTCGCCAAGTCCCCCGACGCCAGCTGCTACCGCCGCGCCGCGCCGGCCCGTGCCTTGGTCCACCACGCCCACCGGTTCTGGCgg                                                                     >  5:398491/1‑132 (MQ=255)
                                                              ggCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGCGCGGCCGCTCGCCAAGTCCCCCGACCCCCGCTGCTACCGCCTCGACGCGCCGGCACGTGACTTGGTCAACCCCCCCCACGTGGTCTGCCgg                                                                     >  2:17802/1‑132 (MQ=255)
                                                              ggCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGATCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCCCTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGt                                                                     >  2:366283/1‑133 (MQ=255)
                                                              ggCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCCCGCGCCCAGTCCCCCTACCCCCGCTGCTACCGCCGCGACGCGCCCGCACTTGACTAGGTCAACCACGCCAACGGGTTCTGCCgg                                                                     >  3:397311/1‑132 (MQ=255)
                                                                cACTGAAGAAGGACCGTTTGCTGCAAGCACCGGTTGCCGAGCGGCCGCGCGCCAAGGCCCCCTACCCCAGCGGATACCGCATCGACGCGCCTGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTTaa                                                                  <  6:420655/134‑1 (MQ=255)
                                                                 aCTGAAAAAGGCCCTTTTGCTGACAGCAACGGTTGCCGAGCTGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAgg                                                                <  6:224379/135‑1 (MQ=255)
                                                                 aCTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAgg                                                                <  5:230671/135‑1 (MQ=255)
                                                                   tGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGCGCGGCCGCTCGCCAACTCCCCCGACCCCCCCAGCTACCGCCCCCACCCGCCCGCCCTTGACTTTGTGAACCCCCCCAACCTGTTGTGCCGTGAAGGcc                                                              >  3:4257/1‑135 (MQ=255)
                                                                   tGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTGCCCCCGCTGCTACCGCCTCGACGCGCCGGCCCTTGACTTGGTCAACCCCGCCACCGTGGTCTGCCGTGAAGGcc                                                              >  2:99221/1‑135 (MQ=255)
                                                                   tGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGCCTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGcc                                                              >  2:174245/1‑135 (MQ=255)
                                                                       aatggCCCGTTTGCGGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGa                                                          <  8:176442/132‑1 (MQ=255)
                                                                                  tGCTGACAGCAACGGTTGCCGCGCGGCCGCTCGCCAAGTCCCCCGACCCCAGCTGCTACCGCCGCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGAcgc                                               >  6:199631/1‑135 (MQ=255)
                                                                                  tGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCGACCCCAGCTGCTACCGCCTCGCCGCGCCGGCACGTGACGTGGTCAACCCCGCCCACGGGTTCTGCCGGGAAGGCCTGGATATCGGGAcgc                                               >  4:96129/1‑135 (MQ=255)
                                                                                   gCTGACAGCAACGGTTGCAGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGACGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGAcgc                                               <  3:355408/134‑1 (MQ=255)
                                                                                     tGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGGCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCggg                                            <  8:445502/135‑1 (MQ=255)
                                                                                        cAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGc                                           <  6:328399/133‑1 (MQ=255)
                                                                                           cAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTtcc                                      <  8:430696/135‑1 (MQ=255)
                                                                                                      gAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGCCGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTGCGCCCGCCTCTa                           >  5:392106/1‑135 (MQ=255)
                                                                                                       aGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTAc                          <  6:223542/135‑1 (MQ=255)
                                                                                                        gCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTAcc                         <  1:173468/135‑1 (MQ=255)
                                                                                                         cGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCa                        <  2:298325/135‑1 (MQ=255)
                                                                                                         cGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCa                        <  3:320745/135‑1 (MQ=255)
                                                                                                         cGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCa                        <  1:400334/135‑1 (MQ=255)
                                                                                                              ggtCGCCAAGTCCCCCTCCCCCCGGTGCCCCCCCCGCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCaa                   <  1:436637/133‑1 (MQ=255)
                                                                                                                 cGCCAAGTCCCCCAACCCCAGCTGCTACCGCCTCGCCGCGCCGGCCCTAGACTTCGTCAACCACGCCAACGGGGTCTGCCGGGAAGGCCTGGATATCGGGCCGCGGGCGGTCCGCCCGCATCTACCACGCAACCt                >  7:369070/1‑135 (MQ=255)
                                                                                                                          ccccTACCCCAGCTGCTCCCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTgc       <  7:458667/135‑1 (MQ=255)
                                                                                                                          ccccTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTgc       <  8:63578/135‑1 (MQ=255)
                                                                                                                          ccccTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTgc       <  6:98306/135‑1 (MQ=255)
                                                                                                                          ccccTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGACAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTgc       <  4:159002/135‑1 (MQ=255)
                                                                                                                           cccTACCCCCGCTGCTACCGCCGCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTgc       <  7:452823/134‑1 (MQ=255)
                                                                                                                               cccccCGCTGCTCCCGCCGCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCggg  <  7:398179/134‑1 (MQ=255)
                                                                                                                               aCCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCggg  >  2:144978/1‑135 (MQ=255)
                                                                                                                                 cccAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCggg  <  5:199631/133‑1 (MQ=255)
                                                                                                                                 cccAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCggg  <  5:391414/133‑1 (MQ=255)
                                                                                                                                 cccAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCggg  <  2:334883/133‑1 (MQ=255)
                                                                                                                                 cccAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCggg  <  1:174245/133‑1 (MQ=255)
                                                                                                                                 cccAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACACCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCggg  <  4:145432/133‑1 (MQ=255)
                                                                                                                                 cccAGCGGCTCCCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCggg  <  1:409693/133‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                 
CACCTGCCGCGCGCGGCCGGCGACCTGATGGAGAAGTACATCGGCCCCGTACAGGCTTCGCTGGCACTGAAAAAGGCCCGTTTGCTGACAGCAACGGTTGCCGAGCGGCCGCTCGCCAAGTCCCCCTACCCCAGCTGCTACCGCCTCGACGCGCCGGCACTTGACTTGGTCAACCACGCCAACGTGTTCTGCCGTGAAGGCCTGGATATCGGGACGCGGGCGTTCCTCCCGCATCTACCACGCAACCTCGGCCGTGCGCGGG  >  ACB122/5267263‑5267524

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: