Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,574,176 A→G 5.7% intergenic (‑73/+19) 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,1760AG 5.7% 440.8 / 9.0 141intergenic (‑73/+19)yicG/orninner membrane protein/oligoribonuclease
Reads supporting (aligned to +/- strand):  ref base A (66/65);  new base G (2/6);  total (68/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 159 total aligned reads displayed.

ATGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGT  >  NC_002947/5574027‑5574324
                                                                                                                                                     |                                                                                                                                                    
aTGCCGCGACGCCCGGCAGACAGCGTGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      <  2:1232105/150‑1 (MQ=255)
aTGCCGCGACGCCCGGCAGACAGCGTGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      <  2:1780737/150‑1 (MQ=255)
aTGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGACAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      >  1:929939/1‑150 (MQ=255)
aTGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      >  1:3233305/1‑150 (MQ=255)
aTGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      <  2:4500980/150‑1 (MQ=255)
aTGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      >  1:3969229/1‑150 (MQ=255)
aTGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                      <  1:4030981/150‑1 (MQ=255)
  gCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAata                                                                                                                                                    >  2:1442610/1‑150 (MQ=255)
  gCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAata                                                                                                                                                    >  2:4782026/1‑150 (MQ=255)
     gcgACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCt                                                                                                                                                 >  2:2963625/1‑150 (MQ=255)
       gACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTgc                                                                                                                                               >  1:1364118/1‑150 (MQ=255)
         cGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATATTGTATCTgcgc                                                                                                                                             <  2:4561087/150‑1 (MQ=255)
           cccGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATATTGTATCTGCGccc                                                                                                                                           <  2:3121440/150‑2 (MQ=255)
           cccGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATATTGTATCTGCGccc                                                                                                                                           <  2:3121462/150‑2 (MQ=255)
                         cgcCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGc                                                                                                                             <  1:3409524/150‑1 (MQ=255)
                             aGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTtga                                                                                                                         <  1:2069696/150‑1 (MQ=255)
                                 aTGGCTTCTGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGaa                                                                                                                     <  2:4775654/150‑1 (MQ=255)
                                 aTGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGaa                                                                                                                     >  1:3818424/1‑150 (MQ=255)
                                 aTGGCTTCGGCGGTAAGGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAACAGGGGGCGCAGATACCATATCTGCGCCGTTTCGACTCACAGCTTGATGaa                                                                                                                     <  2:236699/150‑1 (MQ=255)
                                  tGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAg                                                                                                                    >  2:1891555/1‑150 (MQ=255)
                                  tGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAg                                                                                                                    <  2:266203/150‑1 (MQ=255)
                                  tGGCTTCGGCGGTAAAGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAg                                                                                                                    >  2:5386953/1‑150 (MQ=255)
                                             gTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATa                                                                                                         <  1:1551886/150‑1 (MQ=255)
                                              tAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAg                                                                                                        <  2:3389900/150‑1 (MQ=255)
                                               aaTGGCGAAAAGGTAGAGCATCAACAACATGGAGCGGGTCCGTGCAGGAAAGGCGCGCAGACTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACACACAGCTTGATGAAGTGCTCGCGATAGt                                                                                                       >  1:103134/1‑150 (MQ=255)
                                                 tGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGa                                                                                                     <  2:5181562/150‑1 (MQ=255)
                                                          ggTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATATTGTATCTgcgc                                                                                                                                             <  2:1460333/101‑1 (MQ=255)
                                                           gTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCgg                                                                                           <  1:3194621/150‑1 (MQ=255)
                                                           gTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCgg                                                                                           <  1:5141940/150‑1 (MQ=255)
                                                           gTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCgg                                                                                           <  1:2506460/150‑1 (MQ=255)
                                                            tAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGc                                                                                          <  2:2131999/150‑1 (MQ=255)
                                                            tAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGc                                                                                          <  1:270070/150‑1 (MQ=255)
                                                             agagCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcg                                                                                         >  1:2913731/1‑150 (MQ=255)
                                                              gagCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcga                                                                                        >  2:1268788/1‑150 (MQ=255)
                                                              gagCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcga                                                                                        <  2:5096595/150‑1 (MQ=255)
                                                                 cATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcgatcg                                                                                     >  1:3725533/1‑150 (MQ=255)
                                                                  aTCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcgatcga                                                                                    >  1:470337/1‑150 (MQ=255)
                                                                   tCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcgatcgat                                                                                   >  1:3060482/1‑150 (MQ=255)
                                                                       caacaTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTcgc                                                                               <  1:3901310/150‑1 (MQ=255)
                                                                         acaTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAACGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCgg                                                                             <  2:297106/150‑1 (MQ=255)
                                                                           aTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGat                                                                           >  2:4307042/1‑150 (MQ=255)
                                                                            tGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGata                                                                          <  2:3526302/150‑1 (MQ=255)
                                                                            tGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGata                                                                          <  1:2491666/150‑1 (MQ=255)
                                                                              gTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATc                                                                        <  2:3550753/150‑1 (MQ=255)
                                                                                 cGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTc                                                                     <  1:831672/150‑1 (MQ=255)
                                                                                    gTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAg                                                                  <  1:3275311/150‑1 (MQ=255)
                                                                                     tCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAACAGGGGGCGCAGATACATTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAgc                                                                 <  2:283072/150‑1 (MQ=255)
                                                                                        gTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGcc                                                              >  2:1808850/1‑150 (MQ=255)
                                                                                        gTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGcc                                                              <  1:4223421/150‑1 (MQ=255)
                                                                                         tGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCa                                                             <  1:2241744/150‑1 (MQ=255)
                                                                                         tGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCa                                                             >  1:3795602/1‑150 (MQ=255)
                                                                                         tGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGAGTCGCGGATATCGTCCAGCGCCa                                                             >  2:3914723/1‑150 (MQ=255)
                                                                                            aGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGt                                                          >  1:4169993/1‑150 (MQ=255)
                                                                                            aGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGt                                                          >  1:1636374/1‑150 (MQ=255)
                                                                                            aGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGt                                                          >  1:1636538/1‑150 (MQ=255)
                                                                                               aaaGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTggg                                                       <  1:1284968/150‑1 (MQ=255)
                                                                                               aaaGGCGCGCAGTCTAACCAGGTGCCCAAGGCACCAACCGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTggg                                                       <  2:4169993/150‑1 (MQ=255)
                                                                                                   gcgcgcAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTgg                                                        >  1:1129861/1‑145 (MQ=255)
                                                                                                   gcgcgcAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTgg                                                        <  2:1129861/145‑1 (MQ=255)
                                                                                                         aGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCctt                                             >  1:548967/1‑150 (MQ=255)
                                                                                                          gTCTAACCAGTTGCCCAAGGCACCAAACGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttc                                            <  1:4401905/150‑1 (MQ=255)
                                                                                                          gTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAGTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttc                                            <  1:2961284/150‑1 (MQ=255)
                                                                                                           tCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTCTTGCCcttct                                           >  2:236058/1‑150 (MQ=255)
                                                                                                             tAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGCTAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTg                                         >  2:5142097/1‑150 (MQ=255)
                                                                                                             tAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTg                                         >  1:2967669/1‑150 (MQ=255)
                                                                                                             tAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTg                                         <  1:3608014/150‑1 (MQ=255)
                                                                                                             tAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTg                                         >  1:4904858/1‑150 (MQ=255)
                                                                                                             tAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGGATTGCCCTTCTTg                                         >  1:2967660/1‑150 (MQ=255)
                                                                                                              aaCCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGa                                        >  1:734351/1‑150 (MQ=255)
                                                                                                               aCCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGaa                                       <  2:2919157/150‑1 (MQ=255)
                                                                                                               aCCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGaa                                       >  2:3116046/1‑150 (MQ=255)
                                                                                                                 cAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATCCAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCGTCTTGAAGc                                     >  2:4039357/1‑150 (MQ=255)
                                                                                                                 cAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGc                                     >  1:5041148/1‑150 (MQ=255)
                                                                                                                 cAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGAGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGc                                     >  2:3445951/1‑150 (MQ=255)
                                                                                                                  aGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCt                                    <  2:929939/150‑1 (MQ=255)
                                                                                                                  aGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCt                                    <  2:3233305/150‑1 (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)
                                                                                                                                 cAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCgg                     >  2:1601452/1‑150 (MQ=255)
                                                                                                                                 cAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCgg                     >  1:1713832/1‑150 (MQ=255)
                                                                                                                                   aaaGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGgcg                   >  1:2785398/1‑150 (MQ=255)
                                                                                                                                     aaacggCGCAGATATTGTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGcc                 >  1:2269807/5‑150 (MQ=255)
                                                                                                                                      gggggCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCa                                                                   >  1:3431072/1‑99 (MQ=255)
                                                                                                                                      gggggCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCa                                                                   <  2:3431072/99‑1 (MQ=255)
                                                                                                                                       ggggCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCa               >  1:2789061/1‑150 (MQ=255)
                                                                                                                                         ggCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTCGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAgc             <  1:1142618/150‑1 (MQ=255)
                                                                                                                                          gcgcAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCcgcg            >  1:3953887/1‑150 (MQ=255)
                                                                                                                                          gcgcAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAgcg            >  1:3954232/1‑150 (MQ=255)
                                                                                                                                              aGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGcgcgcg        <  1:1002738/150‑1 (MQ=255)
                                                                                                                                              aGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGcgcgcg        <  1:4782026/150‑1 (MQ=255)
                                                                                                                                                  aCAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCa    <  2:2784355/150‑1 (MQ=255)
                                                                                                                                                  aCAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCa    >  2:3890577/1‑150 (MQ=255)
                                                                                                                                                  aCAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCa    <  1:124220/150‑1 (MQ=255)
                                                                                                                                                    aaTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGt  >  2:2953652/1‑150 (MQ=255)
                                                                                                                                                    aaTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGt  >  2:3326781/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGT  >  NC_002947/5574027‑5574324

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