Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,673,783 G→A 7.1% intergenic (+15/+33) yphA → / ← yphB putative inner membrane protein/putative aldose 1‑epimerase YphB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,673,7830GA 7.1% 352.4 / 10.9 114intergenic (+15/+33)yphA/yphBputative inner membrane protein/putative aldose 1‑epimerase YphB
Reads supporting (aligned to +/- strand):  ref base G (51/53);  new base A (5/3);  total (57/57)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 119 total aligned reads displayed.

CACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTC  >  NC_000913/2673634‑2673932
                                                                                                                                                     |                                                                                                                                                     
cACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATg                                                                                                                                                       >  2:468260/1‑150 (MQ=255)
cACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATg                                                                                                                                                       >  2:657884/1‑150 (MQ=255)
cACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCACGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATg                                                                                                                                                       >  2:797517/1‑150 (MQ=255)
    aTTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcggc                                                                                                                                                   <  2:609000/150‑1 (MQ=255)
       aCTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTg                                                                                                                                                <  1:1216352/150‑1 (MQ=255)
        cTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGa                                                                                                                                               >  2:42502/1‑150 (MQ=255)
        cTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGa                                                                                                                                               <  2:698883/150‑1 (MQ=255)
        cTGGGATATGACCGGCGATGCGGTTGAGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGa                                                                                                                                               >  1:324826/1‑150 (MQ=255)
             atatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                                                                          >  1:732093/1‑150 (MQ=255)
             atatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCACGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                                                                          <  1:797517/150‑1 (MQ=255)
              tatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCTCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATTCGGCGTGAACGCCt                                                                                                                                         <  2:954042/150‑1 (MQ=255)
              tatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCt                                                                                                                                         >  1:568075/1‑150 (MQ=255)
               atGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                        >  2:1238728/1‑150 (MQ=255)
                  aCCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGGAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                     >  2:1090011/1‑150 (MQ=255)
                  aCCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                     >  1:400421/1‑150 (MQ=255)
                     ggCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                  >  1:1156100/1‑150 (MQ=255)
                     ggCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                  >  2:849782/1‑150 (MQ=255)
                        gATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                               <  2:162552/150‑1 (MQ=255)
                         aTGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCAtt                                                                                                                              <  1:633742/150‑1 (MQ=255)
                          tGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCAttt                                                                                                                             <  1:685626/150‑1 (MQ=255)
                          tGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCAttt                                                                                                                             <  2:131535/150‑1 (MQ=255)
                          tGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCAttt                                                                                                                             <  1:1366889/150‑1 (MQ=255)
                          tGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCAttt                                                                                                                             <  2:876661/150‑1 (MQ=255)
                          tGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAtct                                                                                                                             <  1:1369024/150‑3 (MQ=39)
                              gTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaata                                                                                                                         <  1:1115080/150‑1 (MQ=255)
                               ttGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATTGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaataa                                                                                                                        <  1:214205/150‑1 (MQ=255)
                                tGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATaaa                                                                                                                       <  1:397105/150‑1 (MQ=255)
                                   gCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTa                                                                                                                    <  1:1485704/150‑1 (MQ=255)
                                          atGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                                                                             <  1:38522/150‑1 (MQ=255)
                                           tGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAct                                                                                                            >  1:1266447/1‑150 (MQ=255)
                                            gATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCActc                                                                                                           >  2:511123/1‑150 (MQ=255)
                                            gATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCActc                                                                                                           >  1:48639/1‑150 (MQ=255)
                                             aTTAATTTCTGAAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCActct                                                                                                          >  2:1535865/1‑150 (MQ=255)
                                                aaTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                                                                       <  1:657884/150‑1 (MQ=255)
                                                aaTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                                                                       <  2:115378/150‑1 (MQ=255)
                                                aaTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                                                                       >  2:1265216/1‑150 (MQ=255)
                                                 aTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAccc                                                                                                      >  2:900862/1‑150 (MQ=255)
                                                   ttCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGc                                                                                                    >  1:321645/1‑150 (MQ=255)
                                                    tCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCa                                                                                                   <  2:456938/150‑1 (MQ=255)
                                                    tCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCa                                                                                                   <  1:487842/150‑1 (MQ=255)
                                                      tGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAAc                                                                                                 >  2:232491/1‑150 (MQ=255)
                                                       ggAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACg                                                                                                <  2:1156100/150‑1 (MQ=255)
                                                       ggAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGAAATTACAGGGCCGGGGGAAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACg                                                                                                >  2:325859/1‑150 (MQ=255)
                                                            aaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtc                                                                                           <  1:924045/150‑1 (MQ=255)
                                                                tCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTg                                                                                       <  2:568075/150‑1 (MQ=255)
                                                                tCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTg                                                                                       <  1:23634/150‑1 (MQ=255)
                                                                tCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTg                                                                                       >  1:333547/1‑150 (MQ=255)
                                                                 cAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGa                                                                                      >  1:1247394/1‑150 (MQ=255)
                                                                 cAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGa                                                                                      <  2:1111107/150‑1 (MQ=255)
                                                                 cAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGa                                                                                      <  2:48639/150‑1 (MQ=255)
                                                                  aGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGAGCGGCGTTAGGAAAGAGGCCGGATGCGGCGGGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGaa                                                                                     <  2:774371/150‑1 (MQ=255)
                                                                   gTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAg                                                                                    >  1:825529/1‑150 (MQ=255)
                                                                      tCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttg                                                                                 >  2:635030/1‑150 (MQ=255)
                                                                       cGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCGCTCGCTCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCCACGGCATCTCTGAAgttgt                                                                                >  1:247670/1‑150 (MQ=255)
                                                                        gCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttgtt                                                                               <  1:656069/150‑1 (MQ=255)
                                                                          tGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGa                                                                             <  1:1192579/150‑1 (MQ=255)
                                                                           ggCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGAt                                                                            >  2:1348437/1‑150 (MQ=255)
                                                                            gcgcGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGAtt                                                                           >  1:528654/1‑150 (MQ=255)
                                                                                  tCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTg                                                                     <  1:849782/150‑1 (MQ=255)
                                                                                  tCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTCCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTg                                                                     <  2:1397039/150‑1 (MQ=255)
                                                                                       cTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCt                                                                <  1:1478517/150‑1 (MQ=255)
                                                                                          ttGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc                                                             <  1:1348437/150‑1 (MQ=255)
                                                                                            ggCAAATACCGGGCCGCGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGACGGATCCGGCGTGAACGCCTTATACGCCATTTAATAAATTACAGCCACTCTACCCGCAACTACATCTCAGAACTTGTTGATTCCCATGGTGATAgcgc                                                           >  2:1510296/1‑150 (MQ=255)
                                                                                                  tACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATgag                                                     >  1:82783/1‑150 (MQ=255)
                                                                                                  tACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATgag                                                     >  2:110124/1‑150 (MQ=255)
                                                                                                   aCCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc                                                             <  1:361456/141‑1 (MQ=255)
                                                                                                   aCCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATgaga                                                    <  1:468260/150‑1 (MQ=255)
                                                                                                    ccgggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGAt                                                   <  1:758428/150‑1 (MQ=255)
                                                                                                    ccgggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGAt                                                   <  2:137491/150‑1 (MQ=255)
                                                                                                    ccgggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGAt                                                   <  1:569160/150‑1 (MQ=255)
                                                                                                      gggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCa                                                 >  1:473808/1‑150 (MQ=255)
                                                                                                        gccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCAcc                                               >  2:1227589/1‑150 (MQ=255)
                                                                                                             gggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCtt                                          >  1:53479/1‑150 (MQ=255)
                                                                                                              ggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTc                                         >  1:44909/1‑150 (MQ=255)
                                                                                                               gCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTcc                                        >  2:264502/1‑150 (MQ=255)
                                                                                                                        ctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatga                               >  2:655219/1‑150 (MQ=255)
                                                                                                                        ctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatga                               >  2:861187/1‑150 (MQ=255)
                                                                                                                           gATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatgatga                            >  2:615549/1‑150 (MQ=255)
                                                                                                                            aTCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatgatgat                           <  2:338535/150‑1 (MQ=255)
                                                                                                                            aTCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatgatgat                           <  1:753018/150‑1 (MQ=255)
                                                                                                                              cGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGatca                         <  2:324826/150‑1 (MQ=255)
                                                                                                                              cGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGatca                         >  2:1257094/1‑150 (MQ=255)
                                                                                                                               ggCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGatcat                        >  1:1236460/1‑150 (MQ=255)
                                                                                                                                 cGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACGCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATcc                      <  1:1467799/150‑1 (MQ=255)
                                                                                                                                    tAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgc                   >  2:678957/1‑150 (MQ=255)
                                                                                                                                     aGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgcg                  >  2:927546/1‑150 (MQ=255)
                                                                                                                                      ggAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgcgc                 >  1:1295152/1‑150 (MQ=255)
                                                                                                                                      ggAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgcgc                 >  2:1492671/1‑150 (MQ=255)
                                                                                                                                        aaaGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCAt               <  1:819507/150‑1 (MQ=255)
                                                                                                                                        aaaGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCAt               <  1:814097/150‑1 (MQ=255)
                                                                                                                                         aaGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATg              <  1:1090011/150‑1 (MQ=255)
                                                                                                                                          aGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGATGATCATCATCGCCCggaaggggg             <  1:1021280/150‑8 (MQ=255)
                                                                                                                                           gATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCGCCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGc            <  1:1238728/150‑1 (MQ=255)
                                                                                                                                            aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCt           <  1:110124/150‑1 (MQ=255)
                                                                                                                                             tGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTc          >  1:400129/1‑150 (MQ=255)
                                                                                                                                             tGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTc          <  1:1227589/150‑1 (MQ=255)
                                                                                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCa         >  2:117170/1‑150 (MQ=255)
                                                                                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCa         >  1:1512344/1‑150 (MQ=255)
                                                                                                                                               ccGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCgg                                      <  2:1536648/120‑1 (MQ=38)
                                                                                                                                               ccGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCgg                                      >  1:1536648/1‑120 (MQ=38)
                                                                                                                                                     |                                                                                                                                                     
CACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTTC  >  NC_000913/2673634‑2673932

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