Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,863,773 C→T L61F (CTT→TTT)  ygbL → putative 3‑oxo‑tetronate 4‑phosphate decarboxylase YgbL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,863,7730CT96.9% 200.6 / ‑1.9 65L61F (CTT→TTT) ygbLputative 3‑oxo‑tetronate 4‑phosphate decarboxylase YgbL
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (32/30);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

ATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCATGAGCGATTTCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCG  >  NC_000913/2863536‑2864019
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
aTCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCATGAGCGATTTCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTcgc                                                                                                                                                                                                                                           >  1:419502/1‑251 (MQ=255)
    ggCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCATGAGCGATTTCGCAAAAGTAGAGCAGTCTTGGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTcgcc                                                                                                                                                                                                                                          <  2:120093/248‑1 (MQ=255)
        aCTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCATGAGCGATTTCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATg                                                                                                                                                                                                                                   <  1:466451/251‑1 (MQ=255)
           tCGGTGATGGCCACTTTTTTTCGCGAGCCCAAATAGAGTTTTTATCATGAGCGATTTCGCAAAAGTAGAGCAGTCTTGGCGAGAGGTGATGACGCGGATTGCCAGTTCATTCCTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCGGACGGGAATGTACTGGCGACACCGACAGGTTCATGCCTGGGCTATCACGATCCGCAGCGTTTTGCCAAAGTCGCCGCGGATGGc                                                                                                                                                                                                                                 <  1:311951/250‑1 (MQ=255)
            cGGTGATGACGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCATGAGCGATTTCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCg                                                                                                                                                                                                                                >  1:455509/1‑250 (MQ=255)
             ggTGATGACGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCATGAGCGATTTCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGa                                                                                                                                                                                                                               >  2:72309/1‑250 (MQ=255)
                 atgaCGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCATGAGCGATTTCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTACAAAGTCGCCGCGGATGGCGAATGGt                                                                                                                                                                                                                          >  1:364722/1‑251 (MQ=255)
                    aCGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCATGAGCGATTTCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTaa                                                                                                                                                                                                                       <  1:89867/251‑1 (MQ=255)
                           ttttgCACGAGTCCAAAGAGAGTGTTTATCATGAGCGATTTCGCAAAAGTAGAGCATTCTTGGCGAGAGGGGATGACGCGGATTGCCAGTTCCTTCTTTCAGCGCGGCGAGGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCGGTCGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCTATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGCTGAATGGTTAAgtggtg                                                                                                                                                                                                                 <  2:554916/250‑1 (MQ=255)
                                             agagTTTTTATCATGAGCGATTTCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAggg                                                                                                                                                                                              >  2:449266/1‑250 (MQ=255)
                                             agagTTTTTATCATGAGCGATTTCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGAGGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGt                                                                                                                                                                                              <  2:258832/251‑1 (MQ=255)
                                                    ttATCATGAGCGATTTCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTTGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACATGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCtt                                                                                                                                                                                        <  2:206418/250‑1 (MQ=255)
                                                        cATGAGCGATTTCGCAAAAGTAGAGCAGTCTTTGCGCGAGGCGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGATTGGTTAAGTGGTGGCAAACCCTCGAAAGAGGTGCTCTTTCATc                                                                                                                                                                                   >  2:284404/1‑251 (MQ=255)
                                                           gAGCGATTTCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTgg                                                                                                                                                                                <  2:28695/251‑1 (MQ=255)
                                                                aTTTCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTg                                                                                                                                                                           >  1:317557/1‑251 (MQ=255)
                                                                   tCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTa                                                                                                                                                                         <  1:72309/250‑1 (MQ=255)
                                                                     gCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTCCTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACCAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTCTCg                                                                                                                                                                      >  2:439634/1‑251 (MQ=255)
                                                                     gCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCg                                                                                                                                                                      >  1:94762/1‑251 (MQ=255)
                                                                          aGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGAGCGGCTATGCAACCGGTTCGGCTGGCCATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGCGGTGACAAACCCTCGACAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaac                                                                                                                                                                  >  2:254308/1‑250 (MQ=255)
                                                                                 cAGTCTTAGAGAGAGGAGATGCGGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATATCGGTCCGCAGCGGTTTTCCAAAGTCCCCGCGGATGGCGAATGGGTAAGTGGTGACCACCCCCCGCACGCAGTGCTCTGTCATCATGCGCTGTATCCCAAAAATCcgcg                                                                                                                                                          >  2:219351/1‑251 (MQ=255)
                                                                                  aGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGCTCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGCAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCcgcac                                                                                                                                                         >  2:7253/1‑249 (MQ=255)
                                                                                       ttGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTa                                                                                                                                                     <  1:274773/250‑1 (MQ=255)
                                                                                           gagaggagagGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTCCCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAAGCTCGGTCCGCAGCGGTTTTCCACAGTCGCCGCGGATGGCGAATGGTTAAGTGGGGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAACCCGCGCCGGAAAGCg                                                                                                                                                >  2:134885/1‑251 (MQ=255)
                                                                                           gagaGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGAt                                                                                                                                                                                                                                    >  1:61882/1‑167 (MQ=255)
                                                                                           gagaGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTGCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGAt                                                                                                                                                                                                                                    <  2:61882/167‑1 (MQ=255)
                                                                                               ggagATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggtgg                                                                                                                                            >  1:456488/1‑251 (MQ=255)
                                                                                                gagaggccgcgGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggtggt                                                                                                                                           <  2:307303/244‑1 (MQ=255)
                                                                                                  gaTGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGc                                                                                                                                         >  1:109111/1‑251 (MQ=255)
                                                                                                     gACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCAtt                                                                                                                                      <  1:410807/251‑1 (MQ=255)
                                                                                                            aTTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCGGGCAATCTGTCGCTGATTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgcacag                                                                                                                               <  2:455509/251‑1 (MQ=255)
                                                                                                            aTTGCCAGTTCATTATTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAATCGCCGCGGATGGCGAATTGTTAAGTGGTGAACACCCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGAAAAGCGGTGGGGAATTTgcac                                                                                                                                 >  2:534247/1‑249 (MQ=255)
                                                                                                                 cAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCTCTGCTTTTACCTGACGGGAATGTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACAt                                                                                                                          <  2:529500/251‑1 (MQ=255)
                                                                                                                   gTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATgg                                                                                                                        <  1:346513/251‑1 (MQ=255)
                                                                                                                       aTTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACCGCCCATGGTCgc                                                                                                                    >  2:225964/1‑250 (MQ=255)
                                                                                                                              cAGCGCGACTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGAGCTCTTTCATCTGGCGCTGTATCGCAACACTCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGAGCt                                                                                                             >  2:254373/1‑251 (MQ=255)
                                                                                                                                           cAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCa                                                                                                 >  1:5129/1‑250 (MQ=255)
                                                                                                                                            aaCCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGACTTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCAGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCCCATGGTCGACGGCGCTTTCCTGCCGGCa                                                                                                 >  2:344834/1‑249 (MQ=255)
                                                                                                                                            aaCCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAg                                                                                               <  2:243987/251‑1 (MQ=255)
                                                                                                                                                ggTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCt                                                                                           <  1:40305/251‑1 (MQ=255)
                                                                                                                                                    cGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGAc                                                                                       <  1:482199/251‑1 (MQ=255)
                                                                                                                                                         ggCAATCTGTCGCTGCTTTTACCTGAAGGGAATTTACTGGCGACCCCGAAAGGTTCATGCATGGGCAATCTCGATCCGCAGCGGTCTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCCTCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTGTCCTGCCTGCAAGGGCTGGAcagcag                                                                                  <  1:498678/251‑1 (MQ=255)
                                                                                                                                                               cTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCTATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACCGCAGCAACGt                                                                            >  1:351059/1‑251 (MQ=255)
                                                                                                                                                                 gTCGCTGCTTTTACCTGACGGGAATTTACTGGCGCCACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTGCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGtt                                                                           <  2:425210/250‑1 (MQ=255)
                                                                                                                                                                      tgctTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGCCAGCAGCAACGTTATTCGt                                                                     >  2:464375/1‑251 (MQ=255)
                                                                                                                                                                            tACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATACCTGGGCAATCTCGATCCGCAGCGGATTTCCAAAGTCGCCGCGGATGGCGAATGGATAAGTGGTGACAAACCCTCGAAAGAGATGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGAACAGCACATGGTCTACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCACAGTTATTCGTCCGtt                                                                >  2:441849/1‑250 (MQ=255)
                                                                                                                                                                              ccTGACGGGAATTTACTGGCGACACCGACCGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGCAAGAGGTGCACTTTCATCTGGCGCTGTATCGCAACAATCCGAGCTGAAACGCGGGGGTGCATTTTCACAGCACATGGTTGACGGCGATTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTCATTCGTCCGTTcac                                                             >  2:516706/1‑251 (MQ=255)
                                                                                                                                                                              ccTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTcac                                                             <  1:492992/251‑1 (MQ=255)
                                                                                                                                                                                   cGGGAATTTACTGGGGACACCGACTGGTGCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGGGACAAACCCTCTAAAGAGTTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTTGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGGTCACACCAt                                                        <  2:346353/251‑1 (MQ=255)
                                                                                                                                                                                      gAATTTCCTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACg                                                     <  1:340027/251‑1 (MQ=255)
                                                                                                                                                                                        aTTTACTGGCGACCCCGACAGGTTCATGCCTGGGCAATCTCTATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCAGCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACgtg                                                   <  2:337662/251‑1 (MQ=255)
                                                                                                                                                                                         tttACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACgtgg                                                  >  1:350540/1‑251 (MQ=255)
                                                                                                                                                                                                gCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCg                                           <  1:376966/251‑1 (MQ=255)
                                                                                                                                                                                                  gACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGGTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGAACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGATGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCggg                                         >  2:565005/1‑250 (MQ=255)
                                                                                                                                                                                                    cacCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCCAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATg                                       >  2:571399/1‑251 (MQ=255)
                                                                                                                                                                                                           aGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGGGATGCGGATGGGAAATg                                >  2:422467/1‑251 (MQ=255)
                                                                                                                                                                                                            ggTTCATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGAGGATGGGAAATGt                               >  2:576874/1‑251 (MQ=255)
                                                                                                                                                                                                             gTTCATGCCTGGGCAATCTCGATCCGCAGCGTTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTc                              <  2:385982/251‑1 (MQ=255)
                                                                                                                                                                                                                cATGCCTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCg                           >  1:492334/1‑251 (MQ=255)
                                                                                                                                                                                                                    ccTGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTgg                       >  1:17962/1‑251 (MQ=255)
                                                                                                                                                                                                                      tGGGCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAAGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTg                     >  2:89742/1‑251 (MQ=255)
                                                                                                                                                                                                                       gggCAATCGCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGc                    <  2:317557/251‑1 (MQ=255)
                                                                                                                                                                                                                        ggCAATCTCGATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGcc                   <  2:94762/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gATCCGCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTaccg          >  2:439675/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGc   <  1:422467/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGAGACAAACCCTCGAAAGAGGAGCTCTTACATCTGGCGCTGTATCGCAAAAATCCGCGCTGTAAAGCGGTGGGGCCTTTGCACAGCACATGGACGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTACCACCATACGTGGAGATGCGGTTGGGAAATTTCCCGCCGGCGCCTTATTTCCGACCGGGCg  >  2:1980/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
ATCTGGCAACTTCGGTGATGACGCCTTTTTTTCACGAGCCCAAAGAGAGTTTTTATCATGAGCGATTTCGCAAAAGTAGAGCAGTCTTTGCGAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGACAGCAGCAACGTTATTCGTCCGTTCACACCATACGTGGTGATGCGGATGGGAAATGTCCCGCTGGTGCCTTATTACCGACCGGGCG  >  NC_000913/2863536‑2864019

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