Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,908,506 C→G 5.3% intergenic (‑140/+43) pstB ← / ← pstA phosphate ABC transporter ATPase/phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,908,5060CG 5.3% 2259.3 / 27.3 671intergenic (‑140/+43)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (375/248);  new base G (4/42);  total (379/292)
Fisher's exact test for biased strand distribution p-value = 3.40e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 705 total aligned reads displayed.

CCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCG  >  NC_000913/3908357‑3908655
                                                                                                                                                       |                                                                                                                                                     
ccATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCctc                                                                                                                                                       <  2:1724426/150‑1 (MQ=255)
 cATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCctct                                                                                                                                                      <  1:7817097/150‑1 (MQ=255)
  aTACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTg                                                                                                                                                     >  2:1236115/1‑150 (MQ=255)
   tACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgg                                                                                                                                                    <  1:6232008/150‑1 (MQ=255)
    aCTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   <  2:6366373/150‑1 (MQ=255)
      tCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                 >  1:7012188/1‑150 (MQ=255)
        aTTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               <  1:7649656/150‑1 (MQ=255)
         ttGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              <  2:162677/150‑1 (MQ=255)
          tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                             <  2:4373624/150‑1 (MQ=255)
          tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                             >  1:2823089/1‑150 (MQ=255)
          tGCAATCTTGTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTCAAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              <  2:1404660/150‑1 (MQ=255)
           gCAATCTGTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                            <  2:3670564/150‑1 (MQ=255)
              aTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                         <  2:2416747/150‑1 (MQ=255)
               tctggtCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                        <  1:2697431/145‑1 (MQ=255)
                 ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                      >  2:793793/1‑150 (MQ=255)
                  tttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                     >  1:1939868/1‑150 (MQ=255)
                    tCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                   >  2:7532899/1‑150 (MQ=255)
                    tCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑ATGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                   >  1:2550899/1‑150 (MQ=255)
                               ccTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCg                                                                                                                        >  2:7439015/1‑150 (MQ=255)
                               ccTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCg                                                                                                                        >  1:7535299/1‑150 (MQ=255)
                               ccTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCg                                                                                                                        >  2:2140488/1‑150 (MQ=255)
                               ccTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTGAGGGTGAGGGGACAACCg                                                                                                                        >  2:6035143/1‑150 (MQ=255)
                                  gATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtgt                                                                                                                     >  2:6531959/1‑150 (MQ=255)
                                   aTGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTc                                                                                                                    >  1:2044151/1‑150 (MQ=255)
                                     gCAGTTCGCTACTCTGGCTTACCGTTGGCCGGATGAGGCGTTCACGCCGCATCGGGCAAAGTTTTAACG‑TTGTACATAGTC‑AAAACGGCGCTCGGTCGGTCCCCTCGCCCCTGTGGGGTGTGGGTTAGGGTGCGGGGAAAACCGTGTCAg                                                                                                                  <  2:7649606/150‑1 (MQ=255)
                                     gCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACA‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCACTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAg                                                                                                                  >  1:5597660/1‑150 (MQ=255)
                                         ttCGCTAATCTGTCTATCCGTAGGCCTTATAAGGCTTACACGCCGAATCGTGCAAA‑TTTTAACCGGTTTACAGAGAA‑AAAATGGGGCTGGGTGGGCCCCCTCGCCCCTGGGGGTTGTTGTTTTGGGTGCGGGGAAAACCGTGTCAGCAat                                                                                                              <  2:7585502/150‑1 (MQ=255)
                                         ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                              >  1:4153107/1‑150 (MQ=255)
                                            gCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATc                                                                                                           >  1:5894349/1‑150 (MQ=255)
                                               aaTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                                                        <  2:5085357/150‑1 (MQ=255)
                                                    ggCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCCTCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgtgt                                                                                                   <  2:6126451/150‑1 (MQ=255)
                                                     gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                                  >  2:679690/1‑150 (MQ=255)
                                                      cTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGCGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGGGCGGGTGAGGGGACACCCGTGTCAGCACTATCAACCGTGttt                                                                                                 >  2:5662318/1‑150 (MQ=39)
                                                          ccGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttc                                                                                             >  1:1318355/1‑150 (MQ=255)
                                                            gTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt                                                                                           >  1:4581255/1‑150 (MQ=255)
                                                              cggCCGTATCTGGGGTTCACGCCGCACCCGGAAAAGTGTTACCC‑GTGTACAGAGAC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGTGGTTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                         <  1:1754617/149‑1 (MQ=255)
                                                                 ccGGATGCGGCGTGAACGCCTGATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaa                                                                                      >  2:8188501/1‑150 (MQ=38)
                                                                  cGGATCAGGCGTTCACGCCGAATCCGGAAAAGTGTTAACC‑GTGTACATAGTA‑AAAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGTAGTGGTTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaa                                                                                     <  1:6759618/150‑1 (MQ=255)
                                                                    gATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               >  1:1203924/1‑90 (MQ=255)
                                                                    gATCAGGCGTTCACGCCGCATCCGGAAATGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCACGCCCCTCTGGGGTGTGGTTTGGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaaa                                                                                   <  1:845047/150‑1 (MQ=255)
                                                                     aTCAGGCGTTCCCGCCGCACCCGGAAATGTTTTACCC‑GTGTACAGAGTA‑AAAACGGCGCTCGGTCGGTCCCCTCGCCCCTGTGGGTTGTGGTTTCGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaac                                                                                  <  1:3485221/150‑1 (MQ=255)
                                                                     aTCAGGCGTACACGCCGCATCCGCAAAAGTCTAAACC‑GTCTACAGACAA‑AAAAAGGCGCCCGGTCGGTCCCCACGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaac                                                                                  <  2:2604439/150‑1 (MQ=37)
                                                                       cAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑TTGTACTGAGTC‑AAATAGGCGCTCGGTCGTTCCCCTCGCCCCTCGGGGTAGTGGTTTAGGGTGAGGGGAAAACCGTGTCAGCTATATCAACCGTGTTTATTCTTCGCAAAaacaa                                                                                <  2:1048903/150‑1 (MQ=255)
                                                                       cAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaa                                                                                >  2:304007/1‑150 (MQ=255)
                                                                       cAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaa                                                                                >  1:7786762/1‑150 (MQ=255)
                                                                         ggCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                                              >  1:532608/1‑150 (MQ=255)
                                                                            gtacACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGAC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgc                                                                           <  1:6173591/147‑1 (MQ=255)
                                                                             ttCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                                                                          >  1:4882874/1‑150 (MQ=255)
                                                                             ttCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                                                                          >  2:3401733/1‑150 (MQ=255)
                                                                              tCACGCCGCATCCGGCAAAGTGTTAACC‑TTTAACAGAGTA‑AAAAAGGCGCTGGTTGGGTCCCCTCGCCCCTCTGGGTAGTGGGTTTGGTTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACgcgcgc                                                                         <  2:1536637/150‑1 (MQ=255)
                                                                                ccgccccACCCGGAAATGTGTTAACG‑GTGTACAGAGCA‑AAAACGGCGCTCGGTCGGCCCCCTCGCCCCTGTGGGGTGTGGTTTTGGGTGTGGGCAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCa                                                                       <  1:5494586/149‑1 (MQ=255)
                                                                                aCGCCGCATCCGGCAATTTTTTAACG‑GTTTCCTGAGTC‑AAAACGGGGCGCGGTGGGGCCCCTCGCCCCTGTGGGGTGTGTGTTTGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCa                                                                       <  2:3115541/150‑1 (MQ=255)
                                                                                aCGCCGCATCCGGCAATGTGTTAACC‑TTTTACTGCGTC‑AAAAGGGCGCTGGGTCGGTCCCCGCGCCCCTCTGGGGTGAGGGTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCa                                                                       <  2:7114972/150‑1 (MQ=255)
                                                                                  gccgcATCCGTAATTGTTTTAACC‑GTGTACTTATAC‑AAAACGGCGCTGGGTGGGTCCCCTCCCCCCTGTGGGGTGTGTGTTTGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                                                     <  1:6371880/150‑1 (MQ=255)
                                                                                       aTCCGGCATTGTTTTACCG‑GTGTACTGATTA‑AACATGGCGCTGGGTGGGACCCCTCGCCCCTGTGGGTTGTGGTTTTGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGtt                                                                <  2:2823845/150‑1 (MQ=255)
                                                                                       aTCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGtt                                                                >  2:3764682/1‑150 (MQ=255)
                                                                                       aTCCGGAAAAGTGTTAACG‑GTTTACTGAGTC‑AAAACGGCGCTGGGTCGGCCCCCTCGCCCCTGTGGGGTTTGGGTTTGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGtt                                                                <  2:2261759/150‑1 (MQ=255)
                                                                                           ggCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagc                                                            >  1:3341771/1‑150 (MQ=255)
                                                                                             cAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcag                                                          >  1:5707896/1‑150 (MQ=255)
                                                                                              aaaGTGATAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAAAGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagc                                                         >  1:7171790/1‑150 (MQ=255)
                                                                                                   gtTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTAc                                                    >  2:1173819/1‑150 (MQ=255)
                                                                                                        cc‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcaca                                               >  1:1580778/1‑150 (MQ=255)
                                                                                                        cc‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATAACAACCGTGTTTATTCTTCGCAACAACAACGCGCGCCAGAAGGGTCAGCCGCAGTACGcacc                                               >  2:5828993/1‑149 (MQ=38)
                                                                                                            tgtACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAggg                                            >  1:8175773/1‑150 (MQ=255)
                                                                                                            tgaaccgagTC‑AAAATGGCGCTCGGTCGGGCCCCGCGCCCCTCGGGGGTGTGGGTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAggg                                            <  1:2425734/144‑1 (MQ=255)
                                                                                                               aCAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGGTTATGCTTCGCACAAACCACGCGCGCCCGAATGGTCCGCAGCAGTCCGCACAGGGtcc                                         >  2:5314552/1‑148 (MQ=37)
                                                                                                                  gagTC‑AAAAAGGCGCTGGGTCGGGCCCCGCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGa                                      <  1:1029089/150‑1 (MQ=255)
                                                                                                                     tC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCa                                   >  2:4720696/1‑150 (MQ=255)
                                                                                                                       aaaaacGGCGCTCGTTCGGTCCCCTCGCCCCTCTGTGTAGTGGTTTTGGGTGAGGGCAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCaa                                  <  2:2900628/149‑1 (MQ=255)
                                                                                                                        aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                 >  2:3245215/1‑150 (MQ=255)
                                                                                                                        aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                 >  2:3994480/1‑150 (MQ=255)
                                                                                                                        aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                 >  2:7215981/1‑150 (MQ=255)
                                                                                                                        aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                 >  2:6889753/1‑150 (MQ=255)
                                                                                                                         caaaGGCGCTCGGTCGGCCCCCTCGCCCCTCGGGGGAGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATa                                <  1:4432096/149‑1 (MQ=255)
                                                                                                                         aaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATa                                >  1:5777164/1‑150 (MQ=255)
                                                                                                                             ggCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcccc                            <  1:790690/150‑1 (MQ=255)
                                                                                                                             ggCGCTCGGTCGGCCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcccc                            <  1:3010518/150‑1 (MQ=255)
                                                                                                                               cgcTCGGCCGGCCCCCCCCCCCCTCTGGGTAGAGGGTTCTGGTGAGGGCAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCgg                          <  1:2177675/150‑1 (MQ=255)
                                                                                                                                gcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGc                         >  2:4137665/1‑150 (MQ=255)
                                                                                                                                 cTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAACGTGTTTATTCTTAGCAAAAAAAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATAAATACCCCGGcc                        >  1:2339157/1‑150 (MQ=255)
                                                                                                                                    gggggggccccGCGCCCCTCTGGGGTGTGGGTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAg                     <  1:3728067/143‑1 (MQ=255)
                                                                                                                                     gtcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAgg                    >  2:7263086/1‑150 (MQ=255)
                                                                                                                                     accctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAgg                    >  2:5003460/8‑150 (MQ=255)
                                                                                                                                       cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGcc                  >  1:6920482/6‑150 (MQ=255)
                                                                                                                                          tcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAAt               >  1:5163826/1‑150 (MQ=255)
                                                                                                                                           ccccTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAAtt              <  2:1948200/150‑1 (MQ=255)
                                                                                                                                           ccccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAAtt              <  2:4274656/150‑1 (MQ=255)
                                                                                                                                           ccccGCGCCCCTCGGGGTAGTGGGTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAAtt              <  2:6237806/150‑1 (MQ=255)
                                                                                                                                            cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtg             >  1:2944509/1‑150 (MQ=255)
                                                                                                                                             ccTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTACGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTCCTGAACAATACCCCGGCCCAGGccccgtgc            >  2:7074767/1‑144 (MQ=39)
                                                                                                                                              cTCGCCCCTGTGGGGTGTGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct           <  1:6430340/150‑1 (MQ=255)
                                                                                                                                              cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct           >  2:5190720/1‑150 (MQ=255)
                                                                                                                                              cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct           >  1:7279969/1‑150 (MQ=255)
                                                                                                                                              cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct           >  1:4731337/1‑150 (MQ=255)
                                                                                                                                              cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct           >  2:806606/1‑150 (MQ=255)
                                                                                                                                              cTCGCCCCACTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct           >  1:3867540/1‑150 (MQ=255)
                                                                                                                                                cGCCCCTCGGGGGAGAGGTTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctgc         <  1:313396/150‑1 (MQ=255)
                                                                                                                                                     ctctGGGGTGCGGGTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAAGGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTc    <  1:139392/150‑1 (MQ=255)
                                                                                                                                                      tctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTcc   >  2:4521958/1‑150 (MQ=255)
                                                                                                                                                      tctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTcc   >  2:3513306/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACC‑GTGTACAGAGTC‑AAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCG  >  NC_000913/3908357‑3908655

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑