Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,610,735 A→T 5.7% G983G (GGT→GGA)  mcoA ← Mn(II) copper oxidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,610,7350AT 5.7% 214.8 / 2.8 70G983G (GGT→GGA) mcoAMn(II) copper oxidase A
Reads supporting (aligned to +/- strand):  ref base A (27/39);  new base T (2/2);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GGGTGACGGCGATGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTG  >  NC_002947/3610590‑3610882
                                                                                                                                                 |                                                                                                                                                   
gggTGACGGCGATGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCTCTGTGCTGCCGTCCGTGCTCACGGTAccc                                                                                                                                                  <  1:220509/149‑1 (MQ=255)
gggTGACGGCGATGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTAccc                                                                                                                                                  >  8:357262/1‑149 (MQ=255)
gggTGACGGCGATGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTAcc                                                                                                                                                   <  2:34136/148‑1 (MQ=255)
 ggTGACGGCGATGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCt                                                                                                                                                 <  2:195948/149‑1 (MQ=255)
 ggTGACGGCGATGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCt                                                                                                                                                 >  7:200383/1‑149 (MQ=255)
   tGACGGCGATGCTGACCGTGGACGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCATgg                                                                                                                                               >  6:231875/1‑149 (MQ=255)
     aCGGCGATGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGcc                                                                                                                                             >  1:145535/1‑149 (MQ=255)
     aCGGCGATGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGcc                                                                                                                                             <  7:169877/149‑1 (MQ=255)
     aCGGCGATGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGcc                                                                                                                                             >  3:140000/1‑149 (MQ=255)
        gCGATGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGa                                                                                                                                          <  4:140000/149‑1 (MQ=255)
         cGATGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAg                                                                                                                                         <  7:202540/149‑1 (MQ=255)
            tGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCgg                                                                                                                                           <  7:271347/144‑1 (MQ=255)
            tGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGGCCCGGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCgg                                                                                                                                           >  8:271347/1‑144 (MQ=255)
                 aCCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCggc                                                                                                                                   <  2:145535/147‑1 (MQ=255)
                  ccGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCggctgg                                                                                                                                >  8:317873/1‑149 (MQ=255)
                    ggggCCGGGGTCACCGAGGCGGCCCCCTTGGTCTCACGCGCCGTGTAATTGAAGCTGGACGTGAACGGCGTGGCGACCTTAGCGGGTGGGGTGTACGTAACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCggcgggct                                                                                                                              <  2:296010/147‑1 (MQ=255)
                    gTGGCCGGGGTCACCGATGCGGCCCCCTTGGTATCACGCGCCTTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCggctggct                                                                                                                              <  4:211304/149‑1 (MQ=255)
                    gTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCggctggct                                                                                                                              <  2:300444/149‑1 (MQ=255)
                    gTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCggctggct                                                                                                                              <  4:216302/149‑1 (MQ=255)
                       gCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAg                                                                                                                           <  8:289689/149‑1 (MQ=255)
                         cGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGc                                                                                                                         <  3:340998/149‑1 (MQ=255)
                               cACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGa                                                                                                                   >  4:265725/1‑149 (MQ=255)
                                 ccGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGc                                                                                                                         <  2:273050/141‑1 (MQ=255)
                                 ccGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGc                                                                                                                         >  1:273050/1‑141 (MQ=255)
                                  cGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCg                                                                                                                >  7:233960/1‑149 (MQ=255)
                                  cGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCg                                                                                                                >  4:68800/1‑149 (MQ=255)
                                  cGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGAGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCg                                                                                                                >  8:87645/1‑149 (MQ=255)
                                   gAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGt                                                                                                               <  5:182155/149‑1 (MQ=255)
                                     gTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGa                                                                                                             <  5:339697/149‑1 (MQ=255)
                                        ggCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCAGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCg                                                                                                          >  8:348156/1‑149 (MQ=255)
                                                  gTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGc                                                                                                 <  5:211137/148‑1 (MQ=255)
                                                  gTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGc                                                                                                 <  5:205909/148‑1 (MQ=255)
                                                   tATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCg                                                                                               >  5:217929/1‑149 (MQ=255)
                                                        cgcgCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCgg                                                                                          <  5:42823/149‑1 (MQ=255)
                                                        cgcgCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGGCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCgg                                                                                          <  5:48595/149‑1 (MQ=255)
                                                          cgcCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGt                                                                                        <  3:366760/149‑1 (MQ=255)
                                                                   aTTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGATACCCTGGCCTGAGTCCggct                                                                                                                                  >  6:308102/1‑98 (MQ=255)
                                                                   aTTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGATACCCTGGCCTGAGTCCggct                                                                                                                                  <  5:308102/98‑1 (MQ=255)
                                                                   aTTGAAGCTGGCCGTGAACGGCGTGGCGAAAGTAGAGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGt                                                                               >  4:288211/1‑149 (MQ=255)
                                                                      gatgCTGGCCTTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTAGGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCTCTGTCGTTg                                                                             <  8:92897/145‑1 (MQ=255)
                                                                      gAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTg                                                                             <  4:299896/148‑1 (MQ=255)
                                                                      gAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTg                                                                             >  3:299896/1‑148 (MQ=255)
                                                                             gCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCa                                                                      <  3:265725/148‑1 (MQ=255)
                                                                                gTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGtt                                                                  <  8:233960/149‑1 (MQ=255)
                                                                                   aaCGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCa                                                                <  8:263090/148‑1 (MQ=255)
                                                                                     cGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGTGTCGCTGTCGTTGGCCAGCACGTTCAggg                                                             <  1:4082/149‑1 (MQ=255)
                                                                                               aCCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGGCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCtt                                                   <  7:348156/149‑1 (MQ=255)
                                                                                                 cGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTTCCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGa                                                          >  5:31892/1‑140 (MQ=255)
                                                                                                 cGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTTCCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGa                                                          >  5:26172/1‑140 (MQ=255)
                                                                                                 cGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTTCCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGa                                                          <  6:26172/140‑1 (MQ=255)
                                                                                                 cGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTTCCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGa                                                          <  6:31892/140‑1 (MQ=255)
                                                                                                 cGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGc                                                 >  5:202337/1‑149 (MQ=255)
                                                                                                         gtggggtgTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGtt                                         <  7:87645/149‑2 (MQ=255)
                                                                                                         gtggggtgTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGCCTTGCCTGCAGTc                                         >  3:227225/1‑149 (MQ=255)
                                                                                                              gtgtACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCtt                                                   <  7:342257/134‑1 (MQ=255)
                                                                                                              gtgtACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCtt                                                   >  8:342257/1‑134 (MQ=255)
                                                                                                              gtgtACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGc                                    >  2:373900/1‑149 (MQ=255)
                                                                                                               tgtACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGCCTTGCCTGCAGTCGTTGCg                                   <  4:227225/149‑1 (MQ=255)
                                                                                                                      cACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGt                            >  1:257579/1‑149 (MQ=255)
                                                                                                                           tgcttcCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCGGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCTGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTg                       <  5:231875/149‑1 (MQ=255)
                                                                                                                            gctgcCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTgt                      >  5:226201/1‑148 (MQ=255)
                                                                                                                                ccgtccgtGCTCTCGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGTCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCa                   <  3:288211/148‑1 (MQ=255)
                                                                                                                                  gtccgtGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGa                <  7:357262/149‑1 (MQ=255)
                                                                                                                                   tccgtGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGa                >  4:220278/1‑148 (MQ=255)
                                                                                                                                   tccgtGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGa                >  2:161127/1‑148 (MQ=255)
                                                                                                                                     cgtGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAg             <  6:202337/149‑1 (MQ=255)
                                                                                                                                         cTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAggtgg         >  5:247011/1‑149 (MQ=255)
                                                                                                                                           cACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGtt       >  3:238398/1‑149 (MQ=255)
                                                                                                                                                tACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTttgttg  <  3:68800/149‑1 (MQ=255)
                                                                                                                                                tACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTttgttg  <  4:238398/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGGTGACGGCGATGCTGACCGTGGCCGGGGTCACCGAGTCGGCCCCCTTGGTATCACGCGCCGTGTAATTGAAGCTGGCCGTGAACGGCGTGGCGACCGTAGCGGGTGGGGTGTACGTCACTGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTG  >  NC_002947/3610590‑3610882

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