Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 640,502 C→A L250L (CTC→CTA ahpF → alkyl hydroperoxide reductase, AhpF component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913640,5020CA98.5% 203.9 / ‑5.5 65L250L (CTC→CTAahpFalkyl hydroperoxide reductase, AhpF component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (33/31);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

TGGGCGTTCCGGCAGTGTTCGTAAACGGGAAAGAGTTTGGTCAGGGCCGCATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCGGTGCAAAATGGCGCAACATGAACGTTCCG  >  NC_000913/640258‑640748
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        
tGGGCGTTCCGGCAGTGTTCGTAAACGGGAAAGAGTTTGGTCAGGGCCGCATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCCATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATAcc                                                                                                                                                                                                                                                    >  1:224317/1‑251 (MQ=255)
 gggCGTTCCGGCAGTGTTCGTAAACGGGAAAGAGTTTGGTCAGGGCCGCATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAATGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATAcc                                                                                                                                                                                                                                                    <  1:261824/250‑1 (MQ=255)
          ggCAGTGTTCGTAAACGGGAAAGAGTTTGGTCAGGGCCGCATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTTCTCCGCACGTATAGGCCTCCGTACCGGTCTGATGGGCGAAAGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCg                                                                                                                                                                                                                                         >  1:168897/1‑251 (MQ=255)
                        aCGGGAAAGAGTTTGGTCAGGGCCGCATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCCATTTACTCCCCACGGAAAAGCCTTCCTACCGGGCTGATGGGGGAACGTTTTGGTGGTCAGATCCTAGATACCGtttgatcggacgagggtcggtgt                                                                                                                                                                                                                            >  2:177614/1‑230 (MQ=255)
                        aCGGGAAAGAGTTTGGTCAGGGCCGCATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGtt                                                                                                                                                                                                                                                 <  1:177614/230‑1 (MQ=255)
                           ggAAAGAGTTTGGTCAGGGCCGCATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTa                                                                                                                                                                                                                        <  1:228541/251‑1 (MQ=255)
                                gaggttGGTCAGGGCCGCATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCAGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTGCCGaa                                                                                                                                                                                                                   <  2:83455/247‑1 (MQ=255)
                                  gTTTGGTCAGGGCCGCATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGa                                                                                                                                                                                                                 <  2:35406/251‑1 (MQ=255)
                                         cAGGGCCGCATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCAATTTACTCCGCACGTAAAAGCATCCGTACCGGTCTGATGGGCGAAAGTTTTTGTGGTCAGATCCTAGATACCGTTGATATCGAAAAATACATTTCTGTACCGAAGACTGAAgg                                                                                                                                                                                                          >  1:91857/1‑251 (MQ=255)
                                          aGGGCCGCATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCCATTTACTCCGCACGTAAAGGCATCCGGACCCGGCTGGTGGGCGAACGTTTTGGTGGGCATATCCCAGATACCGGTGCTATCGGAAAATACATTTCTGTACCGaaacatggaggg                                                                                                                                                                                                         >  2:324274/1‑241 (MQ=255)
                                                     aCGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGc                                                                                                                                                                                              <  1:119338/251‑1 (MQ=255)
                                                                                    cTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGGAGCAATTTACTCCGCAAGTAAAGGCATCCGGACCGGTCTGATGGGCGAAAGGTTTGGGGGTCAGATCCCAAATACCGGTGATATCGAAAACTCCATTTCTTTACCGAAGACTTAAGGGCAGAAGCTGGCAGGCGCCCTGAAAGGTCACGTtggtgacaac                                                                                                                                                               >  2:298910/1‑247 (MQ=255)
                                                                                     tGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGCTGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGCCGCGGGAGACATTTACTCCGCACGTAAAAGCATCCGTACCGGGCTCATGGGCGCACGCTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAAATACATTTCTGTACCGAAGACTGAAGGGGCGAAGATGGCAGGCGCACTGAAAGTTTACGTTGATCAAtta                                                                                                                                                               >  2:237893/1‑248 (MQ=255)
                                                                                           ggAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTg                                                                                                                                                        >  1:100827/1‑251 (MQ=255)
                                                                                               aaaCGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATTGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCAATTTTCTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGCtgt                                                                                                                                                    >  2:126615/1‑251 (MQ=255)
                                                                                               aaaCGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGtt                                                                                                                                                         <  2:314495/246‑1 (MQ=255)
                                                                                               aaaCGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGtt                                                                                                                                                         >  1:314495/1‑246 (MQ=255)
                                                                                                  cGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGCCGTATTAATCGTCGGTTACGGCCCGGCGGCACCAGGGGCAGCAATTTCCTCCGACCGTAAAGGCACCCGTCCCGGTCTGATGGGCGAACGATTTGGTGGTCGGATTCTAGATACCGTTGATATCGAAAACTACAtt                                                                                                                                                                                                                              >  1:402793/1‑174 (MQ=255)
                                                                                                  cGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGGCTGATGGGCGAACGTTTTGGTGGTCAGAGCCTAGATACCGTTTATATCGAAAACTACAtt                                                                                                                                                                                                                              <  2:402793/174‑1 (MQ=255)
                                                                                                            aagaGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCCATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTTATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACCGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACCGCCa                                                                                                                                        >  1:132778/1‑250 (MQ=255)
                                                                                                               agCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAgc                                                                                                                                    <  2:160592/251‑1 (MQ=255)
                                                                                                               agCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACGGAAAGTTGACGTGGAGTAATACGACGTTGATGTGATCGACAGCCAGAgc                                                                                                                                    <  1:237893/251‑1 (MQ=255)
                                                                                                                   gAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGcc                                                                                                                                 <  2:366821/250‑1 (MQ=255)
                                                                                                                     aCAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCCAACGTTTTGGGGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTAc                                                                                                                                                                                                                       >  2:349486/1‑162 (MQ=255)
                                                                                                                     aCAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCAAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTAc                                                                                                                                                                                                                       <  1:349486/162‑1 (MQ=255)
                                                                                                                      cAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGCAAAGGCAGCAATTTACTCCGCACGTAAAGGCATCCGGACCGGTCTGATGGGCGAACGATTTGGTGGTCAGAACCTAGATACCGTTGATATCGAAAAATACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCAAGATGATGAATAAGACGTTGATGTGATCGAAAGCCAGAGCGCCAGc                                                                                                                              >  1:398783/1‑250 (MQ=255)
                                                                                                                       aaGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCaa                                                                                                                           <  2:89205/251‑1 (MQ=255)
                                                                                                                       aaGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCCTACCGGTCTGATGGGCGAACGTTTTGGGGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAAAACGAAGGTGATGGGATCGACAGCCAGAGCGCCAGCa                                                                                                                            >  2:181986/1‑250 (MQ=255)
                                                                                                                        aGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGGAGCCGGAGCGATTTACTCCCCGCGTAAAAGCATTCGTCCCGGTCTGCTGGGGGAACGCTTTGGTGGGCAGATCCCAGGTACCGTTGCTATCCGAAAATACATTTCTGTACCGAAGACTGCAAGGCAGAAGCCGGCAGGGGCCCTGACAGTTCACGTTGATGAATACCACGTTGCTTTGAGCGAAAACCTAAGCGCCAGCaa                                                                                                                           >  2:331685/1‑250 (MQ=255)
                                                                                                                                        cGTATTAAGCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCcagcagcag                                                                                                          <  1:77244/251‑1 (MQ=255)
                                                                                                                                           aTTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTAATCCGCCCGGAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCCCTGAAAGTTCACGTTGATGAATACGACGTTTATGGGATCGACAGCCAGAGCGCCCGCAAACTGATCCcagcagccggtg                                                                                                       >  1:329664/1‑246 (MQ=255)
                                                                                                                                             tAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAAGGATCCGTACCGGGCTGGTGGGCGAACGGTTTGGGGGGCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGGAGACTGAAGGGCAGAAGCTGGCAGGCGCAATGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACCTATCCCAGCAGCAGTTGa                                                                                                      >  2:78554/1‑250 (MQ=255)
                                                                                                                                              aaTCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGaa                                                                                                     <  2:10956/250‑1 (MQ=255)
                                                                                                                                                   tcgGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCa                                                                                                                                        <  1:235472/211‑1 (MQ=255)
                                                                                                                                                         ccggcccggcGGGTGCAGCGGCCGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGGGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTTAAAGTTCACGTTGATGAATACGACGTTTATGTGATCGACAGCCAGAGCGCCCGCCAAACGATCCCCGCAGCAGTTGAAGGGGGTTTGCAt                                                                                         >  2:202899/1‑251 (MQ=255)
                                                                                                                                                         ccggcccggcGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCAt                                                                                         <  1:117715/251‑1 (MQ=255)
                                                                                                                                                              ccggcGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGAt                                                                                    <  2:100827/251‑1 (MQ=255)
                                                                                                                                                                ggcggGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTg                                                                                  <  2:16714/251‑1 (MQ=255)
                                                                                                                                                                   gggTGCAGCGGCAGCAATTTACTCCGCACGTAACGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGaat                                                                               >  2:380059/1‑250 (MQ=255)
                                                                                                                                                                    ggTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATTGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAAc                                                                              >  1:301238/1‑251 (MQ=255)
                                                                                                                                                                           cGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCt                                                                       <  1:351110/251‑1 (MQ=255)
                                                                                                                                                                              cagcaATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGgc                                                                    >  1:199318/1‑251 (MQ=255)
                                                                                                                                                                                  aaTTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCTAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTATGGCGCgg                                                                >  2:51078/1‑251 (MQ=255)
                                                                                                                                                                                     ttACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTAc                                                             >  1:108331/1‑251 (MQ=255)
                                                                                                                                                                                          ccGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATc                                                                                                                    >  2:240754/1‑191 (MQ=255)
                                                                                                                                                                                          ccGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATc                                                                                                                    <  1:240754/191‑1 (MQ=255)
                                                                                                                                                                                               cGTAAAGGCATCCGTACCGGTCTGATGGGCGCACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACg                                                   >  2:173301/1‑251 (MQ=255)
                                                                                                                                                                                               cGTAAAGGCATCCGTACCGGTCCGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGAGCCCAGCAGCAGTTGAAGGGGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACg                                                   <  2:354230/251‑1 (MQ=255)
                                                                                                                                                                                                          ccGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATc                                         <  2:389410/250‑1 (MQ=255)
                                                                                                                                                                                                              aCCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTgg                                     <  1:187272/250‑1 (MQ=255)
                                                                                                                                                                                                              aCCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGc                                    <  2:94689/251‑1 (MQ=255)
                                                                                                                                                                                                              aCCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGc                                    <  1:72558/251‑1 (MQ=255)
                                                                                                                                                                                                                  gTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGAcc                                <  2:91857/251‑1 (MQ=255)
                                                                                                                                                                                                                    cTGATGGGCGAACGTTTTGGTGGTCAGATACTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCGCGCAGCATTATCGTGGCGACCgg                              >  1:116315/1‑251 (MQ=255)
                                                                                                                                                                                                                      gATGGGCGAACGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCgg                              <  1:244309/249‑1 (MQ=255)
                                                                                                                                                                                                                       aTGGGCGACCGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCGGTGc                           >  2:48890/1‑251 (MQ=255)
                                                                                                                                                                                                                                cGTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCGGTGCAAAATGgcg                  >  1:410949/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACCTTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCGGTGCAAAATGgcgc                 >  1:34622/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gTTTTGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCGGTGCAAAATGgcgc                 >  2:217637/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tttGGTGGTCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCGGTGCAAAATGGCGCa                >  1:309805/1‑250 (MQ=255)
                                                                                                                                                                                                                                          gtCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCGGTGCAAAATGGCGCAACATGAAc        <  1:202899/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gtCAGATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGACTGAAACAGCTTCTGGCGCGGTACGGAAAGCACGCAGCATTATCGTGGCGACCGGTGCAAAGTGGCGCAACATGAAc        <  2:199318/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gATCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCGGTGCAAAATGGCGCAACATGAACGTTc    <  1:169139/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aTCCTAGATACCGTTGATATCGCAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCGGTGCAAAATGGCGCAACATGAAAGTTcc   >  2:16802/1‑251 (MQ=255)
                                                                                                                                                                                                                                                tCCTAGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCGGTGCAAACTGGCGCAACATGAACGTTCCg  >  2:307724/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        
TGGGCGTTCCGGCAGTGTTCGTAAACGGGAAAGAGTTTGGTCAGGGCCGCATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCGGTGCAAAATGGCGCAACATGAACGTTCCG  >  NC_000913/640258‑640748

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