Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 292,433 G→A F172F (TTC→TTT yagK ← uncharacterized protein YagK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913292,4330GA98.2% 183.6 / ‑6.4 57F172F (TTC→TTTyagKuncharacterized protein YagK
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (28/27);  total (30/27)
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.99e-01

GAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCA  >  NC_000913/292190‑292678
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
gAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAg                                                                                                                                                                                                                                                <  2:367936/251‑1 (MQ=255)
   aCGTTCTTTAGCAGCCGCCTTGCTGATCACTATCAGTCATCTATAGACAATTGCGTAACGACACTGCATTTTAGTAATTTTTTGACACAATACAACTAAGGCCTTCTCACTGCGAAGGCCATTACAGCCTTACAGTTAGATTTGGCTACAGCCGAAGTTACGGTCTCATTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAATTCAGTAc                                                                                                                                                                                                                                             >  2:484022/1‑251 (MQ=255)
          ttAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGgt                                                                                                                                                                                                                                      >  2:574395/1‑251 (MQ=255)
             gtagCGGACTTGCTGATGACTTTCATTAATCTATAGAAAATTGCGTAACGAAACTTCATTGTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACTCTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGgtgcc                                                                                                                                                                                                                                   <  2:432004/251‑3 (MQ=255)
                             tgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTGCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTGTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTGCTGCTGAAAGTc                                                                                                                                                                                                                                                  <  2:290045/220‑1 (MQ=255)
                             tgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTCAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTc                                                                                                                                                                                                                                                  >  1:290045/1‑220 (MQ=255)
                                      gTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGCCCTTCTCACTGAGAAGGCCCTTACCGCCTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGACTCCTTCACCACACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGCACGTCAGTACTATTGGTGCCTAGCACCTACTTACAGTTTTCAgg                                                                                                                                                                                                           >  1:586447/1‑250 (MQ=255)
                                             ataGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAAt                                                                                                                                                                                                      <  2:447030/248‑1 (MQ=255)
                                             ataGAAAATTGCGTAAAGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTGCAGGAAAATg                                                                                                                                                                                                     <  2:366909/249‑1 (MQ=255)
                                               aGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTTATGTAGTTCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTTAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTCCAGGAAAATGgac                                                                                                                                                                                                  <  2:309456/250‑1 (MQ=255)
                                                  aaaTTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGgaccaga                                                                                                                                                                                              <  1:153843/251‑1 (MQ=255)
                                                                 aCTGCATTTTAGTAAGTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACGTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGtct                                                                                                                                                                               <  2:434572/251‑1 (MQ=255)
                                                                     actggtAGTAATTGTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGTACATTACCGAATTACAGTTCGATTTGGCTACAGCAGAAGTTACGGTCTCCTTCCCCAAACACTTTGCTTTCAACTTGGGTCAGGTATTCCAGTCGGGTTAACAGCGCCGGGTAGTACTGCTGAAAGTCAGTCCTCTTGGTGTCTAGCACAGTCTTACCGTTTTACTGAAAATGGaccagaccagcttgatcttgtctct                                                                                                                                                                           >  2:172578/6‑236 (MQ=255)
                                                                           aGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGtc                                                                                                                                                                     >  1:251612/1‑251 (MQ=255)
                                                                            gTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGtct                                                                                                                                                                    >  1:320837/1‑251 (MQ=255)
                                                                                 tttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTCGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAggg                                                                                                                                                               >  2:417090/1‑251 (MQ=255)
                                                                                      gAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTg                                                                                                                                                          >  1:379424/1‑251 (MQ=255)
                                                                                        aaCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGt                                                                                                                                                         >  1:54910/1‑250 (MQ=255)
                                                                                         aCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGCTTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTAc                                                                                                                                                       >  2:343481/1‑251 (MQ=255)
                                                                                           aaTACAAATAAGTCCTTCTCACTGAGAAGGCCATTACCGACTTCCAGTTCGAATTGGCTACAGCCGAAGTTACGGTCTACTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTGAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTTGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCGGGGTGATCGTCTCGCTCCAGTCTCAGGGTGCTGTACCa                                                                                                                                                     <  2:404509/251‑1 (MQ=255)
                                                                                           aaTACAAATAAGGCCTTCTCACGGAGAAGGCCATTACCGTCTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGTGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCa                                                                                                                                                     <  2:404550/251‑1 (MQ=255)
                                                                                            aTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTAcc                                                                                                                                                      >  2:411660/1‑249 (MQ=255)
                                                                                             tACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGACGCTGTACCAgg                                                                                                                                                   >  2:123625/1‑251 (MQ=255)
                                                                                                 aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTcc                                                                                                                                               <  1:411660/251‑1 (MQ=255)
                                                                                                    aaGGCCTTCTCACTGAGAAGGACATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCtgt                                                                                                                                            >  2:107512/1‑251 (MQ=255)
                                                                                                     aGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCtgtg                                                                                                                                           >  1:257463/1‑251 (MQ=255)
                                                                                                     aGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGGGATCGTCTCTCTCCCGTCTCAGGGCGCTGTACCAGGCTCCtgtg                                                                                                                                           >  2:399974/1‑251 (MQ=255)
                                                                                                                    gaAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAgg                                                                                                                            >  2:333903/1‑251 (MQ=255)
                                                                                                                      aGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGATCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTACGGtt                                                                                                                          >  1:235565/1‑251 (MQ=255)
                                                                                                                        gCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTAGTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTAt                                                                                                                        <  2:136972/251‑1 (MQ=255)
                                                                                                                            ttACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTc                                                                                                                    >  1:218283/1‑251 (MQ=255)
                                                                                                                             tACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTcc                                                                                                                   >  1:140822/1‑251 (MQ=255)
                                                                                                                             tACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTc                                                                                                                    <  2:169493/250‑1 (MQ=255)
                                                                                                                                        aGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATc                                                                                                        >  1:427578/1‑251 (MQ=255)
                                                                                                                                        aGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGGCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGGCCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATc                                                                                                        <  1:70626/251‑1 (MQ=255)
                                                                                                                                                   gCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAggtggt                                                                                             <  1:343481/251‑1 (MQ=255)
                                                                                                                                                               gTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCtt                                                                                 >  1:69862/1‑251 (MQ=255)
                                                                                                                                                                     gTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGGGGTGGTAGTAAGCGTCTttgttg                                                                            >  2:126702/1‑250 (MQ=255)
                                                                                                                                                                           ttCACCACACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaacaac                                                                      >  1:378849/1‑250 (MQ=255)
                                                                                                                                                                            tCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaacaac                                                                      <  1:482913/249‑1 (MQ=255)
                                                                                                                                                                              aCCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGc                                                                  <  1:107512/251‑1 (MQ=255)
                                                                                                                                                                                 aaaCACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGCTCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGAAGTAAGCGTCTTTGTTGAACAACAGGCaa                                                                >  2:209721/1‑250 (MQ=255)
                                                                                                                                                                                    cacTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAAtgt                                                            <  1:333903/251‑1 (MQ=255)
                                                                                                                                                                                       tttGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAgg                                                                   <  2:99044/241‑1 (MQ=255)
                                                                                                                                                                                       tttGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAgg                                                                   >  1:99044/1‑241 (MQ=255)
                                                                                                                                                                                        ttGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGAt                                                         <  2:54910/250‑1 (MQ=255)
                                                                                                                                                                                         tGCTTTCAACTTTGGTCAGGAAGTCCAGTCGGGATAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGCCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTAAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCaactgggat                                                         >  2:379337/1‑242 (MQ=255)
                                                                                                                                                                                               cAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGGGATCGTCTCTCTCCAGTCTCAGGGCTCTGTACCAGGCTCCTGAGATCATCCCTCTAAGGTTATCGTTCTGATCGTAATCCCCGAGGAGGTAGTAAGCGTCTTTGTTGAACAACAGCCAAAGGGGATAATGGCAt                                                 >  1:444696/1‑251 (MQ=255)
                                                                                                                                                                                                     tGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCttt                                                                                <  2:318504/214‑1 (MQ=255)
                                                                                                                                                                                                     tGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCttt                                                                                >  1:318504/1‑214 (MQ=255)
                                                                                                                                                                                                           ggTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGGTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTTATAGTGGCATTTACCTGATTCa                                     <  1:370562/251‑1 (MQ=255)
                                                                                                                                                                                                           ggTAGTCCAGTCGGGGTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACGATTGGTGGCGAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTGTCTCTCCAGTCGCAGGGCGCTGTGCCAGGCTCCTGGGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTTGTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAtgca                                     <  2:384186/251‑4 (MQ=255)
                                                                                                                                                                                                              aGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGaa                                  <  2:462976/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gcCTGGTAGTTCTGCTGAAAGGTAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCGCTCCAGTCTCAGGGCGCGGTTCCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTTTTTGTTGAACAACAGGCAAAGGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGTCCATataat              <  2:586347/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ggTAGTTCTGCTGAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATaataaat          <  1:18037/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gctgAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAggg    <  1:574395/248‑1 (MQ=255)
                                                                                                                                                                                                                                               gctgAAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCa  >  1:121841/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
GAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCA  >  NC_000913/292190‑292678

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