Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 715,227 A→G 6.4% intergenic (+29/‑185) pgm → / → ybfP phosphoglucomutase/lipoprotein YbfP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913715,2270AG 6.4% 284.7 / 7.4 93intergenic (+29/‑185)pgm/ybfPphosphoglucomutase/lipoprotein YbfP
Reads supporting (aligned to +/- strand):  ref base A (42/45);  new base G (4/2);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

TCCGTCAGGCACGGAAGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACA  >  NC_000913/715078‑715375
                                                                                                                                                     |                                                                                                                                                    
tccgtcAGGCACGGAAGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCa                                                                                                                                                      <  2:477515/150‑1 (MQ=255)
   gtcAGGCACGGAAGACACATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGa                                                                                                                                                   <  1:478012/150‑1 (MQ=255)
     cAGGCACGGAAGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATc                                                                                                                                                 >  2:1067993/1‑150 (MQ=255)
             gAAGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCtttt                                                                                                                                         >  2:699109/1‑150 (MQ=255)
              aaGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCttttt                                                                                                                                        <  1:987432/150‑1 (MQ=255)
               aGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCtttttt                                                                                                                                       <  1:645189/150‑1 (MQ=255)
                gACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGATCTTGCGACCGCCCttttttt                                                                                                                                      <  1:334415/150‑2 (MQ=255)
                   gCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGATCTTGCGACCGCCCtttttttatt                                                                                                                                   <  2:1230124/150‑5 (MQ=39)
                    cATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTa                                                                                                                                  <  1:27993/150‑1 (MQ=255)
                        taAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGAc                                                                                                                              >  1:788571/1‑150 (MQ=255)
                              cTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAacac                                                                                                                        >  2:100184/1‑150 (MQ=255)
                                 cTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACATACACAGa                                                                                                                     <  2:599413/150‑1 (MQ=255)
                                 cTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGa                                                                                                                     >  1:1011423/1‑150 (MQ=255)
                                  tGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGaa                                                                                                                    >  1:1401298/1‑150 (MQ=255)
                                     gAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTg                                                                                                                 <  1:530722/150‑1 (MQ=255)
                                     gAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTg                                                                                                                 <  1:394968/150‑1 (MQ=255)
                                              cTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATgcg                                                                                                        >  2:120677/1‑150 (MQ=255)
                                              cTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCATGATgcg                                                                                                        >  2:1216444/1‑150 (MQ=255)
                                               tCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATgcgc                                                                                                       <  2:788571/150‑1 (MQ=255)
                                                  gTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTAc                                                                                                    <  2:155423/150‑1 (MQ=255)
                                                    gaagaaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGc                                                                                                  <  2:941298/150‑1 (MQ=255)
                                                      agaaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACTTATGACAAACAAAGAATTGCATGATGCGATACGAtt                                                                                                >  2:274440/1‑150 (MQ=255)
                                                      agaaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCtt                                                                                                >  1:564325/1‑150 (MQ=255)
                                                      agaaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCtt                                                                                                >  1:708564/1‑150 (MQ=255)
                                                           aTCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAg                                                                                           >  2:1195431/1‑150 (MQ=255)
                                                           aTCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAg                                                                                           >  1:1120879/1‑150 (MQ=255)
                                                           aTCGCAAGCAGATGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAgg                                                                                          >  2:380268/1‑150 (MQ=255)
                                                            tCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAGGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAgg                                                                                          >  2:1120772/1‑150 (MQ=255)
                                                            tCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACGCATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAgg                                                                                          <  2:1218210/150‑1 (MQ=255)
                                                            tCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAgg                                                                                          <  2:1011423/150‑1 (MQ=255)
                                                             cGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGc                                                                                         >  2:436946/1‑150 (MQ=255)
                                                             cGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGc                                                                                         <  2:887033/150‑1 (MQ=255)
                                                              gCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTATGCTTATCAGGcc                                                                                        >  2:278259/1‑150 (MQ=255)
                                                              gCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGcc                                                                                        <  1:934784/150‑1 (MQ=255)
                                                                aaGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTa                                                                                      <  2:803915/150‑1 (MQ=255)
                                                                 aGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                     >  1:1140080/1‑150 (MQ=255)
                                                                 aGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                     <  2:174139/150‑1 (MQ=255)
                                                                  gCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACg                                                                                    <  1:1067993/150‑1 (MQ=255)
                                                                  gCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACg                                                                                    >  1:1285694/1‑150 (MQ=255)
                                                                   cAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGa                                                                                   <  2:1177299/150‑1 (MQ=255)
                                                                       ttGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGAt                                                                               >  1:481282/1‑150 (MQ=255)
                                                                         gagaAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTATTTTACGTATGACAAACACAGAATTCCCTGATGCGCTACGCTTATCAGGCCTACGAGGATgg                                                                             >  2:1217435/1‑150 (MQ=255)
                                                                           gaaagaaGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTg                                                                           >  2:312442/1‑150 (MQ=255)
                                                                            aaagaaGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGc                                                                          <  1:100184/150‑1 (MQ=255)
                                                                             aagaagCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCa                                                                         >  1:807885/1‑150 (MQ=255)
                                                                              agaagCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCaa                                                                        <  2:1274098/150‑1 (MQ=255)
                                                                               gaagCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAat                                                                       >  1:383837/1‑150 (MQ=255)
                                                                                 agCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAatat                                                                     >  2:721431/1‑150 (MQ=255)
                                                                                    ggTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAatatat                                                                   >  1:929380/1‑149 (MQ=255)
                                                                                    ggTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAatatat                                                                   <  2:929380/149‑1 (MQ=255)
                                                                                    ggTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATAtt                                                                  <  1:440147/150‑1 (MQ=255)
                                                                                      ttGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCTCTACGCTTATCAGGCCTACGAGGATGGTGCAATATAttta                                                                <  2:1012542/150‑3 (MQ=255)
                                                                                       tGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACg                                                                                    <  2:517607/129‑1 (MQ=255)
                                                                                       tGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACg                                                                                    >  1:517607/1‑129 (MQ=255)
                                                                                         agaTTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAAtt                                                             <  2:420650/150‑1 (MQ=255)
                                                                                           aTTGTTAGCGAAGTTCCGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTa                                                           <  1:790139/150‑1 (MQ=255)
                                                                                           aTTGTTAGCGAAGTTCCGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTa                                                           <  1:785149/150‑1 (MQ=255)
                                                                                            ttGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTaa                                                          >  1:1384889/1‑150 (MQ=255)
                                                                                                tAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGa                                                      >  2:721588/1‑150 (MQ=255)
                                                                                                 aGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGAt                                                     >  1:641064/1‑150 (MQ=255)
                                                                                                  gCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGAtt                                                    >  2:1306983/1‑150 (MQ=255)
                                                                                                  gCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGAtt                                                    >  2:1312196/1‑150 (MQ=255)
                                                                                                        ttCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGAt                                                                               >  2:707399/1‑117 (MQ=255)
                                                                                                        ttCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGAt                                                                               <  1:707399/117‑1 (MQ=255)
                                                                                                        ttCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAg                                              <  1:380268/150‑1 (MQ=255)
                                                                                                         tCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAgg                                             >  1:922424/1‑150 (MQ=255)
                                                                                                          cTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGc                                            >  2:161107/1‑150 (MQ=255)
                                                                                                           tGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGcc                                           >  2:1311447/1‑150 (MQ=255)
                                                                                                            gAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCg                                          <  2:1285694/150‑1 (MQ=255)
                                                                                                              aaaaCGCGTAAACTCATTTAATAAAAATAGGGCTTTCGCAAGATCTCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATTTTGAATTTAAGCGATTTTGTAGGCCGGa                                        <  1:1021429/150‑1 (MQ=37)
                                                                                                                   gcgTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAgg                                   <  1:120677/150‑1 (MQ=255)
                                                                                                                   gcgTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAgg                                   <  1:1216444/150‑1 (MQ=255)
                                                                                                                     gTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCg                                 <  2:241194/150‑1 (MQ=255)
                                                                                                                     gTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCg                                 <  1:274440/150‑1 (MQ=255)
                                                                                                                      tAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGt                                <  2:67084/150‑1 (MQ=255)
                                                                                                                      tAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGt                                <  2:807885/150‑1 (MQ=255)
                                                                                                                        aaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTc                              >  2:1149196/1‑150 (MQ=255)
                                                                                                                          cacaTTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCAc                            <  1:699109/150‑1 (MQ=255)
                                                                                                                            caTTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCAcgc                          >  1:1377676/1‑150 (MQ=255)
                                                                                                                             aTTTAATAAAAATAGGCAGGTCGCAAGATCGCACTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCTATTTTGTAGGCCGGATAAGGCGTTCAcgcc                         <  2:522492/150‑1 (MQ=25)
                                                                                                                             aTTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATa                                      >  2:809025/1‑137 (MQ=255)
                                                                                                                             aTTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATa                                      <  1:809025/137‑1 (MQ=255)
                                                                                                                               ttAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCAcgccgc                       <  2:1401298/150‑1 (MQ=255)
                                                                                                                                  cgtaaaaaaGGGCGATCTTGCGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                    >  2:862976/4‑150 (MQ=21)
                                                                                                                                  cgtaaaaaaGGGCGATCTTGCGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                    >  1:355244/4‑150 (MQ=21)
                                                                                                                                  ataaAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCGTATCAGGCCTACGGGGCTGGTGCAATATTTTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                    >  1:372561/1‑150 (MQ=255)
                                                                                                                                   taaAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATa                                      <  2:455842/131‑1 (MQ=255)
                                                                                                                                   taaAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATa                                      >  1:455842/1‑131 (MQ=255)
                                                                                                                                    taaaaaaGGGCGATCTTGCGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                  >  2:648545/2‑150 (MQ=21)
                                                                                                                                    taaaaaaGGGCGATCTTGCGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                  >  2:912990/2‑150 (MQ=21)
                                                                                                                                             gCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctacaaa         <  1:721588/150‑8 (MQ=25)
                                                                                                                                               ggTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAaacaac       <  2:383837/150‑1 (MQ=255)
                                                                                                                                                 tCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGa     <  2:641064/150‑1 (MQ=255)
                                                                                                                                                    cAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAaca  >  2:213405/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCCGTCAGGCACGGAAGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACA  >  NC_000913/715078‑715375

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: