Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,574,175 A→T 5.3% intergenic (‑72/+20) yicG ← / ← orn inner membrane protein/oligoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,574,1750AT 5.3% 481.5 / 10.5 143intergenic (‑72/+20)yicG/orninner membrane protein/oligoribonuclease
Reads supporting (aligned to +/- strand):  ref base A (68/66);  new base T (2/6);  total (70/73)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

CATGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGT  >  NC_002947/5574026‑5574324
                                                                                                                                                     |                                                                                                                                                     
cATGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACa                                                                                                                                                       <  2:5264907/150‑1 (MQ=255)
cATGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACa                                                                                                                                                       <  1:87385/150‑1 (MQ=255)
cATGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACa                                                                                                                                                       >  2:1435264/1‑150 (MQ=255)
 aTGCCGCGACGCCCGGCAGACAGCGTGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      <  2:1232105/150‑1 (MQ=255)
 aTGCCGCGACGCCCGGCAGACAGCGTGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      <  2:1780737/150‑1 (MQ=255)
 aTGCCGCGACGCCCGGCAGACAGCGTGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      <  2:1780890/150‑1 (MQ=255)
 aTGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGACAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      >  1:929939/1‑150 (MQ=255)
 aTGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      >  1:2784355/1‑150 (MQ=255)
 aTGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      >  1:3969229/1‑150 (MQ=255)
 aTGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      <  1:4030981/150‑1 (MQ=255)
 aTGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      >  1:3233305/1‑150 (MQ=255)
 aTGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      <  2:4500980/150‑1 (MQ=255)
 aTGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      <  1:133294/150‑1 (MQ=255)
   gCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAata                                                                                                                                                    >  2:1442610/1‑150 (MQ=255)
      gcgACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCt                                                                                                                                                 >  1:471126/1‑150 (MQ=255)
        gACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTgc                                                                                                                                               >  1:1364118/1‑150 (MQ=255)
            cccGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATATTGTATCTGCGccc                                                                                                                                           <  2:3121332/150‑2 (MQ=255)
            cccGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATATTGTATCTGCGccc                                                                                                                                           <  2:3121462/150‑2 (MQ=255)
                  aGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACCTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGa                                                                                                                                     <  2:5168859/150‑1 (MQ=255)
                          cgcCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGc                                                                                                                             <  1:3409524/150‑1 (MQ=255)
                           gcCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCt                                                                                                                            >  2:4566244/1‑150 (MQ=255)
                               gTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTtgat                                                                                                                        <  1:268805/150‑1 (MQ=255)
                                  aTGGCTTCTGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGaa                                                                                                                     <  2:4775654/150‑1 (MQ=255)
                                  aTGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGaa                                                                                                                     >  1:3818424/1‑150 (MQ=255)
                                  aTGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGaa                                                                                                                     <  1:3165651/150‑1 (MQ=255)
                                  aTGGCTTCGGCGGTAAGGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAACAGGGGGCGCAGATACCATATCTGCGCCGTTTCGACTCACAGCTTGATGaa                                                                                                                     <  2:236699/150‑1 (MQ=255)
                                   tGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAg                                                                                                                    <  2:266203/150‑1 (MQ=255)
                                   tGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAg                                                                                                                    >  2:1891555/1‑150 (MQ=255)
                                    ggCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGACTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGt                                                                                                                   <  1:504735/150‑1 (MQ=255)
                                              gTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATa                                                                                                         <  1:1551886/150‑1 (MQ=255)
                                               tAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAg                                                                                                        <  2:3389900/150‑1 (MQ=255)
                                                aaTGGCGAAAAGGTAGAGCATCAACAACATGGAGCGGGTCCGTGCAGGAAAGGCGCGCAGACTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACACACAGCTTGATGAAGTGCTCGCGATAGt                                                                                                       >  1:103134/1‑150 (MQ=255)
                                                            gTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCTCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCgg                                                                                           >  1:4782744/1‑150 (MQ=255)
                                                            gTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCgg                                                                                           >  1:2131999/1‑150 (MQ=255)
                                                            gTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCgg                                                                                           >  2:3961210/1‑150 (MQ=255)
                                                             tAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGc                                                                                          <  2:2131999/150‑1 (MQ=255)
                                                               gagCATCAACAACCTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcga                                                                                        <  2:2397221/150‑1 (MQ=255)
                                                               gagCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcga                                                                                        >  2:1268788/1‑150 (MQ=255)
                                                               gagCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcga                                                                                        <  2:5096595/150‑1 (MQ=255)
                                                                  cATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcgatcg                                                                                     >  1:3725533/1‑150 (MQ=255)
                                                                    tCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcgatcgat                                                                                   >  1:3060482/1‑150 (MQ=255)
                                                                      aacaacATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTc                                                                                 >  2:4812109/1‑150 (MQ=255)
                                                                        caacaTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTcgc                                                                               <  1:3901310/150‑1 (MQ=255)
                                                                          acaTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAACGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCgg                                                                             <  2:297106/150‑1 (MQ=255)
                                                                             tGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGata                                                                          <  2:3526302/150‑1 (MQ=255)
                                                                             tGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGata                                                                          <  1:2491666/150‑1 (MQ=255)
                                                                               gTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATc                                                                        <  2:3550753/150‑1 (MQ=255)
                                                                                tGCGGGTGCGTGCAGGAAAGGCGCGCTGTCTAACCAGGTGCCCACGGCACCAAAAGGGGGCGCAGATTCTCTATCTGCGCCGTTTCGACTCCCAGCTTGATGAAGTGCTCGCGATAGTGAGGCAGTTCGGCGATCGATTCGCGGATATCg                                                                       <  2:2913700/150‑1 (MQ=255)
                                                                                     gTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAg                                                                  <  1:3275311/150‑1 (MQ=255)
                                                                                      tCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAACAGGGGGCGCAGATACATTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAgc                                                                 <  2:283072/150‑1 (MQ=255)
                                                                                          tGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCa                                                             <  1:2241744/150‑1 (MQ=255)
                                                                                          tGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCa                                                             >  1:3795602/1‑150 (MQ=255)
                                                                                          tGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGAAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCa                                                             >  1:2404530/1‑150 (MQ=255)
                                                                                             aGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGt                                                          >  1:1879036/1‑150 (MQ=255)
                                                                                             aGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGt                                                          >  1:1636538/1‑150 (MQ=255)
                                                                                             aGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGt                                                          >  1:4169993/1‑150 (MQ=255)
                                                                                             aGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGt                                                          >  1:1636526/1‑150 (MQ=255)
                                                                                                aaaGGCGCGCAGTCTAACCAGGTGCCCAAGGCACCAACCGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTggg                                                       <  2:4169993/150‑1 (MQ=255)
                                                                                                    gcgcgcAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTgg                                                        <  2:1129861/145‑1 (MQ=255)
                                                                                                    gcgcgcAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTgg                                                        >  1:1129861/1‑145 (MQ=255)
                                                                                                        gcAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGccc                                               <  2:2809771/150‑1 (MQ=255)
                                                                                                         cAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCct                                              <  2:5182085/150‑1 (MQ=255)
                                                                                                          aGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCctt                                             >  1:2895947/1‑150 (MQ=255)
                                                                                                          aGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCctt                                             >  1:548967/1‑150 (MQ=255)
                                                                                                           gTCTAACCAGTTGCCCAAGGCACCAAACGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttc                                            <  1:4401905/150‑1 (MQ=255)
                                                                                                           gTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAGTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttc                                            <  1:2961284/150‑1 (MQ=255)
                                                                                                           gTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttc                                            <  1:4401836/150‑1 (MQ=255)
                                                                                                              tAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTg                                         <  1:3608014/150‑1 (MQ=255)
                                                                                                              tAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTg                                         >  1:2967669/1‑150 (MQ=255)
                                                                                                              tAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTg                                         >  1:4904858/1‑150 (MQ=255)
                                                                                                              tAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTg                                         >  1:578800/1‑150 (MQ=255)
                                                                                                               aaCCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGa                                        <  2:4401615/150‑1 (MQ=255)
                                                                                                                aCCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGaa                                       >  2:3116046/1‑150 (MQ=255)
                                                                                                                aCCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGaa                                       <  2:2954374/150‑1 (MQ=255)
                                                                                                                  cAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATCCAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCGTCTTGAAGc                                     >  2:4039357/1‑150 (MQ=255)
                                                                                                                  cAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGAGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGc                                     >  2:3445951/1‑150 (MQ=255)
                                                                                                                   aGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCt                                    <  2:3233305/150‑1 (MQ=255)
                                                                                                                   aGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCt                                    <  2:929939/150‑1 (MQ=255)
                                                                                                                    gTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGATTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTg                                   >  1:791026/1‑150 (MQ=255)
                                                                                                                        cccAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTcac                               <  2:548967/150‑1 (MQ=255)
                                                                                                                          cAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTcacgc                             >  1:3458031/1‑150 (MQ=255)
                                                                                                                               cACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTc                        >  2:4652099/1‑150 (MQ=255)
                                                                                                                                aCCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTcc                       <  1:3254949/150‑1 (MQ=255)
                                                                                                                                 ccAAAAGGGGGCGCTGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCATGTGGGTATTTCCCTTCTTGAAGCTGTCACGCACGTCCg                      <  2:782222/150‑1 (MQ=255)
                                                                                                                                 ccAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCg                      <  2:781960/150‑1 (MQ=255)
                                                                                                                                  cAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTcgc                                                                               >  1:3587730/1‑92 (MQ=255)
                                                                                                                                   aaacGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGgc                    <  2:2338291/150‑1 (MQ=255)
                                                                                                                                    aaaGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGgcg                   >  1:2785398/1‑150 (MQ=255)
                                                                                                                                        ggggCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCa               <  2:1476177/150‑1 (MQ=255)
                                                                                                                                          ggCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTCGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAgc             <  1:1142618/150‑1 (MQ=255)
                                                                                                                                           gcgcAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAgcg            >  1:3954232/1‑150 (MQ=255)
                                                                                                                                           gcgcAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAgcg            >  1:3954220/1‑150 (MQ=255)
                                                                                                                                               aGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGcgcgcg        <  1:1002738/150‑1 (MQ=255)
                                                                                                                                                   aCAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCa    >  2:3890577/1‑150 (MQ=255)
                                                                                                                                                   aCAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCa    <  1:124220/150‑1 (MQ=255)
                                                                                                                                                   aCAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCa    <  1:2963625/150‑1 (MQ=255)
                                                                                                                                                     aaTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGt                                                      <  1:3415847/98‑1 (MQ=255)
                                                                                                                                                     aaTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGt                                                      >  2:3415847/1‑98 (MQ=255)
                                                                                                                                                     aaTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGt  >  2:3326781/1‑150 (MQ=255)
                                                                                                                                                     aaTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGt  >  2:2953652/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGT  >  NC_002947/5574026‑5574324

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: