Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,189,606 T→G 6.8% intergenic (+33/‑180) rpoC → / → yjaZ RNA polymerase subunit beta'/conserved protein YjaZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,189,6060TG 6.8% 647.4 / 32.3 218intergenic (+33/‑180)rpoC/yjaZRNA polymerase subunit beta'/conserved protein YjaZ
Reads supporting (aligned to +/- strand):  ref base T (105/97);  new base G (5/10);  total (111/107)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 220 total aligned reads displayed.

ATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCGATAC  >  NC_000913/4189457‑4189754
                                                                                                                                                     |                                                                                                                                                    
    gTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTcggc                                                                                                                                                  >  1:2250586/1‑150 (MQ=255)
           tgctgcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGtttt                                                                                                                                           <  1:452719/150‑1 (MQ=255)
           tgctgcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCCGAAGCGGGtttt                                                                                                                                           <  2:1509425/150‑1 (MQ=255)
             ctgcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGTCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCCGAAGCGGGttttta                                                                                                                                         <  2:1458493/150‑2 (MQ=255)
             ctgcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGTGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGtttttt                                                                                                                                         >  1:1743315/1‑150 (MQ=255)
             ctgcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGtttttt                                                                                                                                         >  1:538346/1‑150 (MQ=255)
              tgcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCCGAAGCGGGtttttac                                                                                                                                        <  1:2080767/150‑3 (MQ=255)
               gcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTTCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTa                                                                                                                                       <  2:801619/150‑1 (MQ=255)
               gcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCCGAAGCGGGtttttacg                                                                                                                                       <  1:2203955/150‑4 (MQ=255)
                cGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTAt                                                                                                                                      <  2:1769934/150‑1 (MQ=255)
                  ggTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgg                                                                                                                                    <  1:2114501/150‑1 (MQ=255)
                   gTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATggg                                                                                                                                   <  1:830354/150‑1 (MQ=255)
                      aaGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGGATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGTGTTTTTTTATgggggg                                                                                                                                <  1:26737/150‑1 (MQ=255)
                       aGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGa                                                                                                                               >  2:794780/1‑150 (MQ=255)
                       aGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGa                                                                                                                               >  1:276037/1‑150 (MQ=255)
                        gCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAg                                                                                                                              <  2:484413/150‑1 (MQ=255)
                            cGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTa                                                                                                                          >  1:1808051/1‑150 (MQ=255)
                              gCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTGGGCGGGTTTTTTTATGGGGGGAGTGTAgg                                                                                                                        <  2:1405463/150‑1 (MQ=255)
                               cTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAggg                                                                                                                       >  1:779553/1‑150 (MQ=255)
                                   aCCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAag                                                                                                                   <  1:775531/150‑1 (MQ=255)
                                    ccGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAaga                                                                                                                  >  2:2111891/1‑150 (MQ=255)
                                      gCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGc                                                                                                                <  1:364057/150‑1 (MQ=255)
                                        aGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGAAt                                                                                                              >  1:975136/1‑150 (MQ=255)
                                             aCTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTc                                                                                                         <  2:490120/150‑1 (MQ=255)
                                              cTGCAGAAGACGCATCTGCCAGCCTGGCAGACCTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGTAGTTTAGGGAAAGAGCAttgggca                                                                                                        >  2:2474492/1‑145 (MQ=255)
                                               tGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAg                                                                                                       <  1:791834/150‑1 (MQ=255)
                                               tGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAg                                                                                                       <  2:2263371/150‑1 (MQ=255)
                                                gCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGa                                                                                                      >  1:1851776/1‑150 (MQ=255)
                                                gCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGa                                                                                                      >  1:1589608/1‑150 (MQ=255)
                                                        cagttCTGCCAGGCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAata                                                                                                   >  1:1923422/5‑145 (MQ=255)
                                                        cagttCTGCCAGGCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTaa                                                                                              >  2:631563/5‑150 (MQ=255)
                                                        cGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTaa                                                                                              >  2:1247119/1‑150 (MQ=255)
                                                         gCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAg                                                                                             <  2:2393542/150‑1 (MQ=255)
                                                         gCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAg                                                                                             <  1:1646986/150‑1 (MQ=255)
                                                           aTCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGa                                                                                           <  1:1329905/150‑1 (MQ=255)
                                                              tGCCAGCCTGGCAGAACTGCTGGACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAAtt                                                                                        >  2:788328/1‑150 (MQ=255)
                                                              tGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAAtt                                                                                        >  1:1687854/1‑150 (MQ=255)
                                                               gCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAAttt                                                                                       >  2:1998857/1‑150 (MQ=255)
                                                                ccAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTc                                                                                      <  1:50648/150‑1 (MQ=255)
                                                                ccAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTc                                                                                      >  1:1104667/1‑150 (MQ=255)
                                                                ccAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTc                                                                                      <  1:2488503/150‑1 (MQ=255)
                                                                   gCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGa                                                                                   <  2:463568/150‑1 (MQ=255)
                                                                           gAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCAt                                                                           >  2:2489371/1‑150 (MQ=255)
                                                                              ctgctgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAAt                                                                        <  1:1353699/150‑1 (MQ=255)
                                                                               tgctgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATc                                                                       <  2:1974768/150‑1 (MQ=255)
                                                                               tgctgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATc                                                                       >  1:1727903/1‑150 (MQ=255)
                                                                                gctgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCa                                                                      >  2:604976/1‑150 (MQ=255)
                                                                                  tgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAAt                                                                    >  1:1162535/1‑150 (MQ=255)
                                                                                   gAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATg                                                                   <  1:1371210/150‑1 (MQ=255)
                                                                                   gAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATg                                                                   >  2:608548/1‑150 (MQ=255)
                                                                                     aCGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTa                                                                 >  1:161356/1‑150 (MQ=255)
                                                                                      cGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTag                                                                <  2:1122155/150‑1 (MQ=255)
                                                                                      cGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTag                                                                <  1:910649/150‑1 (MQ=255)
                                                                                       gCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTaga                                                               <  2:227205/150‑1 (MQ=255)
                                                                                        cAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTagag                                                              >  1:1409323/1‑150 (MQ=255)
                                                                                          ggTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTa                                                                 >  1:2067430/1‑145 (MQ=255)
                                                                                              tGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACt                                                        <  2:595297/150‑1 (MQ=255)
                                                                                                gtcGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTaa                                                      <  1:1469194/148‑1 (MQ=255)
                                                                                                  gggTTCTGATAACGAGTAATCGTAAATCCGCAAATATCGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGTAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAata                                                    <  2:1289413/149‑1 (MQ=255)
                                                                                                  cGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAata                                                    <  1:1890200/150‑1 (MQ=255)
                                                                                                  cGGTACTGATAACGAGAAATCGTTAAACCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTGATGGGGGGAGTTTAGGGAAAGAGCATATGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGAAGAGATGACTAata                                                    >  2:1356281/1‑150 (MQ=255)
                                                                                                       cTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTg                                               <  2:308976/150‑1 (MQ=255)
                                                                                                         gATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGaa                                             <  2:17280/150‑1 (MQ=255)
                                                                                                         gATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGaa                                             >  2:2475458/1‑150 (MQ=255)
                                                                                                         gATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTCTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGaa                                             >  1:1212164/1‑150 (MQ=255)
                                                                                                           tAACCAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAac                                           <  2:538346/150‑1 (MQ=255)
                                                                                                             accaGTAATCGTTAATCCGCAAATAACGTAATAACCCGCTGCGGAGGGTTTTTTTATGGGGGGAGTTGAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGTTGACTAATATCCTGAAactg                                         >  2:1429186/4‑150 (MQ=255)
                                                                                                               gAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAACCAATGTAGAGATGACTAATATCCTGAAactgac                                       <  1:1098333/150‑1 (MQ=255)
                                                                                                                  tgaTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAActgactga                                    <  2:1033791/148‑1 (MQ=255)
                                                                                                                   aaTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGaa                                   >  1:632141/1‑150 (MQ=255)
                                                                                                                    aTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGTGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAAc                                  <  1:2399476/150‑1 (MQ=255)
                                                                                                                     tCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGATTGAACt                                 <  1:592195/150‑1 (MQ=255)
                                                                                                                     tCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACt                                 <  1:368503/150‑1 (MQ=255)
                                                                                                                      cGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTa                                >  2:2227845/1‑150 (MQ=255)
                                                                                                                        ttAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAAt                              <  2:2319263/150‑1 (MQ=255)
                                                                                                                        ttAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAAt                              >  1:1075629/1‑150 (MQ=255)
                                                                                                                          aaTCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTg                            >  2:1635059/1‑150 (MQ=255)
                                                                                                                          aaTCCGCAAATAACCTAAAAACCCGCTACGGCGGGTTTTTTTATGGGGAGAGTTTAGGCAAAGAGCATATCACAAAAAAATTAAGGAAATTCTCAATACTCATAATCAATGTAGAGAAGACTAATATCCTGAAACTCACTGAACTAAtag                            >  2:1847740/1‑148 (MQ=255)
                                                                                                                           aTCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGa                           >  1:1283776/1‑150 (MQ=255)
                                                                                                                           aTCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGa                           <  1:1537606/150‑1 (MQ=255)
                                                                                                                            tCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAg                          <  2:2482673/150‑1 (MQ=255)
                                                                                                                             ccGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTc                                                                                      <  2:779063/89‑1 (MQ=255)
                                                                                                                             ccGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGt                         <  1:2257908/150‑1 (MQ=255)
                                                                                                                             ccGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGt                         >  1:1444002/1‑150 (MQ=255)
                                                                                                                               gCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCa                       <  2:1687854/150‑1 (MQ=255)
                                                                                                                                cAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAATCTGACTGAACTAATTGAGTCaa                      >  2:1606930/1‑150 (MQ=255)
                                                                                                                                cAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCaa                      <  1:2465471/150‑1 (MQ=255)
                                                                                                                                      aCGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCgg                >  1:920539/1‑150 (MQ=255)
                                                                                                                                      aCGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCgg                >  1:787218/1‑150 (MQ=255)
                                                                                                                                      aCGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCgg                <  1:62024/150‑1 (MQ=255)
                                                                                                                                      aCGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCgg                >  2:349679/1‑150 (MQ=255)
                                                                                                                                      aCGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCgg                <  2:1726787/150‑1 (MQ=255)
                                                                                                                                      aCGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATCTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCgg                >  2:119988/1‑150 (MQ=255)
                                                                                                                                        aaaaaaaCCCGCCGAAGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCa              >  1:473542/3‑150 (MQ=255)
                                                                                                                                           aaaaCCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAgg           >  2:1930622/1‑150 (MQ=255)
                                                                                                                                            aaaCCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGa          >  1:2159197/1‑150 (MQ=255)
                                                                                                                                            aaaCCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGa          <  1:336108/150‑1 (MQ=255)
                                                                                                                                             aaCCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGAt         <  2:584733/150‑1 (MQ=255)
                                                                                                                                               cccGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTc       >  2:1539087/1‑150 (MQ=255)
                                                                                                                                                  gCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCGAt    <  2:1441923/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCGATAC  >  NC_000913/4189457‑4189754

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