Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 641,355 G→A 7.1% intergenic (+37/+84) ahpF → / ← uspG alkyl hydroperoxide reductase, AhpF component/universal stress protein G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913641,3550GA 7.1% 437.5 / 14.7 141intergenic (+37/+84)ahpF/uspGalkyl hydroperoxide reductase, AhpF component/universal stress protein G
Reads supporting (aligned to +/- strand):  ref base G (73/58);  new base A (5/5);  total (78/63)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

Only 100 of 155 total aligned reads displayed.

GCGGGTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAG  >  NC_000913/641208‑641504
                                                                                                                                                   |                                                                                                                                                     
gCGGGTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCggg                                                                                                                                                     <  2:346102/150‑1 (MQ=255)
gCGGGTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCggg                                                                                                                                                     <  1:2205643/150‑1 (MQ=255)
          gTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTg                                                                                                                                           <  2:359382/150‑1 (MQ=255)
            acgacgGTTCCGTACAAGCAGATCATCATGGCCTCTGGCGATGGTGCCTAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGGATAAGAAGAAGTAAGATTCACCTGGATTTGCTTAGCCGCCGGGGTCAAACCTggc                                                                                                                                         <  2:569140/150‑1 (MQ=255)
             cgacgGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTggcg                                                                                                                                        >  1:1150378/1‑150 (MQ=255)
               acgGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCAGGTTTGACCCCggcggc                                                                                                                                      <  2:2475603/150‑1 (MQ=255)
               acgGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCAGGTTTGACCCCggcggc                                                                                                                                      <  1:109662/150‑1 (MQ=255)
               acgGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCAGGTTTGACCCCggcggc                                                                                                                                      <  2:842058/150‑1 (MQ=255)
                cgGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCt                                                                                                                                     <  2:2097413/150‑1 (MQ=255)
                  gTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCttt                                                                                                                                   >  2:1546983/1‑150 (MQ=255)
                   ttCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCAGGTTTGACCCCGGCGGCtaag                                                                                                                                  <  2:2401068/150‑4 (MQ=255)
                    tCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTATGAGTGCTTTTGAATACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCttttt                                                                                                                                 >  2:1055851/1‑150 (MQ=255)
                      cGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTa                                                                                                                               >  1:1298234/1‑150 (MQ=255)
                       gTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTAt                                                                                                                              >  1:1821636/1‑150 (MQ=255)
                        tACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATg                                                                                                                             <  2:1298234/150‑1 (MQ=255)
                          cAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGc                                                                                                                           >  1:423756/1‑150 (MQ=255)
                                gATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTaaa                                                                                                                     >  2:1876741/1‑150 (MQ=255)
                                gATCATCATCGCCACTGGCGAAGGTACCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATgg                                                                                                                            <  2:1065576/143‑1 (MQ=255)
                                   catcatcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAg                                                                                                                  >  2:1205861/1‑150 (MQ=255)
                                   catcatcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAg                                                                                                                  <  1:1289227/150‑1 (MQ=255)
                                    atcatcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGc                                                                                                                 >  1:2330506/1‑150 (MQ=255)
                                     tcatcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGcc                                                                                                                <  2:1223736/150‑1 (MQ=255)
                                          gCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCccgcc                                                                                                           <  2:1453350/150‑1 (MQ=255)
                                              cTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCTGGa                                                                                                       >  1:2389188/1‑150 (MQ=255)
                                               tGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGAt                                                                                                      <  1:1124588/150‑1 (MQ=255)
                                               tGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGAt                                                                                                      >  1:1389791/1‑150 (MQ=255)
                                                ggCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATg                                                                                                     <  2:2370500/150‑1 (MQ=255)
                                                  cGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGc                                                                                                   <  2:866167/150‑1 (MQ=255)
                                                        tGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCggg                                                                                             >  1:1524662/1‑150 (MQ=255)
                                                        tGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCggg                                                                                             >  1:465641/1‑150 (MQ=255)
                                                        tGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCggg                                                                                             >  1:1201550/1‑150 (MQ=255)
                                                           cAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAg                                                                                          >  1:909893/1‑150 (MQ=255)
                                                              aGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGc                                                                                       <  2:624414/150‑1 (MQ=255)
                                                               gCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCg                                                                                      <  2:423481/150‑1 (MQ=255)
                                                                      tGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACt                                                                               >  1:2135531/1‑150 (MQ=255)
                                                                       gAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACtt                                                                              >  1:948469/1‑150 (MQ=255)
                                                                        aGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGTAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTa                                                                             <  2:146132/150‑1 (MQ=255)
                                                                         gTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTAt                                                                            >  1:1978256/1‑150 (MQ=255)
                                                                            cTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCt                                                                         >  1:842196/1‑150 (MQ=255)
                                                                                tGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCa                                                                     >  2:1309939/1‑150 (MQ=255)
                                                                                 gACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCaa                                                                    >  1:2090186/1‑150 (MQ=255)
                                                                                   cTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAAtt                                                                  >  2:52179/1‑150 (MQ=255)
                                                                                     aCCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTaa                                                                >  2:2345058/1‑150 (MQ=255)
                                                                                     aCCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCCCGGGCAGTGCGGGATACTTATTCTGGCAATTAAc                                                               >  2:1040511/1‑150 (MQ=255)
                                                                                      ccTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAAc                                                               <  1:123761/150‑1 (MQ=255)
                                                                                       cTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACg                                                              <  1:1624642/150‑1 (MQ=255)
                                                                                        tGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCa                                                                     >  1:1008768/1‑142 (MQ=255)
                                                                                        tGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGc                                                             <  1:1447551/150‑1 (MQ=255)
                                                                                        tGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGc                                                             <  2:2123506/150‑1 (MQ=255)
                                                                                             cGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAAc                                                        >  1:1218031/1‑150 (MQ=255)
                                                                                             cGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAAc                                                        >  2:746783/1‑150 (MQ=255)
                                                                                             cGCACCAAAACTGCATAAGAAGAAGTAAGAATCACCTGCAATTGCTTAGCCGCCGCGGTCAAACCTGCCGGCTTTTATATCGCATTAAAAAGCCCCGCCGGGAAGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAAc                                                        >  1:786317/1‑150 (MQ=255)
                                                                                               cACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCa                                                      >  1:1742167/1‑150 (MQ=255)
                                                                                               aaCCAAAACTGCATAAGAAGAAGTAAGATTCACATGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATAATTATTCTGGCAATTAACGCACAACCa                                                      >  2:878604/2‑150 (MQ=255)
                                                                                                 ccAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGc                                                    >  2:2453662/1‑150 (MQ=255)
                                                                                                    aaaCTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCAcc                                                 <  1:1524164/150‑1 (MQ=255)
                                                                                                     aaCTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGAt                                                                                  <  1:1061605/116‑1 (MQ=255)
                                                                                                         gCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAg                                            >  2:615765/1‑150 (MQ=255)
                                                                                                            tAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGAtt                                         <  1:1546983/150‑1 (MQ=255)
                                                                                                              agaagaagTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtgg                                       >  2:1977679/1‑150 (MQ=255)
                                                                                                                  gaagTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtg                                   <  1:1840973/150‑1 (MQ=255)
                                                                                                                   aagTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtgg                                  >  2:1760515/1‑150 (MQ=255)
                                                                                                                    agTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggc                                 >  2:1343222/1‑150 (MQ=255)
                                                                                                                     gTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggcg                                >  1:1791234/1‑150 (MQ=255)
                                                                                                                     gTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggcg                                >  1:2382788/1‑150 (MQ=255)
                                                                                                                     gTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggcg                                >  2:1812644/1‑150 (MQ=255)
                                                                                                                     gTAAGATTCACCTGCAATTGCGTAGCCGCCGGGGGCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACGTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggcg                                <  1:1644048/150‑1 (MQ=255)
                                                                                                                       aaGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGa                              >  2:886719/1‑150 (MQ=255)
                                                                                                                         gATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATc                            >  2:7365/1‑150 (MQ=255)
                                                                                                                           ttCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCAc                          <  2:1998764/150‑1 (MQ=255)
                                                                                                                           ttCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCAc                          <  1:1670645/150‑1 (MQ=255)
                                                                                                                            tCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACg                         <  1:615765/150‑1 (MQ=255)
                                                                                                                             cACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGc                        <  2:987317/150‑1 (MQ=255)
                                                                                                                             cACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGc                        >  1:2010914/1‑150 (MQ=255)
                                                                                                                             cACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGc                        >  1:1036179/1‑150 (MQ=255)
                                                                                                                               ccTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCATGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTc                      <  2:110737/150‑1 (MQ=255)
                                                                                                                                cTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCg                     >  1:1070461/1‑150 (MQ=255)
                                                                                                                                     ataaaaaaGCCGCCAGGTTTGACCCCGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCg                >  2:1327759/8‑150 (MQ=255)
                                                                                                                                     ataaaaaaGCCGCCAGGTTTGACCCCGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCg                >  1:1456350/8‑150 (MQ=255)
                                                                                                                                       tGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGtt              <  2:1524662/150‑1 (MQ=255)
                                                                                                                                        gCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGtca             <  1:2178174/150‑3 (MQ=255)
                                                                                                                                        gCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTa             <  2:423756/150‑1 (MQ=255)
                                                                                                                                        gCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTa             >  1:1059347/1‑150 (MQ=255)
                                                                                                                                        gCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTACGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTa             <  1:1553705/150‑1 (MQ=255)
                                                                                                                                         cTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAg            >  1:2449909/1‑150 (MQ=255)
                                                                                                                                         cTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAg            >  1:1484063/1‑150 (MQ=255)
                                                                                                                                         cTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAg            >  1:342894/1‑150 (MQ=255)
                                                                                                                                          aaaGCCGCCAGGTTTGACCCCGGCGGCTTTTTTATGCATTAAAAAGCCCCGCCGGGATGGCTACGCGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGaa          >  2:1223919/3‑150 (MQ=255)
                                                                                                                                            aGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTTGCAATTAACGCACAACGAGCACCGGCAGATTGGCGTGGCGGATCACTCTCGAGGCGTTAGAAc         <  2:1069965/150‑1 (MQ=255)
                                                                                                                                            aGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAAc         <  2:465641/150‑1 (MQ=255)
                                                                                                                                            aGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAAc         <  1:1876741/150‑1 (MQ=255)
                                                                                                                                            aGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAAc         <  1:1417579/150‑1 (MQ=255)
                                                                                                                                             gccgccGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAAcc        >  1:1403532/1‑150 (MQ=255)
                                                                                                                                             gccgccGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAAcc        <  1:1387053/150‑1 (MQ=255)
                                                                                                                                              ccgccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCt       >  1:1659683/1‑150 (MQ=255)
                                                                                                                                              ccgccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCt       <  1:2453662/150‑1 (MQ=255)
                                                                                                                                               cgccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTa      <  2:1536527/150‑1 (MQ=255)
                                                                                                                                                gccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTaa     >  2:1839303/1‑150 (MQ=255)
                                                                                                                                                gccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTaa     >  2:1103987/1‑150 (MQ=255)
                                                                                                                                                  cgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACa   <  1:368488/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGGGTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAG  >  NC_000913/641208‑641504

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