Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 715,222 G→A 7.2% intergenic (+24/‑190) pgm → / → ybfP phosphoglucomutase/lipoprotein YbfP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913715,2220GA 7.2% 301.2 / 8.8 98intergenic (+24/‑190)pgm/ybfPphosphoglucomutase/lipoprotein YbfP
Reads supporting (aligned to +/- strand):  ref base G (47/43);  new base A (4/3);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

CGTCCGTCAGGCACGGAAGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAA  >  NC_000913/715076‑715366
                                                                                                                                                  |                                                                                                                                                
cgtccgtcAGGCACGGAAGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCg                                                                                                                                               >  1:458902/1‑150 (MQ=255)
 gtccgtcAGGCACGGAAGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGc                                                                                                                                              >  2:90841/1‑150 (MQ=255)
    cgtcAGGCACGGAAGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAg                                                                                                                                           >  1:557509/1‑150 (MQ=255)
       cAGGCACGGAAGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATc                                                                                                                                        >  2:458657/1‑150 (MQ=255)
         ggCACGGAAGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGc                                                                                                                                      >  1:1541289/1‑150 (MQ=255)
                  gACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGATCTTGCGACCGCCCttttttt                                                                                                                             <  2:389278/150‑2 (MQ=255)
                  gACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGATCTTGCGACCGCCCttttttt                                                                                                                             <  2:1103196/150‑2 (MQ=255)
                    cGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACGTTTAATAAAAAAAGGGCGGTCGCAAGATCg                                                                                                                                       >  2:427781/1‑138 (MQ=255)
                    cGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACGTTTAATAAAAAAAGGGCGGTCGCAAGATCg                                                                                                                                       <  1:427781/138‑1 (MQ=255)
                    cGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACg                                                                                                                           <  2:708693/150‑1 (MQ=255)
                       atatAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGATCTTGCGACCGCCCtttttttattaa                                                                                                                        <  2:1391587/150‑7 (MQ=38)
                        tataAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATg                                                                                                                       <  2:1050691/150‑1 (MQ=255)
                        tataAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATg                                                                                                                       >  1:1014408/1‑150 (MQ=255)
                         ataAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGa                                                                                                                      >  2:7845/1‑150 (MQ=255)
                            aGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACaa                                                                                                                   <  1:1163245/150‑1 (MQ=255)
                             gATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACaaa                                                                                                                  >  1:500025/1‑150 (MQ=255)
                              aTCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAac                                                                                                                 <  2:1005419/150‑1 (MQ=255)
                               tCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAaca                                                                                                                >  1:819582/1‑150 (MQ=255)
                               tCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAaca                                                                                                                >  1:117037/1‑150 (MQ=255)
                                cTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAacac                                                                                                               >  2:67855/1‑150 (MQ=255)
                                 tACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAAATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAAcaca                                                                                                              >  1:827258/1‑150 (MQ=255)
                                     gCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAAt                                                                                                          >  2:1590914/1‑150 (MQ=255)
                                     gCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAAt                                                                                                          >  2:143057/1‑150 (MQ=255)
                                      cGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAAtt                                                                                                         <  1:465087/150‑1 (MQ=255)
                                       gAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTg                                                                                                        >  2:554270/1‑150 (MQ=255)
                                         aaGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGcc                                                                                                      <  1:660714/150‑1 (MQ=255)
                                         aaGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGcc                                                                                                      <  1:1312017/150‑1 (MQ=255)
                                          aGCTTCCTCGTTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCt                                                                                                     <  2:1133127/150‑1 (MQ=255)
                                           gCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTg                                                                                                    >  2:727541/1‑150 (MQ=255)
                                             ttCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGAt                                                                                                  <  2:699989/150‑1 (MQ=255)
                                                cTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATgcg                                                                                               <  1:1117730/150‑1 (MQ=255)
                                                cTCGGTGAAGAACATCGAAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTTAAAAACGCTTAAACACATTTAATAAAAAAAGGGCGGTCTCAAGATCGACCTTTTTTACGTATGACAAACACATAATTGCCTGATgcg                                                                                               <  2:237955/150‑1 (MQ=255)
                                                 tCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATgcgc                                                                                              >  2:1081918/1‑150 (MQ=255)
                                                      gaagaaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGc                                                                                         <  2:458902/150‑1 (MQ=255)
                                                       aagaaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCt                                                                                        <  2:6982/150‑1 (MQ=255)
                                                        agaaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCtt                                                                                       >  1:1449883/1‑150 (MQ=255)
                                                         gaaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTa                                                                                      >  2:1628207/1‑150 (MQ=255)
                                                           aCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATc                                                                                    <  2:1070336/150‑1 (MQ=255)
                                                              tCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAgg                                                                                 >  1:1077279/1‑150 (MQ=255)
                                                               cGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGc                                                                                <  2:526260/150‑1 (MQ=255)
                                                                gCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCCGGcc                                                                               <  2:981228/150‑1 (MQ=255)
                                                                gCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGcc                                                                               <  1:563448/150‑1 (MQ=255)
                                                                   aGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                            >  2:23243/1‑150 (MQ=255)
                                                                          tGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATg                                                                     >  2:669638/1‑150 (MQ=255)
                                                                             gaaagaaGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTg                                                                  >  1:1572654/1‑150 (MQ=255)
                                                                              aaagaaGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGc                                                                 <  1:424548/150‑1 (MQ=255)
                                                                                agaagCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCaa                                                               <  1:1590914/150‑1 (MQ=255)
                                                                                  aagCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAata                                                             >  1:1209476/1‑150 (MQ=255)
                                                                                  aagCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAata                                                             <  1:67855/150‑1 (MQ=255)
                                                                                      ggTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATAtt                                                         >  1:278512/1‑150 (MQ=255)
                                                                                       gTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTg                                                        <  1:859813/150‑1 (MQ=255)
                                                                                        ttGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGa                                                       <  1:706334/150‑1 (MQ=255)
                                                                                        ttGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGa                                                       <  1:1218323/150‑1 (MQ=255)
                                                                                         tGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGaa                                                      <  1:7845/150‑1 (MQ=255)
                                                                                         tGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGaa                                                      >  1:936390/1‑150 (MQ=255)
                                                                                         tGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGaa                                                      <  2:1014408/150‑1 (MQ=255)
                                                                                          gagaTTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGTCGGTCGCAAGATCGCCCTTTTTTACGTATGACATACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAAt                                                     <  1:722983/150‑1 (MQ=255)
                                                                                           agaTTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAAtt                                                    <  1:1337069/150‑1 (MQ=255)
                                                                                               tGTTAGCGAAGTTCTTAAAAACTCGTAAAGACTTTTTATAAAAAAAGGGCGGTCTCAAGATCGGCGTTTTTTACGTATGACAAAGACTTTATTGCCTGTTGGGCTACGCTTATCAGGCCTTCGAGGATGGTGCAATATATTGAATTTatg                                                <  2:834284/150‑3 (MQ=255)
                                                                                               tGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAg                                                <  2:113891/150‑1 (MQ=255)
                                                                                                  tAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGa                                             >  1:44071/1‑150 (MQ=255)
                                                                                                  tAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGa                                             >  1:38061/1‑150 (MQ=255)
                                                                                                   aGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGAt                                            >  1:344409/1‑150 (MQ=255)
                                                                                                   aGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGAt                                            >  2:74129/1‑150 (MQ=255)
                                                                                                    gCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGAtt                                           >  2:260611/1‑150 (MQ=255)
                                                                                                      gAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGAtttt                                         >  1:610200/1‑150 (MQ=255)
                                                                                                            cTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGc                                   >  1:1578880/1‑150 (MQ=255)
                                                                                                             tGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGcc                                  >  2:743954/1‑150 (MQ=255)
                                                                                                                  aaCGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATa                             <  2:101178/150‑1 (MQ=255)
                                                                                                                  aaCGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATa                             <  2:315366/150‑1 (MQ=255)
                                                                                                                    cgcgTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAg                           <  2:344409/150‑1 (MQ=255)
                                                                                                                      cgTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGc                         >  2:450430/1‑150 (MQ=255)
                                                                                                                        tAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGt                       >  1:1550599/1‑150 (MQ=255)
                                                                                                                         aaaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCGGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGtt                      >  1:1520769/1‑150 (MQ=255)
                                                                                                                         aaaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGtt                      >  2:930021/1‑150 (MQ=255)
                                                                                                                         aaaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGtt                      >  2:1227799/1‑150 (MQ=255)
                                                                                                                         aaaCACATTTAATAAAAAAAGGGCGCTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCCGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGtt                      >  2:7354/1‑150 (MQ=255)
                                                                                                                          aaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTc                     <  1:260611/150‑1 (MQ=255)
                                                                                                                          aaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTc                     <  2:117037/150‑1 (MQ=255)
                                                                                                                              caTTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATTTATTGAATTTTAGCGATTTTGTTGGCCGGGTAAGGCGTTAAcgg                 >  2:711132/1‑149 (MQ=37)
                                                                                                                               aTTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCAcgcc                >  1:158288/1‑150 (MQ=255)
                                                                                                                               aTTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCAcgcc                <  2:492366/150‑1 (MQ=255)
                                                                                                                                tttAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCAcgccg               >  1:1249613/1‑150 (MQ=255)
                                                                                                                                 ttAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCAcgccgc              >  2:668300/1‑150 (MQ=255)
                                                                                                                                 atacgtaaaaaaGGGCGATCTTGCGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCAcgccgc              >  2:134865/7‑150 (MQ=21)
                                                                                                                                  taataaAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCa             <  1:1628207/150‑1 (MQ=255)
                                                                                                                                    cgtaaaaaaGGGCGATCTTGCGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc           >  2:1504036/4‑150 (MQ=21)
                                                                                                                                     taaAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc          >  1:664440/1‑150 (MQ=255)
                                                                                                                                      taaaaaaGGGCGATCTTGCGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg         >  2:1630942/2‑150 (MQ=21)
                                                                                                                                      aaaaaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTT‑TAGGCCGGATAAGGCGTTCACGCCGCATCCgg        >  2:149786/1‑150 (MQ=255)
                                                                                                                                       aaaaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg        <  1:458657/150‑1 (MQ=25)
                                                                                                                                       aaaaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg        <  2:1077279/150‑1 (MQ=255)
                                                                                                                                        aaaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc       <  1:1301905/150‑1 (MQ=25)
                                                                                                                                        aaaaaGGGCGATCTTGCGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGATTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc       >  2:1228345/1‑150 (MQ=21)
                                                                                                                                         aaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa      <  1:824677/150‑1 (MQ=255)
                                                                                                                                          aaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa     <  1:1081918/150‑1 (MQ=255)
                                                                                                                                            aGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctac   <  2:610200/150‑5 (MQ=25)
                                                                                                                                             gggCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaaaaa  >  1:610973/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CGTCCGTCAGGCACGGAAGACGCATATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAA  >  NC_000913/715076‑715366

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: ‑