Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,497,004 T→A 11.3% L26Q (CTG→CAG)  cobL → precorrin‑6Y C(5,15)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,497,0040TA11.3% 161.0 / 10.7 62L26Q (CTG→CAG) cobLprecorrin‑6Y C(5,15)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base A (3/4);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

TACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGCCGGTGTGCGTGCTG  >  NC_002947/5496861‑5497146
                                                                                                                                                 |                                                                                                                                              
tACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCctgct                                                                                                                                             <  2:58163/147‑1 (MQ=255)
tACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTgg                                                                                                                                           <  5:248216/149‑1 (MQ=255)
 aCCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCTAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTgg                                                                                                                                           <  4:25134/148‑1 (MQ=255)
   cTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGgcg                                                                                                                                        >  5:165891/1‑149 (MQ=255)
        ggAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCCAGCGCCCAGCAGGgcgcgc                                                                                                                                     <  8:8433/149‑3 (MQ=255)
         gAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGcc                                                                                                                                  >  3:37009/1‑149 (MQ=255)
           gCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCgg                                                                                                                                <  6:306218/149‑1 (MQ=255)
           gCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGACCCGGCCAGCGCCCAGCAGGgc                                                                                                                                         >  5:64303/1‑142 (MQ=255)
           gCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGACCCGGCCAGCGCCCAGCAGGgc                                                                                                                                         <  6:64303/142‑1 (MQ=255)
                cgtGCCCGCGAAGGGGCGCGAAGCGGCCCCAAACATGGCCGACTGAAATACCTGACACCCTCGCTGACACTAGTAGGCATCGCCGAAGCCGGCGTCCGCGGCCTGCGCAACCAGCCCCGGC‑‑GCCCCCTGCCGCGCGCTGGCCGGGTCtt                                                                                                                           <  2:127705/147‑1 (MQ=255)
                   gCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCgg                                                                                                                        >  8:110854/1‑149 (MQ=255)
                     ccGTGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGCCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCa                                                                                                                      <  1:98029/149‑1 (MQ=255)
                     ccGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCCAGCGCCCAGCAGGgcgcgccg                                                                                                                                   <  4:126760/138‑5 (MQ=255)
                     ccGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCCAGCGCCCAGCAGGgcgcgccg                                                                                                                                   >  3:126760/1‑134 (MQ=255)
                     ccGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCCAGCGCCCAGCAGGgcgcgcc                                                                                                                                    >  6:145646/1‑134 (MQ=255)
                     ccGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCCAGCGCCCAGCAGGgcgcgcc                                                                                                                                    <  5:145646/137‑4 (MQ=255)
                     ccGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCgg                                                                                                                        >  1:35041/1‑147 (MQ=255)
                     ccGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCgg                                                                                                                        >  4:270916/1‑147 (MQ=255)
                     ccGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCgg                                                                                                                        >  8:269290/1‑147 (MQ=255)
                     ccGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCa                                                                                                                      >  1:224163/1‑149 (MQ=255)
                               gcgcGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCgcgc                                                                                                              >  1:229636/1‑147 (MQ=255)
                                 gcgAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGc                                                                                                          >  3:169558/1‑149 (MQ=255)
                                          ccccAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTtgc                                                                                                 <  6:165891/149‑1 (MQ=255)
                                           cccAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCAGCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTtgct                                                                                                <  7:269290/149‑1 (MQ=255)
                                            ccAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGc                                                                                                                       >  5:251284/1‑125 (MQ=255)
                                            ccAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGc                                                                                                                       <  6:251284/125‑1 (MQ=255)
                                             cAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTtgctgc                                                                                              <  7:110854/149‑1 (MQ=255)
                                              aaaaaTGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTtgctgc                                                                                              >  6:129615/1‑148 (MQ=255)
                                              aaaaaTGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGACGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTtgctgc                                                                                              >  5:118765/1‑148 (MQ=255)
                                                aaaTGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCAc                                                                                           <  2:229636/149‑1 (MQ=255)
                                                     gCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCTGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACTTTGc                                                                                      <  3:61333/149‑1 (MQ=255)
                                                     gCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGc                                                                                      <  5:129615/149‑1 (MQ=255)
                                                        gATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTg                                                                                   <  2:35041/149‑1 (MQ=255)
                                                             aaaTACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGc                                                                              <  2:41253/149‑1 (MQ=255)
                                                              aaTACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGA‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCAAGCTGGCCTTGCTGCCACGTTGCGTGACCGCt                                                                             >  8:19915/1‑149 (MQ=255)
                                                               aTACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCt                                                                             >  6:48087/1‑148 (MQ=255)
                                                                       cacCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTtg                                                                    <  3:89463/149‑1 (MQ=255)
                                                                         cccTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTtggc                                                                  <  3:267607/149‑1 (MQ=255)
                                                                         cccTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTtggc                                                                  <  4:37009/149‑1 (MQ=255)
                                                                          ccTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTtggcc                                                                 >  8:256135/1‑149 (MQ=255)
                                                                          ccTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTtggcc                                                                 >  6:185639/1‑149 (MQ=255)
                                                                             ggCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCgg                                                                                                                        >  4:266363/1‑91 (MQ=255)
                                                                             ggCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCgg                                                                                                                        <  3:266363/91‑1 (MQ=255)
                                                                                     gtagtaGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGccc                                                      <  2:227072/149‑1 (MQ=255)
                                                                                         tagGCATCGACGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTtc                                                  <  6:118765/149‑1 (MQ=255)
                                                                                             cATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGc                                              <  7:256135/149‑1 (MQ=255)
                                                                                              aTCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGccc                                                           <  2:90133/135‑1 (MQ=255)
                                                                                              aTCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGccc                                                           >  1:90133/1‑135 (MQ=255)
                                                                                                 ggCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTgg                                           <  3:270916/148‑1 (MQ=255)
                                                                                                 ggCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGc                                          >  1:40602/1‑149 (MQ=255)
                                                                                                          ggCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCTCACTGCTGGGCGCTGGACGGGTCTTCGGCAGCGCTCGCCAGCTTGCCTTGCTGCCACGTTGCGTGAGCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCt                                 <  7:19915/149‑1 (MQ=255)
                                                                                                            cTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTgg                               >  8:98627/1‑149 (MQ=255)
                                                                                                                aGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCt                           >  1:170384/1‑149 (MQ=255)
                                                                                                                  cGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCCGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTc                         >  8:77705/1‑149 (MQ=255)
                                                                                                                        gggCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCgcggc                    >  2:94980/1‑148 (MQ=255)
                                                                                                                            aaaCAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGcc               <  6:29182/149‑1 (MQ=255)
                                                                                                                            aaaCAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGcc               >  5:275668/1‑149 (MQ=255)
                                                                                                                            aaaCAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGcc               >  2:1791/1‑149 (MQ=255)
                                                                                                                             aaCAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCAAGCCg              <  2:40602/149‑1 (MQ=255)
                                                                                                                              aCAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGCACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGCCgg             >  7:239114/1‑149 (MQ=255)
                                                                                                                                 ggCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGCCGgtgt          >  6:276067/1‑149 (MQ=255)
                                                                                                                                  gCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGCCGGtgtg         >  2:156646/1‑149 (MQ=255)
                                                                                                                                           gcgcCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGCCGGTGTGCGTGCTg  <  1:1791/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGCCGGTGTGCGTGCTG  >  NC_002947/5496861‑5497146

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: ‑