Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,629,498 T→G 5.9% L472F (TTA→TTC ydfI ← putative NAD‑dependent D‑mannonate oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,629,4980TG 5.9% 1832.8 / 6.6 589L472F (TTA→TTCydfIputative NAD‑dependent D‑mannonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (188/344);  new base G (50/2);  total (243/346)
Fisher's exact test for biased strand distribution p-value = 6.19e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 621 total aligned reads displayed.

AACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACTTTGTACAGCCTTCTGAATAACAGGTAA  >  NC_000913/1629349‑1629647
                                                                                                                                                     |                                                                                                                                                     
 aCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCAAGCGCta                                                                                                                                                      >  1:1400052/1‑150 (MQ=255)
  ccGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaa                                                                                                                                                     >  2:1274467/1‑150 (MQ=255)
   cGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaat                                                                                                                                                    <  2:886632/150‑1 (MQ=255)
      aGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATaaa                                                                                                                                                 <  1:6091964/150‑1 (MQ=255)
      aGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATaaa                                                                                                                                                 <  2:5871840/150‑1 (MQ=255)
      aGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATaaa                                                                                                                                                 <  2:5328940/150‑1 (MQ=255)
       gttgttGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGCATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAg                                                                                                                                                <  2:1114416/150‑1 (MQ=255)
       gttgttGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAg                                                                                                                                                >  1:9758760/1‑150 (MQ=255)
        ttgttgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGa                                                                                                                                               <  1:3143545/150‑1 (MQ=255)
        ttgttgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGa                                                                                                                                               <  1:3327295/150‑1 (MQ=255)
         tgttgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGAc                                                                                                                                              <  2:7831918/150‑1 (MQ=255)
         tgttgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGAc                                                                                                                                              <  2:3217120/150‑1 (MQ=255)
              aaGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTa                                                                                                                                         >  2:1984583/1‑150 (MQ=255)
               aGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTAc                                                                                                                                        <  2:1782263/150‑1 (MQ=255)
                     tgAAGAGAAACTGGGACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTcc                                                                                                                                  <  2:2585809/150‑1 (MQ=255)
                     tgAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGTCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTcc                                                                                                                                  <  2:4130232/150‑1 (MQ=255)
                      gAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCg                                                                                                                                 >  2:7281661/1‑150 (MQ=255)
                       aaGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGt                                                                                                                                <  2:8632054/150‑1 (MQ=255)
                       aaGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTAGCGCCAAGCGCTAATAAAGAAAACTAAGTTTACGt                                                                                                                                >  1:8562366/1‑150 (MQ=255)
                         gagaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCa                                                                                                                              <  2:1469354/150‑1 (MQ=255)
                          agaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCCTGCGCGAATAAAGACAAGTACGTTTACGGCAc                                                                                                                             >  1:2986588/1‑150 (MQ=255)
                          agaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATAGCTTTTTTACTTCACGGAATATGTTGCCACGGTCGCTTTCGAGCCATGCGCTAATAAAGACCAGTACGTTTCCGTCAc                                                                                                                             >  2:5660468/1‑150 (MQ=255)
                            aaaCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTATTTACTTCACGGAATAATATGCCACGGGCGCGTTCGCGCCCACCGCGACTAAAGACAAGTCCCATTCCGCCAcac                                                                                                                           >  1:541053/1‑148 (MQ=255)
                              aCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCtt                                                                                                                         <  2:10290169/150‑1 (MQ=255)
                               cTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTg                                                                                                                        <  2:3348833/150‑1 (MQ=255)
                               cTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTg                                                                                                                        >  2:3718471/1‑150 (MQ=255)
                                 ggCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCa                                                                                                                      <  2:4572233/150‑1 (MQ=255)
                                  gCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAg                                                                                                                     >  2:10156580/1‑150 (MQ=255)
                                  gCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCCTCACTCTTGCAg                                                                                                                     >  2:6262783/1‑150 (MQ=255)
                                    accacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACACTTGCAGTa                                                                                                                   >  2:3521730/1‑150 (MQ=255)
                                      cacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaa                                                                                                                 <  1:9230242/150‑1 (MQ=255)
                                        cTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaaca                                                                                                               >  2:837087/1‑150 (MQ=255)
                                          ggACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaaa                                                                                                             <  2:614675/150‑1 (MQ=255)
                                           gACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaaac                                                                                                            <  1:9504190/150‑1 (MQ=255)
                                              cTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTAt                                                                                                         >  2:6295740/1‑150 (MQ=255)
                                              cTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCGTTCGCGCCATGCGCTAATAAAGACAATTAAGTTTCCGTCACTCTTGCAGTAAAAAAACTAt                                                                                                         >  2:9462298/1‑150 (MQ=255)
                                               tGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTAGCGCCATGCGCGAATAAAGACAAATAAGTTTGCGTCAAGAGTGCAGTAAAAAAACAAtt                                                                                                        >  1:4958713/1‑150 (MQ=255)
                                                 aTTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAAAAAACAATTGt                                                                                                      >  1:7144173/1‑150 (MQ=255)
                                                 aTTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTAGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGGCACTCTTGCAGTaaacaaaaaaatgt                                                                                                      >  1:6960579/1‑143 (MQ=255)
                                                  ttCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTc                                                                                                     <  1:8854318/150‑1 (MQ=255)
                                                  ttCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTc                                                                                                     <  1:4340033/150‑1 (MQ=255)
                                                   tCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCt                                                                                                    <  1:9064905/150‑1 (MQ=255)
                                                    cGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTg                                                                                                   <  1:5777110/150‑1 (MQ=255)
                                                       ctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCa                                                                                                <  1:2704211/150‑1 (MQ=255)
                                                        tcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCGTTCGCGCCCAGCGCGAATAAAAACAAATACGGTTGCGGCACTCATGCAGTAAAAAAAAAATTTTATGGCaa                                                                                               >  1:2864168/1‑150 (MQ=255)
                                                         cGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCaaa                                                                                              <  2:7488902/150‑1 (MQ=255)
                                                          ggCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAat                                                                                             <  2:4201916/150‑1 (MQ=255)
                                                            cATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatca                                                                                           <  1:5228524/150‑1 (MQ=255)
                                                            cATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACACTTGCAGTAAACCAACTATTGTCTGGCAAatcc                                                                                           >  2:2330882/1‑149 (MQ=255)
                                                                aaCTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCAc                                                                                       <  1:3487949/150‑1 (MQ=255)
                                                                aaCTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCCCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCAc                                                                                       >  2:5052345/1‑150 (MQ=255)
                                                                   tAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTAGCGCCCTGCGCTAATAAATAAAACTAAGTTTGCGTCACACGTTCAGTAAACAAACAAATTTATGGCAAATAATAAACaa                                                                                    >  1:6477788/1‑150 (MQ=255)
                                                                       aaaagcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCGGTCGCGCCAAGCGCGAATAAAGACAACTAAGTTTCCGGCCCTCTTGCCGGAAACCAACAAATGTCTGGCCCATCATCACCAAAGaa                                                                                >  1:1654928/7‑149 (MQ=255)
                                                                          gtgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGc                                                                             <  2:2634085/150‑1 (MQ=255)
                                                                          gtgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGc                                                                             <  1:8197053/150‑1 (MQ=255)
                                                                           gcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCGTTAGCGCCATGCGCTTATAAATACAAGTACGTTTCCGGCACTCATTCAGTAAACAAAAAATTTTAATGCAAAACATCAACAAAGATAGcc                                                                            >  1:10118207/3‑150 (MQ=255)
                                                                           gcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCGTTAGCGCCAAGCGCGAAAAAATACAAGTAACGTTACGTCACTCTTTCAGTTAAAAAAAAAATGTGTGGCCAAACATCACCAAAGAtaggc                                                                            >  1:7910241/3‑146 (MQ=255)
                                                                            gAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTGGCGCCATTGGCTAATAAAAAAAAGTACCTTTACGTCACGCTTTCAGTTAACAAAAAAATTTCTGGGAAATCATCACCAAAGATAGccg                                                                           >  1:390873/1‑149 (MQ=255)
                                                                             aaTGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCtt                                                                          <  2:7608489/150‑1 (MQ=255)
                                                                                   gccgCTGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGc                                                                    <  1:7511788/150‑1 (MQ=255)
                                                                                   gccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTAGCGCCATGCGCGAATAAAGAAAAGTAAGTTTCCGTCAAGCATTCAGTAAACAAACTATTTTGTGGCAAAAAATAAACAAAGATAGCCTTAATAGc                                                                    >  1:6255569/1‑150 (MQ=255)
                                                                                    ccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGcc                                                                   <  2:9859804/150‑1 (MQ=255)
                                                                                           cggcATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATg                                                            >  2:4359538/1‑150 (MQ=255)
                                                                                              cATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACt                                                         >  2:6398190/1‑150 (MQ=255)
                                                                                              cATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACt                                                         <  2:1593106/150‑1 (MQ=255)
                                                                                                 tGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGa                                                      <  1:1252190/150‑1 (MQ=255)
                                                                                                 tGCATTTTTACTTCACGGAATATATTGCCACGGTCGCTTTCGCGCCATGCGATAATAAACAAAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACAATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGa                                                      <  1:10347771/150‑1 (MQ=255)
                                                                                                   cTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGAcg                                                    >  2:2485072/1‑150 (MQ=255)
                                                                                                    ttttttACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGAcgc                                                   >  2:7137135/1‑150 (MQ=255)
                                                                                                    ttttttACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGAcgc                                                   <  2:4409992/150‑1 (MQ=255)
                                                                                                    ttttttACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGAcgc                                                   <  1:9917143/150‑1 (MQ=255)
                                                                                                    ttttttACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGAcgc                                                   >  2:7918575/1‑150 (MQ=255)
                                                                                                      ttttACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGAcgcgc                                                 >  2:7012949/1‑150 (MQ=255)
                                                                                                      ttttACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGAcgcgc                                                 <  1:4601124/150‑1 (MQ=255)
                                                                                                       tttACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGAAAAGTAAGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATAATCACCAAAGATAGCCATAATCGCCAGCAAAGACAGGAcccgcg                                                >  2:10640353/1‑150 (MQ=255)
                                                                                                           cTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCttt                                            <  2:2131479/150‑1 (MQ=255)
                                                                                                           cTTCACGGAATATTGTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCttt                                            <  2:5503433/150‑1 (MQ=255)
                                                                                                            ttCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTc                                           >  2:10806972/1‑150 (MQ=255)
                                                                                                            ttCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTc                                           >  2:9227346/1‑150 (MQ=255)
                                                                                                            ttCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTATGGCAAATCATAACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTc                                           >  2:5751338/1‑150 (MQ=255)
                                                                                                             tCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTcc                                          <  1:5496606/150‑1 (MQ=255)
                                                                                                              cACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTGCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTccc                                         >  2:9120516/1‑150 (MQ=255)
                                                                                                              cACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTAAGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTccc                                         >  2:4292411/1‑150 (MQ=255)
                                                                                                               aCGGAATATTTTGCCACGGTCGCTTTCGCGCCGTGCGCTAATAAATACAAGTAAGTTTGCGTCACTCTTTCAGTAAACAAACTATTGTCTGGAAAATCATCACCAAATATCGCCTTTATCGCCAGCAATGACTGGACGCGCGCtttagct                                        >  2:3979282/1‑146 (MQ=255)
                                                                                                                 ggAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTc                                      <  2:8503782/150‑1 (MQ=255)
                                                                                                                     tatTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCa                                  <  2:6922335/150‑1 (MQ=255)
                                                                                                                          tGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCactact                             <  2:5170159/150‑1 (MQ=255)
                                                                                                                           gCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACtt                            <  2:2777424/150‑1 (MQ=255)
                                                                                                                                 gTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACTTTGTACa                      <  2:6643143/150‑1 (MQ=255)
                                                                                                                                 gTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACTTTGTACa                      <  2:3014901/150‑1 (MQ=255)
                                                                                                                                      tttCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTGAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACGTTGTACCGCCtt                 >  2:2190140/1‑150 (MQ=255)
                                                                                                                                         cgcgCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACTTTGTACAGCCTTCTg              >  1:10686396/1‑150 (MQ=255)
                                                                                                                                          gcgcCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACTTTGTACAGCCTTCTGa             >  2:9653567/1‑150 (MQ=255)
                                                                                                                                           cgcCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACTTTGTACAGCCTTCTGaa            >  1:2408662/1‑150 (MQ=255)
                                                                                                                                            gcCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACTTTGTACAGCCTTCTGAAt           <  1:1884058/150‑1 (MQ=255)
                                                                                                                                            gcCATGCGCTAATAAAGACAAGTACGTGTCCGTCACTCTGGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCGTTCCCTTCGGCACTACTTTGTACAGCCTTCTGAAt           <  2:10069894/150‑1 (MQ=255)
                                                                                                                                              cATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACTTTGTACAGCCTTCTGAATaa         >  1:721626/1‑150 (MQ=255)
                                                                                                                                               aTGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACTTTGTACAGCCTTCTGAATAAc        <  2:6792470/150‑1 (MQ=255)
                                                                                                                                                tGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACTTTGTACAGCCTTCTGAATAACa       <  1:2254317/150‑1 (MQ=255)
                                                                                                                                                tGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACTTTGTACAGCCTTCTGAATAACa       <  1:10512952/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCAATGACTGGACGCGCGCTTTCCCTTCGGCACTACTTTGTACAGCCTTCTGAATAACAGGTAA  >  NC_000913/1629349‑1629647

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: