Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,168,351 C→T 8.7% intergenic (+151/‑21) 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,3510CT 8.7% 195.5 / 9.3 69intergenic (+151/‑21)rrsB/gltT16S ribosomal RNA of rrnB operon/tRNA‑Glu
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base T (4/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCA  >  NC_000913/4168203‑4168496
                                                                                                                                                    |                                                                                                                                                  
ttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAc                                                                                                                                                    >  3:779673/1‑149 (MQ=255)
ttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAc                                                                                                                                                    >  3:141640/1‑149 (MQ=255)
ttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAc                                                                                                                                                    >  4:762546/1‑149 (MQ=39)
      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                                                                                                                                       >  4:141640‑M2/1‑148 (MQ=255)
             cTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAa                                                                                                                                       >  3:1868062‑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                                                                                                                                            <  2:1786380/131‑1 (MQ=255)
                          acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCGTTGCGTGATa                                                                                                                                            >  1:1786380/1‑131 (MQ=255)
                          acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCGTTGCGTGATa                                                                                                                                            <  4:1310530/131‑1 (MQ=255)
                          acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCGTTGCGTGATa                                                                                                                                            >  3:1310530/1‑131 (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                                                                                                                     >  2:636197/1‑149 (MQ=39)
                               gATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGt                                                                                                                     >  3:416499/1‑149 (MQ=38)
                                    tCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgt                                                                                                                                >  4:1876204‑M2/1‑125 (MQ=255)
                                    tCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgt                                                                                                                                <  3:1876204‑M2/133‑9 (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                                                                                                                                >  1:998219‑M2/1‑92 (MQ=255)
                                                                     gTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAattttcgt                                                                                                                                <  2:998219‑M2/100‑9 (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                                  <  4:784605‑M2/132‑103 (MQ=255)
                                                                                                                                   acGAAAATATCACGCAACGCGTGATAAGCAattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgcttttcactttctatcag                                  >  3:784605‑M2/1‑30 (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        <  1:1680486/142‑1 (MQ=32)
                                                                                                                                                  aaCGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCaa        >  2:1680486/1‑142 (MQ=32)
                                                                                                                                                   aCGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCa  <  2:827828/147‑1 (MQ=34)
                                                                                                                                                   aCGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCa  >  1:827828/1‑147 (MQ=35)
                                                                                                                                                    |                                                                                                                                                  
TTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCA  >  NC_000913/4168203‑4168496

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: