Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,168,351 C→T 8.3% 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.3% 198.6 / 4.1 71intergenic (+151/‑21)rrsB/gltT16S ribosomal RNA of rrnB operon/tRNA‑Glu
Reads supporting (aligned to +/- strand):  ref base C (30/35);  new base T (3/3);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTC  >  NC_000913/4168203‑4168495
                                                                                                                                                         |                                                                                                                                                
ttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAAc                                                                                                                                                  >  2:465154/1‑149 (MQ=39)
ttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAAc                                                                                                                                                  >  4:351825/1‑149 (MQ=39)
  aaaGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgc                                                                                                                                                <  3:72345‑M2/149‑3 (MQ=255)
  aaaGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACg                                                                                                                                                 >  6:7725‑M2/1‑147 (MQ=255)
    agaagCGTACTTTGTAGTGCTCACACAGATTGTTTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              <  3:234344‑M2/149‑5 (MQ=255)
    agaagCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAGGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              <  6:74627‑M2/149‑5 (MQ=255)
    agaagCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAGGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              >  5:74627‑M2/1‑145 (MQ=255)
    agaagCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCCACCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              <  7:412586‑M2/149‑5 (MQ=255)
    agaagCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              <  4:175548‑M2/149‑5 (MQ=255)
      aagCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCGTTGCGTGa                                                                                                                                            <  3:255357/149‑1 (MQ=255)
          gTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAg                                                                                                                                        <  5:131271/149‑1 (MQ=39)
          gTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAg                                                                                                                                        >  4:75757/1‑149 (MQ=39)
          gTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAg                                                                                                                                        >  1:321382/1‑149 (MQ=39)
            aCTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGc                                                                                                                                       >  5:20390/1‑148 (MQ=39)
             cTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCaa                                                                                                                                     <  2:321382/149‑1 (MQ=39)
                 gTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAAtttt                                                                                                                                 >  5:25870/1‑149 (MQ=38)
                 gTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAAttt                                                                                                                                  <  2:5587/148‑1 (MQ=38)
                  tAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTc                                                                                                                                <  4:388794/149‑1 (MQ=38)
                  tAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTc                                                                                                                                <  8:206701/149‑1 (MQ=38)
                        tCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTcc                                                                                                                          >  6:263735/1‑149 (MQ=37)
                         cacacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTccc                                                                                                                         >  1:175845/1‑149 (MQ=37)
                          acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTcccc                                                                                                                        >  1:292581/1‑149 (MQ=37)
                          acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTcccc                                                                                                                        >  6:42152/1‑149 (MQ=37)
                          acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTcc                                                                                                                          <  6:217224/147‑1 (MQ=37)
                              aGATTGTCTGATAGAAAGTGAGAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCg                                                                                                                    <  7:80556/149‑1 (MQ=37)
                               gATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGt                                                                                                                   <  1:465154/149‑1 (MQ=37)
                                aTTGTCTGATATAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTTAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACTCGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTc                                                                                                                  <  3:80370/149‑1 (MQ=255)
                                    tttGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTaga                                                                                                              <  4:234344/147‑1 (MQ=37)
                                    tCTGATGAAAAGTGAATAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCGTTgcgagatattttcgtgtagggtgagctttccttaa                                                                                                               <  2:14857‑M2/148‑36 (MQ=255)
                                    tCTGATGAAAAGTGAATAGCAAGGCGATTACGCGTTGGGAGTCAGGCTGAAGAGTATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATGACGCGTTgcgtgatattttcgtgtagggtgagctttcattaa                                                                                                               >  1:14857‑M2/1‑113 (MQ=255)
                                    tCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTaga                                                                                                              <  5:256833/149‑1 (MQ=37)
                                      tGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTcc                                                                                                                          >  3:318810/1‑135 (MQ=37)
                                      tGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTcc                                                                                                                          <  4:318810/135‑1 (MQ=37)
                                      tGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAgg                                                                                                            >  8:425397/1‑149 (MQ=37)
                                       gATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGc                                                                                                           >  1:185276/1‑149 (MQ=35)
                                        aTAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTagag                                                                                                             <  8:354391/146‑1 (MQ=35)
                                         tAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGcc                                                                                                          >  4:318258/1‑148 (MQ=34)
                                          aGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              >  2:55316‑M2/1‑107 (MQ=255)
                                          aGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              <  1:55316‑M2/111‑5 (MQ=255)
                                           gAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCa                                                                                                        <  6:270338/148‑1 (MQ=32)
                                           gAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAg                                                                                                       <  5:416725/149‑1 (MQ=32)
                                             aaGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGa                                                                                                     <  2:185276/149‑1 (MQ=32)
                                             aaGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGa                                                                                                     <  3:351825/149‑1 (MQ=32)
                                                 gAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACAcc                                                                                                 >  1:213992/1‑149 (MQ=32)
                                                  aaaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACAcc                                                                                                 <  6:25870/148‑1 (MQ=32)
                                                  aaaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCg                                                                                                <  7:261431/149‑1 (MQ=32)
                                                   aaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCTATTTTCGTGTCCCCTTCGTCTTGAGGCCCAGCACAcccc                                                                                               >  2:438817/1‑147 (MQ=9)
                                                         agtcGTTTACGCGTTGGGAGTGAGGCTGAAGATAATAAGGCCGTTCGCTTTCTATTAATTAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTc                                                                                         <  5:96449/146‑1 (MQ=0)
                                                                 aCGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGATCACCCTACACGAAAATTGCTTATCACGCGTTgcgtgatattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctc                                                                                       >  6:415359‑M2/1‑89 (MQ=255)
                                                                   gcgTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgtgat                                                                                                                                           >  1:49124‑M2/1‑82 (MQ=255)
                                                                   gcgTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgtgat                                                                                                                                           <  2:49124‑M2/89‑8 (MQ=255)
                                                                                   tGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgtgata                                                                                                                                          >  4:181171‑M2/1‑66 (MQ=255)
                                                                                       gagaATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCGTTgcgtgatattttcgtgtagggtgagctttcattaa                                                                                                               >  3:398115‑M2/1‑62 (MQ=255)
                                                                                       gagaATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCGTTgcgtgatattttcgtgtagggtgagctttcattaa                                                                                                               <  4:398115‑M2/97‑36 (MQ=255)
                                                                                                            ctcttcaTGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCACCGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGttt                                       <  8:412586/142‑1 (MQ=0)
                                                                                                                                tACACGAAAAT‑‑‑‑‑ATCACGCAACgcgtgataagcaattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgctattcacttttcatcagacaatctgtgt                     <  1:434900‑M2/146‑126 (MQ=255)
                                                                                                                                  cacGAAAAT‑‑‑‑‑ATCACGCAACgcgtgataagcaattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgctattcacttttcatcagacaatctgtgtga                   >  2:434900‑M2/1‑19 (MQ=255)
                                                                                                                                   acGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTaa               >  7:54138/1‑149 (MQ=32)
                                                                                                                                     gAAAAT‑‑‑‑‑ATCACGCAACgcgtgataagcaattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgcttttcactttc                                      >  3:375386‑M2/1‑16 (MQ=255)
                                                                                                                                     gAAAAT‑‑‑‑‑ATCACGCAACgcgtgataagcaattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgcttttcactttc                                      <  4:375386‑M2/124‑109 (MQ=255)
                                                                                                                                     gAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAata             >  6:357819/1‑149 (MQ=32)
                                                                                                                                     gAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAata             >  7:164199/1‑149 (MQ=32)
                                                                                                                                     gAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAata             >  5:399049/1‑149 (MQ=32)
                                                                                                                                     gAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAata             >  5:375500/1‑149 (MQ=32)
                                                                                                                                        aaT‑‑‑‑‑ATCACGCAACgcgtgatattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgctattcacttttc                                          >  6:245074‑M2/1‑13 (MQ=255)
                                                                                                                                        aaT‑‑‑‑‑ATCACGCAACgcgtgatattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgctattcacttttc                                          <  5:245074‑M2/117‑105 (MQ=255)
                                                                                                                                         at‑‑‑‑‑atCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATAtctc         <  5:357819/149‑1 (MQ=32)
                                                                                                                                                     gCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACt   >  5:412010/1‑148 (MQ=32)
                                                                                                                                                     gCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACt   >  4:365721/1‑148 (MQ=32)
                                                                                                                                                     gCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTc  <  7:180116/149‑1 (MQ=33)
                                                                                                                                                     gCAACGCGTGATAAGCAATTTTAGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTc  <  3:318258/149‑1 (MQ=18)
                                                                                                                                                       aaCGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCaaaa     <  6:187765/144‑1 (MQ=32)
                                                                                                                                                       aaCGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCaaaa     >  5:187765/1‑144 (MQ=32)
                                                                                                                                                         |                                                                                                                                                
TTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTC  >  NC_000913/4168203‑4168495

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