Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,457,149 T→C 6.8% R110R (CGA→CGG)  PP_1275 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,457,1490TC 6.8% 170.3 / 2.4 60R110R (CGA→CGG) PP_1275hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base C (2/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

TGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGT  >  NC_002947/1457001‑1457288
                                                                                                                                                     |                                                                                                                                           
tGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTAt                                                                                                                                             >  3:36641/1‑149 (MQ=255)
tGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTAt                                                                                                                                             >  3:110378/1‑149 (MQ=255)
tGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTAt                                                                                                                                             >  2:166018/1‑149 (MQ=255)
  gaggaAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATcg                                                                                                                                           >  1:287453/1‑149 (MQ=255)
   aggaAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATcgc                                                                                                                                          <  1:166018/149‑1 (MQ=255)
    ggaAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCGATTCACCTGGCCCG‑CGGTATcgcg                                                                                                                                         >  7:283217/1‑149 (MQ=255)
      taaCACATCAATTGGCAAACGGACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCgg                                                                                                                                        <  2:248179/147‑1 (MQ=255)
          acaTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATcgcg                                                                                                                                         >  8:277693/1‑143 (MQ=255)
          acaTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATcgcg                                                                                                                                         <  7:277693/143‑1 (MQ=255)
          acaTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGt                                                                                                                                   >  5:170799/1‑149 (MQ=255)
          acaTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGt                                                                                                                                   >  5:364990/1‑149 (MQ=255)
           caTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCGATTCACCTGGCCCG‑CGGTATCGCGGGGCGtt                                                                                                                                  <  8:283217/149‑1 (MQ=255)
           caTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGtt                                                                                                                                  >  8:202667/1‑149 (MQ=255)
               aaTTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAg                                                                                                                              <  4:168945/149‑1 (MQ=255)
                aTTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGATGTAgg                                                                                                                             >  4:215457/1‑149 (MQ=255)
                  tGGCAAACGGACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCGGCAATTCACCTGGCCCG‑CGATACCGCgggccag                                                                                                                                   <  3:282053/141‑5 (MQ=255)
                    gCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGtt                                                                                                                         <  6:364990/149‑1 (MQ=255)
                    gCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGtt                                                                                                                         >  4:50371/1‑149 (MQ=255)
                    gCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCATTGTAGGTGtt                                                                                                                         <  7:31031/149‑1 (MQ=255)
                        aCGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCg                                                                                                                     <  3:6801/149‑1 (MQ=255)
                          gTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATc                                                                                                                                            >  7:340568/1‑124 (MQ=255)
                          gTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATc                                                                                                                                            <  8:340568/124‑1 (MQ=255)
                             ccTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCa                                                                                                                 <  4:36641/148‑1 (MQ=255)
                             ccTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAg                                                                                                                >  3:340438/1‑149 (MQ=255)
                              cTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAgg                                                                                                               >  4:206361/1‑149 (MQ=255)
                               ttGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGt                                                                                                              >  8:134988/1‑149 (MQ=255)
                                  cAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCa                                                                                                           >  3:123796/1‑149 (MQ=255)
                                       gctcggctcgGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTa                                                                                                      >  2:93999/1‑149 (MQ=255)
                                           ggctcggAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGc                                                                                                  >  1:67210/1‑149 (MQ=255)
                                              tcggAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGTCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CCGTATCGCGGGGCGTTGTAGGAGTTCTCGGCCAGGTTCACGGTATCGCCCg                                                                                               >  2:127642/1‑149 (MQ=255)
                                              tcggAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGTAGTTCTAGGCCAGGTTCACGGTATCGCCCg                                                                                               >  8:345115/1‑149 (MQ=255)
                                                 gAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGt                                                                                             <  1:93999/149‑1 (MQ=255)
                                                  aaTTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGt                                                                                             >  4:27435/1‑147 (MQ=255)
                                                  aaTTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGt                                                                                             <  3:27435/147‑1 (MQ=255)
                                                   aTTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCg                                                                                          <  4:110378/149‑1 (MQ=255)
                                                       cAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGc                                                                                                  <  8:111674/137‑1 (MQ=255)
                                                       cAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGc                                                                                                  >  7:111674/1‑137 (MQ=255)
                                                        aGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTAcgc                                                                                     <  4:340438/149‑1 (MQ=255)
                                                                 ggCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGTAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCAcgcg                                                                            >  2:118881/1‑149 (MQ=255)
                                                                    cggcccccGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACgcgc                                                                           <  8:169925/141‑1 (MQ=255)
                                                                          ccGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAg                                                                   <  6:170799/149‑1 (MQ=255)
                                                                              aaCGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCg                                                               <  3:206361/149‑1 (MQ=255)
                                                                               aCGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCgg                                                              >  4:54377/1‑149 (MQ=255)
                                                                                cGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGt                                                             <  7:134988/149‑1 (MQ=255)
                                                                                   aGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAg                                                          >  5:220305/1‑149 (MQ=255)
                                                                                    gcggcgGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGc                                                         <  3:54377/149‑1 (MQ=255)
                                                                                            aaTGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGt                                                                                             <  3:279714/105‑1 (MQ=255)
                                                                                            aaTGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGt                                                                                             >  4:279714/1‑105 (MQ=255)
                                                                                                       gTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGttct                                      <  3:50371/149‑1 (MQ=255)
                                                                                                          ttGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCg                                          >  3:355946/1‑142 (MQ=255)
                                                                                                          ttGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCg                                          <  4:355946/142‑1 (MQ=255)
                                                                                                            gCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCAt                                  <  7:345115/148‑1 (MQ=255)
                                                                                                            gCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCAt                                  <  2:287453/148‑1 (MQ=255)
                                                                                                              cACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCg                               >  3:159276/1‑149 (MQ=255)
                                                                                                               agtCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCa                                                                    <  8:229396/109‑1 (MQ=255)
                                                                                                                        ccGGCTTCAATTCACCTGGCCCG‑CGGTAGCGCGGGGCGGTGTCGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGGCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAg                     <  3:215457/149‑1 (MQ=255)
                                                                                                                             tGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGAGTACCTGCAAGGCACg                >  8:16264/1‑149 (MQ=255)
                                                                                                                              gCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACgg               >  5:107797/1‑149 (MQ=255)
                                                                                                                                aaTTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTc             <  4:123796/149‑1 (MQ=255)
                                                                                                                                       acgccccG‑CGATACCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTcggcgg        <  5:41705/143‑1 (MQ=255)
                                                                                                                                       acgccccG‑CGATACCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGATGCCGTTCTTCATGCGTACCTGCAAGGCACGGTcggcgg        >  6:41705/5‑147 (MQ=255)
                                                                                                                                         gccccG‑CGATACCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGa      >  3:309880/3‑147 (MQ=255)
                                                                                                                                           cccG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAgt  >  5:169172/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
TGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCG‑CGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGT  >  NC_002947/1457001‑1457288

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