Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,794,640 T→A A218A (GCT→GCA)  rfaD → ADP‑L‑glycero‑D‑mannoheptose 6‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,794,6400TA97.3% 254.0 / NA 75A218A (GCT→GCA) rfaDADP‑L‑glycero‑D‑mannoheptose 6‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (38/35);  minor base C (2/0);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCTGATGTGAATCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCGGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCT  >  NC_000913/3794402‑3794874
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
tCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGAtgtaaacctg                                                                                                                                                                                                                                >  1:109625/1‑244 (MQ=255)
 cTACGGATACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCtgt                                                                                                                                                                                                                               >  1:127770/1‑251 (MQ=255)
   aCGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTTTTTGGGCGACGTGGCAGATGTAAACCTGTgg                                                                                                                                                                                                                             >  2:113690/1‑251 (MQ=255)
   aCGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGAtgtaaacctg                                                                                                                                                                                                                                <  1:83637/248‑8 (MQ=255)
      gTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTc                                                                                                                                                                                                                          >  1:400511/1‑251 (MQ=255)
             aaaaTTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGaaa                                                                                                                                                                                                                   <  2:380609/251‑1 (MQ=255)
              aaaTTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTTAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCCGATGTAAACCTGTGGTTCCTGGaaa                                                                                                                                                                                                                   >  1:142176/1‑250 (MQ=255)
                      gTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGttt                                                                                                                                                                                                          <  2:313521/251‑1 (MQ=255)
                           atgaATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGc                                                                                                                                                                                                     <  2:234911/251‑1 (MQ=255)
                               atatGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGACGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCt                                                                                                                                                                                                 >  1:295898/1‑251 (MQ=255)
                                 atGTTCGTCAAATCCTGCCGGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTGCCGGCATCTTc                                                                                                                                                                                               <  2:71251/251‑1 (MQ=255)
                                     tCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAAt                                                                                                                                                                                            <  2:253481/250‑1 (MQ=255)
                                          aaaTCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGt                                                                                                                                                                                      >  1:272539/1‑251 (MQ=255)
                                                 gcctgAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTc                                                                                                                                                                               >  1:35656/1‑251 (MQ=255)
                                                  cctgAAGCGATCTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCg                                                                                                                                                                              >  1:6890/1‑251 (MQ=255)
                                                           aaCTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATc                                                                                                                                                                     <  1:371635/251‑1 (MQ=255)
                                                             cTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATCGCGAGCGTCGCTTTCCATCTCAACAACCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACATGTGGATCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCt                                                                                                                                                                   >  2:271855/1‑251 (MQ=255)
                                                                   gATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAgg                                                                                                                                                             >  1:6719/1‑251 (MQ=255)
                                                                          ggCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGc                                                                                                                                                      <  1:330023/251‑1 (MQ=255)
                                                                          ggCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGc                                                                                                                                                      <  1:178375/251‑1 (MQ=255)
                                                                           gCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCt                                                                                                                                                     <  1:60772/251‑1 (MQ=255)
                                                                                cGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGc                                                                                                                                                >  2:373849/1‑251 (MQ=255)
                                                                                    aTTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACGTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTTGAAAATGGCGTTTCCGGCATCGTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTAc                                                                                                                                             <  1:273416/250‑1 (MQ=255)
                                                                                          aCGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTGGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCt                                                                                                                                      <  2:34230/251‑1 (MQ=255)
                                                                                            gTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCtt                                                                                                                                     >  2:261699/1‑250 (MQ=255)
                                                                                             tttATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGATAACTTCAAACGCGTCTTCGTCTATGTGGGCGCCGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTa                                                                                                                                    <  2:109625/250‑1 (MQ=255)
                                                                                              ttATGGACCGAGTGAAGGCCATAAAGGCAGTATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTAt                                                                                                                                   <  2:142176/250‑1 (MQ=255)
                                                                                                aTGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGTCTTATcac                                                                                                                                <  1:113690/251‑1 (MQ=255)
                                                                                                aTGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATcac                                                                                                                                >  2:378443/1‑251 (MQ=255)
                                                                                                   gACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCATAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAAGGGCGTTTCCGGCATCTTCAATCTCGGTACTTGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTGCGCTGGCGTATCACAAg                                                                                                                             <  1:130714/251‑1 (MQ=255)
                                                                                                           gAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGcc                                                                                                                      <  1:353819/250‑1 (MQ=255)
                                                                                                            aaGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTCGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGcc                                                                                                                      >  1:405211/1‑249 (MQ=255)
                                                                                                               gCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGcc                                                                                                                      >  1:296273/1‑246 (MQ=255)
                                                                                                                ccATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGcc                                                                                                                      >  1:167150/1‑245 (MQ=255)
                                                                                                                ccATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCg                                                                                                                <  1:256821/251‑1 (MQ=255)
                                                                                                                ccATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGCTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACACGaaacacc                                                                                                                      >  2:302839/1‑241 (MQ=255)
                                                                                                                  aTAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCCCCAACCAGCTTAACACCGGAGAATCACCGAAGCTGTTTGAAGGTAGCGAGAAATTAAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGAAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTCCTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGa                                                                                                               >  2:226983/1‑250 (MQ=255)
                                                                                                                  aTAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGACGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGCCGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGa                                                                                                               >  2:311053/1‑250 (MQ=255)
                                                                                                                   tAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAAt                                                                                                             >  2:191285/1‑251 (MQ=255)
                                                                                                                          agcaTGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTcc                                                                                                      >  2:433806/1‑251 (MQ=255)
                                                                                                                               ggCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTcc                                                                                                 >  1:4531/1‑251 (MQ=255)
                                                                                                                                 cGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCg                                                                                               >  2:318977/1‑251 (MQ=255)
                                                                                                                                     cGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGCAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGAt                                                                                            >  2:371928/1‑250 (MQ=255)
                                                                                                                                      gTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTGCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTTCGCTGGCTTATCACAAGAAAGGCCAGATCGAATGCTTTCCGTTCCGGGATaa                                                                                          <  2:205347/251‑1 (MQ=255)
                                                                                                                                        cGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAAc                                                                                        <  1:373849/251‑1 (MQ=255)
                                                                                                                                         gCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACt                                                                                       <  1:354300/251‑1 (MQ=255)
                                                                                                                                           tttCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTg                                                                                      >  1:174683/1‑250 (MQ=255)
                                                                                                                                              ccATCTCAACACCCAGTTTAACAACGGTGAATCCCCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGaaa                                                                                   <  2:108466/250‑1 (MQ=255)
                                                                                                                                              ccATCTCAACACCCAGCTTACCAACGGTGAATCGCCGAAGCTGTTTGAAAGTATCGAGAACTTCAAACGCGACTTCGTCTATGTGTGCGACGTGGCAGATGTAAACCTGGGTTTCCTGGAAAATGGCGTTACCGGCATCTTCAACCTCGGTACTGGTCGTGCGGAATCATTCCATGCTGTAGCTGTTGCTACGCTGGCTTATCACAAGAAAGGCCAGATCAAATGCTTTCCGTCCCCGGATAAACTGAAAg                                                                                  <  2:282498/251‑1 (MQ=255)
                                                                                                                                               cATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAgg                                                                                 >  2:34385/1‑251 (MQ=255)
                                                                                                                                                         cccAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTGCCTGGAAAATGGCGTTTCCGGCATGTGCAATCGCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCGACGCTGGCTTATCACAAGAAAGGGCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAg                                                                       <  1:226983/251‑1 (MQ=255)
                                                                                                                                                            aGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCt                                                                                                                                          >  1:369121/1‑181 (MQ=255)
                                                                                                                                                            aGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGGAGCTGAGGCTACGCt                                                                                                                                          <  2:369121/181‑1 (MQ=255)
                                                                                                                                                              cTTAACACCGGTGAATCACCGAAGATGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCAATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGCTCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTCTCAGGCGtt                                                                  >  2:21798/1‑251 (MQ=255)
                                                                                                                                                                      cGGTGAATCACCGCAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAACCCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAgg                                                          >  1:305212/1‑251 (MQ=255)
                                                                                                                                                                      cGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAgg                                                          >  1:6636/1‑251 (MQ=255)
                                                                                                                                                                            aTCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATc                                                    <  2:167150/251‑1 (MQ=255)
                                                                                                                                                                             tCACCGAAGCTGTTTGACGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCt                                                   >  2:261329/1‑251 (MQ=255)
                                                                                                                                                                                 cGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGAcg                                               >  1:58297/1‑250 (MQ=255)
                                                                                                                                                                                    aGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGCCTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAAt                                            >  1:145751/1‑251 (MQ=255)
                                                                                                                                                                                     gCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATc                                           <  1:415910/251‑1 (MQ=255)
                                                                                                                                                                                        gTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTgc                                        <  1:313833/251‑1 (MQ=255)
                                                                                                                                                                                          ttGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTg                                                                                                                                                                            <  2:213290/117‑1 (MQ=255)
                                                                                                                                                                                          ttGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTg                                                                                                                                                                            >  1:213290/1‑117 (MQ=255)
                                                                                                                                                                                            gAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGTGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGa                                                                                                               <  1:59408/176‑1 (MQ=255)
                                                                                                                                                                                            gAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGa                                                                                                               >  2:59408/1‑176 (MQ=255)
                                                                                                                                                                                            gAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTgcgcg                                     <  1:14614/250‑1 (MQ=255)
                                                                                                                                                                                                  aGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCggg                              <  1:261329/251‑1 (MQ=255)
                                                                                                                                                                                                             cAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACaaa                    <  1:138722/250‑1 (MQ=255)
                                                                                                                                                                                                             cAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAAc                   <  1:311053/251‑1 (MQ=255)
                                                                                                                                                                                                               aaCGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCg                 >  1:349455/1‑251 (MQ=255)
                                                                                                                                                                                                                      cTTCGTCTAAGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGATCaaa           >  1:97619/1‑250 (MQ=255)
                                                                                                                                                                                                                        tCGTCTATGTGGGCGACGTGGCCGATGTAAACCTGTGGTTCCAGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCTTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGACAGGCCGCTATCAGGCGTTCACTCAGGCCGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAAcc        >  2:384868/1‑251 (MQ=255)
                                                                                                                                                                                                                          gTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAATGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGCCAAACCGTTCAAAACCGt      >  1:354342/1‑251 (MQ=255)
                                                                                                                                                                                                                              aTGTGGGCGACGTGGTAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCt  <  2:127770/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
TCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCTGATGTGAATCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCGGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCT  >  NC_000913/3794402‑3794874

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