Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,161,198 C→A 8.5% P23P (CCC→CCA)  PP_2776 → homocysteine S‑methyltransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,161,1980CA 8.5% 180.2 / 6.0 59P23P (CCC→CCA) PP_2776homocysteine S‑methyltransferase family protein
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base A (3/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CCAAAGCCTGGGTCGCCGCCAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGAA  >  NC_002947/3161054‑3161342
                                                                                                                                                  |                                                                                                                                                
ccAAAGCCTGGGTCGCCGCCAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTcc                                                                                                                                              <  1:5648/149‑1 (MQ=255)
 cAAAGCCTGGGTCGCCGCCAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCg                                                                                                                                             >  5:334802/1‑149 (MQ=255)
  aaaGCCTGGGTCGCCGCCAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCgg                                                                                                                                            <  6:43530/149‑1 (MQ=255)
          gggcgccgcCAATCAGGCGACCGTCGATCGTTGGGTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCGCGATGGTGGGATGGGACGAGAACTCACGCGGATCGGCGCACCCTTCCGGCAGCCg                                                                                                                                      <  4:160240/144‑1 (MQ=255)
           gTCGCCGCCAATCAGGCGACCGTCGCTCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGaa                                                                                                                                   <  7:151272/149‑1 (MQ=255)
              gccgccAATCAGGCGTCCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCTGCCGGaagg                                                                                                                                 <  3:4302/148‑3 (MQ=255)
               ccgccAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGt                                                                                                                               <  5:227509/149‑1 (MQ=255)
                cgccAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTc                                                                                                                              <  8:297281/149‑1 (MQ=255)
                 gccAATCAGGCGACCGGCGATCGTTGGTTTTTTTCCACTGGCACCCATGTGGCGTCCCCCAACCCCCCCCCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGt                                                                                                                               <  1:308234/149‑1 (MQ=255)
                    aaTCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCa                                                                                                                          >  4:355268/1‑149 (MQ=255)
                          gCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGgcg                                                                                                                    >  2:99910/1‑149 (MQ=255)
                                   gATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGaa                                                                                                           <  6:334802/149‑1 (MQ=255)
                                       gttggttgAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGATTGGTCGGCACTGGCGCTGTTCGAAgcgc                                                                                                       <  1:99910/149‑1 (MQ=255)
                                       gttggttgAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTg                                                                                                                 <  2:334917/139‑1 (MQ=255)
                                       gttggttgAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTg                                                                                                                 >  1:334917/1‑139 (MQ=255)
                                        ttggttgAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGcc                                                                                                      >  5:357227/1‑149 (MQ=255)
                                             tgggTGCCAGTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGaa                                                                                                  >  8:176546/4‑148 (MQ=255)
                                                   ccACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGa                                                                                           <  7:87414/149‑1 (MQ=255)
                                                        ggCAGCCACCTGGG‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTGGATGGGGGGATGGGGCGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGGCCGAATACGTGACGCAg                                                                                      >  1:269453/1‑149 (MQ=255)
                                                         gCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAgg                                                                                     >  8:226835/1‑149 (MQ=255)
                                                             ccAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCa                                                                                 >  4:340358/1‑149 (MQ=255)
                                                              cAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCAc                                                                                <  3:340358/149‑1 (MQ=255)
                                                                       t‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGt                                                                         >  4:184379/1‑148 (MQ=255)
                                                                       t‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGt                                                                         >  2:182547/1‑148 (MQ=255)
                                                                                  cTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTgccgcc                                                              <  3:355268/149‑1 (MQ=255)
                                                                                   tCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGccgccg                                                             >  7:181834/1‑149 (MQ=255)
                                                                                   tCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGccgccg                                                             >  7:203847/1‑149 (MQ=255)
                                                                                     ggCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCgg                                                            <  7:176546/148‑1 (MQ=255)
                                                                                     ggCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCgg                                                            <  3:266092/148‑1 (MQ=255)
                                                                                      gCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCgggg                                                          >  8:221425/1‑149 (MQ=255)
                                                                                           ccTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGTGCCAgg                                                     <  3:184379/149‑1 (MQ=255)
                                                                                           ccTGGACTTTGCTCGATGGTGGGAGGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAgg                                                     <  2:242178/149‑1 (MQ=255)
                                                                                           ccTGGACGTTGCTCGATGGTGGTATGGGACGAGAACTCAAGCGGAGCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGACGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAgg                                                     <  4:33671/149‑1 (MQ=255)
                                                                                            cTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAggg                                                    >  5:10316/1‑149 (MQ=255)
                                                                                                 cTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATc                                               <  1:182547/149‑1 (MQ=255)
                                                                                                   acgCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATcac                                             <  7:64999/147‑1 (MQ=255)
                                                                                                     gCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATcacc                                            <  7:221425/148‑1 (MQ=255)
                                                                                                      cTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATcacca                                           <  4:351022/148‑1 (MQ=255)
                                                                                                      cTCGATGGTGGGATGGGACAAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATc                                               >  5:257693/1‑144 (MQ=255)
                                                                                                      cTCGATGGTGGGATGGGACAAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATc                                               <  6:257693/144‑1 (MQ=255)
                                                                                                       tCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATcaccac                                          >  4:339069/1‑148 (MQ=255)
                                                                                                              tgggatgggaCGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGc                                   >  1:99608/1‑148 (MQ=255)
                                                                                                                ggagcggaCGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCCGCCCGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACTTTACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGTTACTCACCACAAACAGCTAc                                >  8:174552/6‑149 (MQ=255)
                                                                                                                 gatgggaCGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTa                                                  <  7:263911/130‑1 (MQ=255)
                                                                                                                 gatgggaCGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTa                                                  >  8:263911/1‑130 (MQ=255)
                                                                                                                         gagaACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGgtggtg                       <  6:357227/149‑1 (MQ=255)
                                                                                                                           gaACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGcc                     >  8:10411/1‑149 (MQ=255)
                                                                                                                             aCTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGt                   <  7:10411/149‑1 (MQ=255)
                                                                                                                               tCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATcaccac                                          >  7:214459/1‑124 (MQ=255)
                                                                                                                               tCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATcaccac                                          <  8:214459/124‑1 (MQ=255)
                                                                                                                                   gCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTcaca             >  2:328963/1‑149 (MQ=255)
                                                                                                                                      gATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTg          <  7:226835/149‑1 (MQ=255)
                                                                                                                                          ggagcACCCTTCCGGCATCCGGAATGGTGGGGACTGGCGCTGTTGGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACGATGATCTACGCGGTGGTGCCGTTTCACATTGGCGa      <  7:174552/146‑1 (MQ=255)
                                                                                                                                             cgaCCATTCCGGCTGCCGGAAGGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAgg    >  7:18688/3‑148 (MQ=255)
                                                                                                                                              gaCCATTCCGGCTGCCGGAAGGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGa      >  1:158149/2‑145 (MQ=255)
                                                                                                                                              gaCCATTCCGGCTGCCGGAAGGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGa      >  1:163779/2‑145 (MQ=255)
                                                                                                                                              gaCCATTCCGGCTGCCGGAAGGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGa      <  2:163779/144‑1 (MQ=255)
                                                                                                                                              gaCCATTCCGGCTGCCGGAAGGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGa      <  2:158149/144‑1 (MQ=255)
                                                                                                                                              cACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGaa  >  6:348644/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CCAAAGCCTGGGTCGCCGCCAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGC‑T‑GCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGAA  >  NC_002947/3161054‑3161342

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