Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,561 G→C 7.6% intergenic (‑579/+24) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,5610GC 7.6% 165.4 / 3.4 52intergenic (‑579/+24)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (25/23);  new base C (2/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

CTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAGCCACCAATGATCTCACCGATGCCCGGAGCCAGAACGTCCATAGCCGCAA  >  NC_000913/987261‑987856
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                       
ttATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCggg                                                                                                                                                                                                                                                                                                         <  2:298030/300‑1 (MQ=255)
 tATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAACAATGCGGATTAGTTTTTGCGCGTAAATGTTGCCATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCCAATAAACTTCTTTCCATAACAAAACAAAAACCTTACAGAAAAGAAAAAAAAAAAACCAACCGgc                                                                                                                                                                                                                                                                                                        >  1:515919/1‑301 (MQ=255)
           tcagacgtgtgctcttccgatcTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtt                                                                                                                                                                                                                                                                                              <  2:132862/280‑1 (MQ=255)
               ccTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAg                                                                                                                                                                                                                                                                                          <  2:351248/301‑1 (MQ=255)
                     atTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAAAAAAAAAGAACAAGCGGGCGCAGGCTCTTGAAGACGCtc                                                                                                                                                                                                                                                                                    >  1:281895/1‑300 (MQ=255)
                          ttcgatcTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTa                                                                                                                                                                                                                                                                               <  1:475966/295‑1 (MQ=255)
                                  aaaCGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTTACTTAAAAAACCGCTAACTTAACAGAAAAGATAAAAAAAAAACCAGCACCCGCTGGCCCttagaacagaagaacgtcatgtaag                                                                                                                                                                                                                                                                       >  1:132737/1‑278 (MQ=255)
                                  aaaCGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCAAGCGGGCGCTGGCTCTTGAAGACGCTAAATTaagaacgga                                                                                                                                                                                                                                                                       >  2:476567/1‑293 (MQ=255)
                                         taTTTTGTTTCTGATGTAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcg                                                                                                                                                                                                                                                                <  2:101618/300‑1 (MQ=255)
                                          aTTTTGTTTCTGACGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgc                                                                                                                                                                                                                                                               >  2:547619/1‑301 (MQ=255)
                                            tttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                                                                                                                                                                                                                             <  2:256663/301‑1 (MQ=255)
                                                        ggAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGCAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCCAATAGATTGCTTGACTTAACACACAGATAAATTGACAGAAAAGATAAAAAAACAGCCAACGGGCCATGGCCCTTGCAGACGCCAAATTAAAAACTGGCATTACGCCGAATACCTGGGaa                                                                                                                                                                                                                                                 >  1:24613/1‑301 (MQ=255)
                                                                      ccTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGCTAACTTGACAGAAAAGAAAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAAAAACTGGCGGTACCCAGAATACGGAGGAACCGaaaccaaaacc                                                                                                                                                                                                                                   >  1:212255/1‑292 (MQ=255)
                                                                                               aTCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCCCTTGAAGACGCCAAATAAGAAACTGGCGTTACGCCGAGTAAGTGGGAACGGAATCACATCACGTCACTTTTGCACGCCCGTTagctt                                                                                                                                                                                                          >  1:3911/1‑297 (MQ=255)
                                                                                                 ccAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAg                                                                                                                                                                                                        <  2:3911/301‑1 (MQ=255)
                                                                                                  cAGAATGCTAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGc                                                                                                                                                                                                       <  2:286857/301‑1 (MQ=255)
                                                                                                       tGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCa                                                                                                                                                                                                  <  2:282198/301‑1 (MQ=255)
                                                                                                              aaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGGTTGAATTAAACAAAAAAAAACTTGACAGAAAAAATAAAAAAAGAACCAGCGGGCGGCGGCCCTTGGAGACGCCAAAATAGAAGCCGGCGGTAAGCCGAAGAAGGTGGAACGGGAACACAACAACTACGTTTTTGAAGCCAAGTACATAAACaaaaaaaaattt                                                                                                                                                                                           >  2:402664/1‑291 (MQ=255)
                                                                                                                                      aGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCGAATAGATTGCTTGACTTAACAAACACATCACTTCACAGAAAAGGTAAAAAAAGAGCCCGCGGGCGCCGGCTCTTGAAGACGCCACAGTAGAAGCAGACGTTACGCGGAGTACGTGGGAACGCAATCCCATCACCTACGTTTCTCACGACAATTAACTAAACAAACAAACGATCAAAACCaaaaaccaaaaccacaaa                                                                                                                                                                   >  1:74726/1‑285 (MQ=255)
                                                                                                                                            aaaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAAc                                                                                                                                                             <  1:546946/300‑1 (MQ=255)
                                                                                                                                              aaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGTGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAAACCGAAACCTGAATGCGAAACaa                                                                                                                                                           >  2:628234/1‑300 (MQ=255)
                                                                                                                                                   cTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACAACACGTACGTTTTTCACGCCAATTACATAAGAAATCAGACGTTCAAAAACAAGACCAAAAcacaaaaaaagaaaaataccg                                                                                                                                                      >  1:47208/1‑281 (MQ=255)
                                                                                                                                                     aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAAAAGTTCAAAACCAAGAACGAAACCAGAAAGCGGAACAGTaacctt                                                                                                                                                    >  1:213412/1‑296 (MQ=255)
                                                                                                                                                       ttcGCAATAAGTATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGc                                                                                                                                                  <  2:472253/299‑1 (MQ=255)
                                                                                                                                                       tCCGAAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATAGTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGc                                                                                                                                                  <  2:624703/301‑1 (MQ=255)
                                                                                                                                                           cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCCGTTACGTAAGCAATAAGAAGTTCAAAACCAAGACCGAAAACTGAAAGCGGAACAGTACCGTAAAGACg                                                                                                                                              >  2:439476/1‑301 (MQ=255)
                                                                                                                                                                aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCATAGCGACGCAGaa                                                                                                                                         >  2:405697/1‑300 (MQ=255)
                                                                                                                                                                                   tcgcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAAtc                                                                                                                      <  2:439346/300‑1 (MQ=255)
                                                                                                                                                                                    tgcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAAtct                                                                                                                     <  2:47261/300‑1 (MQ=255)
                                                                                                                                                                                             gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAGAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTc                                                                                                            <  1:119971/301‑1 (MQ=255)
                                                                                                                                                                                                  aaTGTATTATTTCGTAATAGGAGGTTCGTTGTGGAGGGTTTGCAAATAGATTTCTTGACTTAACAAACAGAGAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGcc                                                                                                       <  1:402150/301‑1 (MQ=255)
                                                                                                                                                                                                             ccGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTCACAAACCGATAACTTGACAGCAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCCCTTGAAGACGCCCAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCCAAACCAAGACCGCAACCCGAATGCGGAACACTACACTAGCACCCCACAGCGCGATACCACCAATAATCTTCATTAtcaacgcccctatcacacc                                                                                            >  1:508108/1‑283 (MQ=255)
                                                                                                                                                                                                               gTAACAGGATGATCGTTATGGCTGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTTACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATa                                                                                          <  1:219090/301‑1 (MQ=255)
                                                                                                                                                                                                                 aaCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATAGCGATACCACCCGTAATCTTCTTTATTCAGGCCCAATTCCAGAATAca                                                                                        >  2:345365/1‑300 (MQ=255)
                                                                                                                                                                                                                  aCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGAACGAAACCTGAAAGCGGAACAGTACCGTAGCGACGCAGAACGCGAAACCAACAGAAATCTTCTTTATTCAGGCCCATTTCCAACAAACGc                                                                                       >  2:148372/1‑301 (MQ=255)
                                                                                                                                                                                                                   cttttccgattgtATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCt                                                                                      <  1:78536/289‑1 (MQ=255)
                                                                                                                                                                                                                   aaGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCt                                                                                      <  2:212498/300‑1 (MQ=255)
                                                                                                                                                                                                                       ttgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTc                                                                                  <  2:443827/300‑1 (MQ=255)
                                                                                                                                                                                                                               tATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCtagaacggaaaaa                                                                          >  2:78626/1‑289 (MQ=255)
                                                                                                                                                                                                                                    cGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACAACCAGTAATCTTCTTTATTCAGGCCCATTTCCACCATACGCTCGTAAACCACGTCCAGa                                                                     >  1:413120/1‑301 (MQ=255)
                                                                                                                                                                                                                                       ggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGt                                                                  <  2:213655/301‑1 (MQ=255)
                                                                                                                                                                                                                                            cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGACCGTCCAGCACGTCCAAACGttcttt                                                             >  1:419929/1‑300 (MQ=255)
                                                                                                                                                                                                                                                      gattGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAg                                                   <  1:405134/299‑1 (MQ=255)
                                                                                                                                                                                                                                                      gCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAg                                                   >  2:79168/1‑301 (MQ=255)
                                                                                                                                                                                                                                                         tGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGACCAGACGTTCTTCACGATGGGAGcca                                                >  2:412627/1‑301 (MQ=255)
                                                                                                                                                                                                                                                               aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAGCCACCAATg                                          <  2:465517/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAGCCACCAATGATCTCAcc                                  <  2:516591/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACATAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAGCCACCAATGATCTCAcc                                  >  2:147435/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                           aCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCCCGTCCAGCACGCCCACACGTTCCTCACCCTGGAAACCACCAATCACCCCCCCGAtc                              >  1:11418/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                      aaGATAAAATAAGAGTCAGCGGGCGCTGGCTCTTGAAGACGCTTAATTTGAAGCTTGCGTTACTCGGAGTACGTGGGAACGGAATCACATCTCGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTAGTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTGCACGCTGGGAGCCACCAATGATCTCACCGATGCCCGGAGCCAg                   <  1:412064/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            aaaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAGCCACCAATGATCTCACCGATGCCCGGAGCAAGAACGTc             >  2:529906/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                aaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAGCCACCAATGATCTCACCGATGCCCGGAGCCAGAACGTCCATa         >  2:495693/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                      cAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTTGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCGGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAGCCACCAATGATCTCACCGATGCCCGGAGCCAGAACGTCCATAGCCGCa   <  1:495046/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                       aGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCCCGTCCAGACGTTCTTCACGCTGGAAGCCACCAATAATCACACCGAAGCCCGGAGCCAGAACATCCAtccccccca  >  1:179349/1‑293 (MQ=255)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                       
CTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTTCACGCTGGGAGCCACCAATGATCTCACCGATGCCCGGAGCCAGAACGTCCATAGCCGCAA  >  NC_000913/987261‑987856

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: