Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,168,349 A→G 8.8% intergenic (+149/‑23) rrsB → / → gltT 16S ribosomal RNA of rrnB operon/tRNA‑Glu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,168,3490AG 8.8% 195.9 / 9.2 68intergenic (+149/‑23)rrsB/gltT16S ribosomal RNA of rrnB operon/tRNA‑Glu
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base G (4/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAA  >  NC_000913/4168203‑4168492
                                                                                                                                                  |                                                                                                                                                
ttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAc                                                                                                                                                >  3:779673/1‑149 (MQ=255)
ttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAc                                                                                                                                                >  4:762546/1‑149 (MQ=39)
ttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAc                                                                                                                                                >  3:141640/1‑149 (MQ=255)
      aagCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGa                                                                                                                                          <  4:779673/149‑1 (MQ=255)
      aagCGTACTTTGTAGTGCTCACACAGATTGACTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGa                                                                                                                                          <  3:303820/149‑1 (MQ=255)
        gCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATa                                                                                                                                        >  4:1868062/1‑149 (MQ=40)
             cTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGa                                                                                                                                          >  1:1483026/1‑142 (MQ=40)
             cTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGa                                                                                                                                          <  2:1483026/142‑1 (MQ=39)
             cTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAa                                                                                                                                   >  3:1868062‑M2/1‑148 (MQ=255)
             cTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAa                                                                                                                                   >  4:141640‑M2/1‑148 (MQ=255)
                  tAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttc                                                                                                                              >  1:1829950‑M2/1‑143 (MQ=255)
                     tGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgtg                                                                                                                           <  2:1426275‑M2/149‑10 (MQ=255)
                       cTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATaa                                                                                                                                       >  4:181464/1‑135 (MQ=39)
                       cTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATaa                                                                                                                                       <  3:181464/135‑1 (MQ=39)
                        tCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCGTTGCGTGATa                                                                                                                                        >  1:1193915/1‑133 (MQ=255)
                        tCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGa                                                                                                                                          >  3:1383206/1‑131 (MQ=39)
                        tCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGa                                                                                                                                          <  4:1383206/131‑1 (MQ=38)
                        tCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTcc                                                                                                                        >  3:809657/1‑149 (MQ=38)
                          acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCGTTGCGTGATa                                                                                                                                        <  4:1310530/131‑1 (MQ=255)
                          acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCGTTGCGTGATa                                                                                                                                        >  1:1786380/1‑131 (MQ=255)
                          acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCGTTGCGTGATa                                                                                                                                        >  3:1310530/1‑131 (MQ=255)
                          acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCGTTGCGTGATa                                                                                                                                        <  2:1786380/131‑1 (MQ=255)
                          acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTcccc                                                                                                                      >  1:1834382/1‑149 (MQ=39)
                           cacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCt                                                                                                                     >  2:677343/1‑149 (MQ=39)
                            acaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATa                                                                                                                                        <  3:117482/129‑1 (MQ=40)
                            acaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATa                                                                                                                                        >  4:117482/1‑129 (MQ=40)
                              aGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCg                                                                                                                  <  1:1426275/149‑1 (MQ=39)
                               gATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGt                                                                                                                 >  3:416499/1‑149 (MQ=38)
                               gATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGt                                                                                                                 >  2:636197/1‑149 (MQ=39)
                                    tCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgt                                                                                                                            <  3:1876204‑M2/133‑9 (MQ=255)
                                    tCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgt                                                                                                                            >  4:1876204‑M2/1‑125 (MQ=255)
                                      tGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttc                                                                                                                              >  2:1713383‑M2/1‑123 (MQ=255)
                                      tGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttc                                                                                                                              <  1:1713383‑M2/129‑7 (MQ=255)
                                        aTAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACTCGAAAATATCACGCAACGCGt                                                                                                                                            <  3:903421/113‑1 (MQ=37)
                                        aTAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACTCGAAAATATCACGCAACGCGt                                                                                                                                            >  4:903421/1‑113 (MQ=38)
                                             aaGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGt                                                                                                                                            <  4:131546/108‑1 (MQ=32)
                                             aaGTGAAAAGCAAGACGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGt                                                                                                                                            >  3:131546/1‑108 (MQ=18)
                                                  aaaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTcc                                                                                                                        <  4:1286376/123‑1 (MQ=32)
                                                  aaaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCATACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTcc                                                                                                                        >  3:1286376/1‑123 (MQ=17)
                                                    aaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGcc                                                                                            <  2:1834382/149‑1 (MQ=11)
                                                         agacGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTATCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGc                                                                                                         <  1:1448365/128‑1 (MQ=1)
                                                                     gTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgt                                                                                                                            <  2:998219‑M2/100‑9 (MQ=255)
                                                                     gTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgt                                                                                                                            >  1:998219‑M2/1‑92 (MQ=255)
                                                                      ttGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgt                                                                                                                            >  2:887494‑M2/1‑91 (MQ=255)
                                                                      ttGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgt                                                                                                                            <  1:887494‑M2/99‑9 (MQ=255)
                                                                                          aaTAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTGGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGATTCGAACCCCtgtt                                                       >  3:1848031/1‑144 (MQ=1)
                                                                                          aaTAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTGGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGATTCGAACCCCtgtt                                                       <  4:1848031/147‑4 (MQ=1)
                                                                                                                       aGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgtg                                                                                                                           <  2:84607‑M2/51‑10 (MQ=255)
                                                                                                                       aGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgtg                                                                                                                           >  1:84607‑M2/1‑42 (MQ=255)
                                                                                                                              ccTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcc                                                                  >  3:1182953‑M2/1‑35 (MQ=255)
                                                                                                                              ccTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcc                                                                  <  4:1182953‑M2/101‑67 (MQ=255)
                                                                                                                               cTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaa                                                        >  2:705062‑M2/1‑34 (MQ=255)
                                                                                                                               cTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaa                                                        <  1:705062‑M2/110‑77 (MQ=255)
                                                                                                                                tACACGAAAATATCACGCAACGCGTGATAAGCAattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgcttttcactttctatcagacaatctgtgtgag                <  1:563572‑M2/149‑117 (MQ=255)
                                                                                                                                 acacGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCtt              >  1:1080621/1‑149 (MQ=11)
                                                                                                                                  cacGAAAATATCACGCAACGCGTGATAAGCAattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgcttttcactttctatcagacaatctgtgtga                 >  1:771508‑M2/1‑31 (MQ=255)
                                                                                                                                  cacGAAAATATCACGCAACGCGTGATAAGCAattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgcttttcactttctatcagacaatctgtgtga                 <  2:771508‑M2/146‑116 (MQ=255)
                                                                                                                                  cacGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGTGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTa             <  1:1342681/149‑1 (MQ=14)
                                                                                                                                   acGAAAATATCACGCAACGCGTGATAAGCAattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgcttttcactttctatcag                              >  3:784605‑M2/1‑30 (MQ=255)
                                                                                                                                   acGAAAATATCACGCAACGCGTGATAAGCAattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgcttttcactttctatcag                              <  4:784605‑M2/132‑103 (MQ=255)
                                                                                                                                   acGAAAATATCACGCAACGCGTGATAAGCAattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgctattcacttttcatcagacaatctgtgtgagcac             >  3:741620‑M2/1‑30 (MQ=255)
                                                                                                                                   acGAAAATATCACGCAACGCGTGATAAGCAattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgctattcacttttcatcagacaatctgtgtgagcac             <  4:741620‑M2/149‑120 (MQ=255)
                                                                                                                                   acGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTaa            >  1:1746247/1‑149 (MQ=32)
                                                                                                                                     gAAAATATCACGCAACGCGTGATA‑‑‑tttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgcttttcactttctatcagacaatctgt                        <  1:1658625‑M2/134‑110 (MQ=255)
                                                                                                                                      aaaaTATCACGCAACGCGTGATAAGCAattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgcttttcactttctatcagacaatctgtgtgagcactac          <  1:1140288‑M2/149‑123 (MQ=255)
                                                                                                                                        aaTATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATAtct       <  2:1746247/149‑1 (MQ=32)
                                                                                                                                         ttatCACGCGTTGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATAtctc      >  1:1679052/2‑149 (MQ=12)
                                                                                                                                             cACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCaaaa  <  2:1405627/149‑1 (MQ=32)
                                                                                                                                                  aaCGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCaa    >  2:1680486/1‑142 (MQ=32)
                                                                                                                                                  aaCGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCaa    <  1:1680486/142‑1 (MQ=32)
                                                                                                                                                  |                                                                                                                                                
TTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAA  >  NC_000913/4168203‑4168492

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