Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 450,923 C→A 7.0% noncoding (117/127 nt) PP_mr05 → RsmY ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947450,9230CA 7.0% 678.1 / 28.5 200noncoding (117/127 nt)PP_mr05RsmY ncRNA
Reads supporting (aligned to +/- strand):  ref base C (68/118);  new base A (3/11);  total (71/129)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 213 total aligned reads displayed.

TCTAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCAGGC  >  NC_002947/450775‑451070
                                                                                                                                                    |                                                                                                                                                   
tCTAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcg                                                                                                                                                    <  2:2866868/150‑1 (MQ=255)
  tAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggc                                                                                                                                                  >  2:3023302/1‑150 (MQ=255)
  tAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggc                                                                                                                                                  >  2:599946/1‑150 (MQ=255)
   aaGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggcg                                                                                                                                                 >  2:1513693/1‑150 (MQ=255)
   aaGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggcg                                                                                                                                                 <  2:2856129/150‑1 (MQ=255)
   aaGATAGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggcg                                                                                                                                                 >  2:4133245/1‑150 (MQ=255)
    aGATTGGATCCACTGTGAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggcgg                                                                                                                                                <  2:3036895/150‑1 (MQ=255)
     gATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCggg                                                                                                                                               <  2:1544472/150‑1 (MQ=255)
         ggATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGttt                                                                                                                                           <  2:2080252/150‑1 (MQ=255)
         ggATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGttt                                                                                                                                           <  1:1879205/150‑1 (MQ=255)
         ggATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGttt                                                                                                                                           <  2:2727009/150‑1 (MQ=255)
         ggATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTAGAAAACCCCGCCGAAGCGGGGttt                                                                                                                                           <  1:793366/150‑1 (MQ=255)
         ggATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTAGAAAACCCCGCCGAAGCGGGGttt                                                                                                                                           <  1:3062848/150‑1 (MQ=255)
          gATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGtttt                                                                                                                                          <  2:4471317/150‑1 (MQ=255)
           aTCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGTTTTc                                                                                                                                         <  2:174219/150‑1 (MQ=255)
             ccACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCtt                                                                                                                                       <  1:1750350/150‑1 (MQ=255)
              cACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCttt                                                                                                                                      <  2:2899716/150‑1 (MQ=255)
               aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtg                                                                                                                                     <  1:1040228/150‑1 (MQ=255)
               aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTAGAAAACCCCGCCGAAGCGGGGTTTTccaga                                                                                                                                     <  2:4063268/150‑5 (MQ=255)
               aCTGTTAATCTTTACCACATGGACGTAGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACTAAGGACCACCTGACAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCACGCTTCGGCGGGGTTTTCTTtg                                                                                                                                     >  1:2785603/1‑150 (MQ=255)
                 tGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgcg                                                                                                                                   >  2:1239064/1‑150 (MQ=255)
                  gtgaaTCCTTACCACAGGGACGTCGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGACACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgt                                                                                                                                  <  2:829137/147‑1 (MQ=255)
                   ttAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAATCCCCGCTTCGGCGGGGTTTTCTTtgcgtg                                                                                                                                 >  2:3028383/1‑150 (MQ=255)
                    tAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgtgc                                                                                                                                <  1:2022657/150‑1 (MQ=255)
                      aTCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTgcgtgcgt                                                                                                                              <  1:446373/150‑1 (MQ=255)
                          ttACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGttt                                                                                                                          >  1:2364217/1‑150 (MQ=255)
                              cacaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                                      <  2:1868473/150‑1 (MQ=255)
                              cacaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                                      <  1:376267/150‑1 (MQ=255)
                              cacaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                                      <  1:3152109/150‑1 (MQ=255)
                               acaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCCAGGACGACGCCGAAGGACCCCCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAg                                                                                                                     <  1:381715/150‑1 (MQ=255)
                               acaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAg                                                                                                                     <  2:1460574/150‑1 (MQ=255)
                                caTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAgg                                                                                                                    <  1:3780439/150‑1 (MQ=255)
                                   ggACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCa                                                                                                                 <  2:2431072/150‑1 (MQ=255)
                                    gACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAg                                                                                                                <  1:522717/150‑1 (MQ=255)
                                    gACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAg                                                                                                                <  2:4391815/150‑1 (MQ=255)
                                          gcgcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATc                                                                                                          <  1:977785/150‑1 (MQ=255)
                                            gcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCa                                                                                                        >  1:2052579/1‑150 (MQ=255)
                                            gcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCa                                                                                                        >  1:3417348/1‑150 (MQ=255)
                                                   gcgcAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGa                                                                                                 <  2:3959539/150‑1 (MQ=255)
                                                    cgcAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGTGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGAt                                                                                                <  2:1995509/150‑1 (MQ=255)
                                                       aaaGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCt                                                                                             <  2:2876723/150‑1 (MQ=255)
                                                             cAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATc                                                                                       <  1:4251140/150‑1 (MQ=255)
                                                             cAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATc                                                                                       <  2:3108584/150‑1 (MQ=255)
                                                              aGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGTCAGCACATCCAGCGGTGATGCTCCCATCa                                                                                      <  1:4075161/150‑1 (MQ=255)
                                                                   aaGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGa                                                                                 <  1:4170062/150‑1 (MQ=255)
                                                                    aGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAg                                                                                >  2:114127/1‑150 (MQ=255)
                                                                       acgacgacgAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGcgt                                                                             <  2:1751689/150‑1 (MQ=255)
                                                                         gacgacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGcgtcc                                                                           >  1:3960770/1‑150 (MQ=255)
                                                                            gacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGcgtccgtc                                                                        <  1:65000/150‑1 (MQ=255)
                                                                            gacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGcgtccgtc                                                                        >  1:2228138/1‑150 (MQ=255)
                                                                            gacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGcgtccgtc                                                                        <  1:3517892/150‑1 (MQ=255)
                                                                              cgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGtctc                                                                      <  2:2254457/150‑1 (MQ=255)
                                                                               gaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTctct                                                                     >  1:2455883/1‑150 (MQ=255)
                                                                               gaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTGTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTctct                                                                     >  2:3432110/1‑150 (MQ=255)
                                                                               gaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGCTGTTCTTTGCGTGCGTGTTTTACAGGCAAGCACATCCAGAGGAGATGCTCCCATCAGGCGAGCGTCAGTctct                                                                     >  2:4472991/1‑150 (MQ=255)
                                                                                 aGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtc                                                                   <  1:4040739/150‑1 (MQ=255)
                                                                                   gACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtcct                                                                 <  1:142496/150‑1 (MQ=255)
                                                                                      caccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCa                                                              >  1:1113714/1‑150 (MQ=255)
                                                                                      caccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCa                                                              <  1:3215556/150‑1 (MQ=255)
                                                                                          tGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCaa                                                          <  1:1440708/150‑1 (MQ=255)
                                                                                           gCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAAc                                                         >  1:3470852/1‑150 (MQ=255)
                                                                                              aGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTg                                                      <  1:2297528/150‑1 (MQ=255)
                                                                                                gATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTgcg                                                    <  1:1492106/150‑1 (MQ=255)
                                                                                                gATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTgcg                                                    <  1:3621512/150‑1 (MQ=255)
                                                                                                     cGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgcagca                                               <  1:4348026/150‑1 (MQ=255)
                                                                                                     cGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgcagca                                               <  1:4475404/150‑1 (MQ=255)
                                                                                                      ggggCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCttt                                                                                                                                      <  2:854251/62‑1 (MQ=255)
                                                                                                      ggggCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCttt                                                                                                                                      >  1:854251/1‑62 (MQ=255)
                                                                                                      ggggCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCCGCAACGTGCGccgcag                                              >  2:2568149/1‑150 (MQ=255)
                                                                                                      ggggCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGcagcag                                              >  1:1184457/1‑150 (MQ=255)
                                                                                                       gggCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCagcagc                                             <  2:3457762/150‑1 (MQ=255)
                                                                                                             tACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCcaca                                       <  1:1552068/150‑1 (MQ=255)
                                                                                                              aCGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACag                                      >  1:258071/1‑150 (MQ=255)
                                                                                                              aCGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACag                                      >  2:3785684/1‑150 (MQ=255)
                                                                                                              aCGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACag                                      >  1:257962/1‑150 (MQ=255)
                                                                                                                ggATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACagag                                    <  2:1058509/150‑1 (MQ=255)
                                                                                                                 gATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAgg                                   >  1:233233/1‑150 (MQ=255)
                                                                                                                  aTGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGc                                  <  1:3851869/150‑1 (MQ=255)
                                                                                                                  aTGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGc                                  <  1:2684840/150‑1 (MQ=255)
                                                                                                                     tcgcaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctg                               <  2:2194744/147‑1 (MQ=255)
                                                                                                                      cacaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctgc                              >  2:730311/1‑150 (MQ=255)
                                                                                                                       acaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctgct                             >  2:2172247/1‑150 (MQ=255)
                                                                                                                         aGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTgg                           <  1:1827560/150‑1 (MQ=255)
                                                                                                                          gggAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGgc                          >  2:1526801/1‑150 (MQ=255)
                                                                                                                           ggAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGgcg                         >  1:2408539/1‑150 (MQ=255)
                                                                                                                           ggAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGgcg                         <  1:732305/150‑1 (MQ=255)
                                                                                                                            gAAACAGTCTGGAACACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCCGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCACCGTGCGCAGCAGCTCCCCAGAGGCctgctgttcg                        <  2:4104288/150‑5 (MQ=255)
                                                                                                                            gAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGgcgc                        >  2:3941643/1‑150 (MQ=255)
                                                                                                                               aCAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTgg                     <  1:579046/150‑1 (MQ=255)
                                                                                                                               aCAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTgg                     <  2:1321505/150‑1 (MQ=255)
                                                                                                                                   tCTGGAAATCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCAt                 >  2:283350/1‑150 (MQ=255)
                                                                                                                                   tCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCAt                 >  1:68735/1‑150 (MQ=255)
                                                                                                                                    cTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATc                <  2:694389/150‑1 (MQ=255)
                                                                                                                                     tGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcg               >  2:3577928/1‑150 (MQ=255)
                                                                                                                                     tGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcg               <  2:2364217/150‑1 (MQ=255)
                                                                                                                                        aaaaCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGa            <  2:3155729/150‑1 (MQ=255)
                                                                                                                                          aaCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGaaa          <  2:1387360/150‑1 (MQ=255)
                                                                                                                                             cccGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATAc       <  1:2649994/150‑1 (MQ=255)
                                                                                                                                              ccGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATAcc      >  1:2848079/1‑150 (MQ=255)
                                                                                                                                               cGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGt                                                       >  2:4221811/1‑100 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCTAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCAGGC  >  NC_002947/450775‑451070

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