Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,869,271 C→T 6.1% T33M (ACG→ATG)  PP_3418 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,869,2710CT 6.1% 193.8 / 2.1 65T33M (ACG→ATG) PP_3418hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base T (2/2);  total (34/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.57e-01

GTAAGGCCAAAGGGAAGTATGAGCCTTACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAATGACAACCGCGATCTGCCGGGTATGGAGGGCAGT  >  NC_002947/3869123‑3869419
                                                                                                                                                    |                                                                                                                                                    
gTAAGGCCAAAGGGAAGTATGAGCCTTACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATAc                                                                                                                                                      >  7:213640/1‑149 (MQ=255)
 gaaGGCCAAAGCGAAGTAGGAGCCTTCCCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGGCGATTGCCCTGCTCGCCTCGCCGCTGGCCTGTGCTGTCGGCACAGGGCCGACCTACCCCGACGACCCGCACAATACg                                                                                                                                                     <  4:270340/148‑1 (MQ=255)
      ccAAAGGGAAGTATGAGCCTTACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGgcgc                                                                                                                                                 >  3:139949/1‑148 (MQ=255)
      ccAAAGGGAAGTATGAGCCTTACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGcc                                                                                                                                                <  6:318214/149‑1 (MQ=255)
       cAAAGGGAAGTATGAGCCTTACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCca                                                                                                                                               >  1:327472/1‑149 (MQ=255)
       cAAAGGGAAGTATGAGCCTTACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCca                                                                                                                                               <  6:282651/149‑1 (MQ=255)
       cAAAGGGAAGTATGAGCCTTACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCca                                                                                                                                               >  5:282651/1‑149 (MQ=255)
       cAAAGGGAAGTATGAGCCTTACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCca                                                                                                                                               >  5:367604/1‑149 (MQ=255)
         aaGGGAAGTATGAGCCTTACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCcac                                                                                                                                              >  2:80827/1‑148 (MQ=255)
              aaGTATGAGCCTTACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCg                                                                                                                                        >  5:351795/1‑149 (MQ=255)
                   tGAGCCTTACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTg                                                                                                                                   >  3:4105/1‑149 (MQ=255)
                       ccTTACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGAGCCGCACAATACGGCGCCACAACCCGGCGAGGa                                                                                                                                 >  6:124868/1‑147 (MQ=255)
                         ttACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGc                                                                                                                             >  7:326912/1‑149 (MQ=255)
                            cccGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCAcc                                                                                                                          <  6:85679/149‑1 (MQ=255)
                            cccGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCAcc                                                                                                                          <  4:157654/149‑1 (MQ=255)
                             ccGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCAccc                                                                                                                         <  6:367604/149‑1 (MQ=255)
                              cGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCt                                                                                                                        >  3:69149/1‑149 (MQ=255)
                               gTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGgcgc                                                                                                                                                 >  7:149890/1‑123 (MQ=255)
                               gTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGgcgc                                                                                                                                                 <  8:149890/123‑1 (MQ=255)
                                  cAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAAc                                                                                                                    <  8:213640/149‑1 (MQ=255)
                                   aGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAAcc                                                                                                                   <  3:354041/149‑1 (MQ=255)
                                      aaaGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGa                                                                                                                <  4:139949/149‑1 (MQ=255)
                                             ccTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGaa                                                                                                         <  2:327472/149‑1 (MQ=255)
                                                 ccATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGcc                                                                                                      <  4:4105/148‑1 (MQ=255)
                                                          gCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCgg                                                                                            <  8:326912/149‑1 (MQ=255)
                                                          gCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCgg                                                                                            <  1:80827/149‑1 (MQ=255)
                                                                 cctcctGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCg                                                                                     <  3:299368/149‑1 (MQ=255)
                                                                 cctcctGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCg                                                                                     >  3:279453/1‑149 (MQ=255)
                                                                     ctGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGAc                                                                                                                                >  7:199596/1‑102 (MQ=255)
                                                                     ctGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGAc                                                                                                                                <  8:199596/102‑1 (MQ=255)
                                                                      tGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGACGCGGGCCACCGACCcg                                                                                >  7:156825/1‑149 (MQ=255)
                                                                          gATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCgcgc                                                                             <  6:351795/148‑1 (MQ=255)
                                                                            ttGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATa                                                                          <  7:251001/149‑1 (MQ=255)
                                                                                cgtGCTCGCCTCGCCGCTGGCGTTTGCTGTAGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATGGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAgg                                                                      <  8:156825/147‑1 (MQ=255)
                                                                                ccTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAgg                                                                      >  7:199510/1‑149 (MQ=255)
                                                                                ccTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAgg                                                                      >  3:194611/1‑149 (MQ=255)
                                                                                 cTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAggg                                                                     >  7:4287/1‑149 (MQ=255)
                                                                                  tGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAggg                                                                     >  7:244045/1‑148 (MQ=255)
                                                                                     tcgcctcgccGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAAc                                                                 >  3:360318/1‑149 (MQ=255)
                                                                                      cgcctcgccGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGCCCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCCTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACg                                                                >  2:363944/1‑149 (MQ=255)
                                                                                       gcctcgccGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGCCAGCACCCTCAACCCGATCGAGAAGCCCTTACCGCGGGCCACCGACCCGCGCATACAGGGCAACGa                                                               >  5:368872/1‑149 (MQ=255)
                                                                                         ctcgccGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAAc                                                             >  8:180997/1‑149 (MQ=255)
                                                                                                   gCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGccc                                                   <  3:184151/149‑1 (MQ=255)
                                                                                                   gCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGccc                                                   <  4:279453/149‑1 (MQ=255)
                                                                                                        tGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCgcc                                                <  4:69149/147‑1 (MQ=255)
                                                                                                          cTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACACGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCgccagc                                             >  6:81129/1‑148 (MQ=255)
                                                                                                           tGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCgccagcc                                            >  1:46931/1‑148 (MQ=255)
                                                                                                               ggCACAGGTCCGACCTACCCCGACGACCCGGACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGTGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACa                                                                        >  1:35079/1‑116 (MQ=255)
                                                                                                               ggCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGTGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACa                                                                        <  2:35079/116‑1 (MQ=255)
                                                                                                               ggCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCa                                       <  8:4287/149‑1 (MQ=255)
                                                                                                                    aGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAAt                                   >  4:106201/1‑148 (MQ=255)
                                                                                                                    aGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAATg                                  <  7:180997/149‑1 (MQ=255)
                                                                                                                            ccTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACTAACCGCAAAGCCCGCCAGCCAAGCAGAATTACAACcgc                          <  8:244045/149‑1 (MQ=255)
                                                                                                                            ccTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCAAACAGGGCAACGAACCGCAAAGCCCGCCACCAAACCAGAATGACAACcgc                          >  8:217565/1‑149 (MQ=255)
                                                                                                                               aCCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAATGACAACCGCGAt                       >  1:149362/1‑149 (MQ=255)
                                                                                                                                 cccGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAgg                                                                      <  4:250773/100‑1 (MQ=255)
                                                                                                                                 cccGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAgg                                                                      >  3:250773/1‑100 (MQ=255)
                                                                                                                                  ccGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAATGACAACCGCGATCTg                    >  6:10857/1‑149 (MQ=255)
                                                                                                                                          ccGCACAATATGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAATGACAACCGCGAt                       >  5:160024/1‑138 (MQ=255)
                                                                                                                                          ccGCACAATATGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAATGACAACCGCGAt                       >  5:165756/1‑138 (MQ=255)
                                                                                                                                          ccGCACAATATGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAATGACAACCGCGAt                       <  6:160024/138‑1 (MQ=255)
                                                                                                                                          ccGCACAATATGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAATGACAACCGCGAt                       <  6:165756/138‑1 (MQ=255)
                                                                                                                                          ccGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGa                                                                                                                <  4:337251/49‑1 (MQ=255)
                                                                                                                                          ccGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGa                                                                                                                >  3:337251/1‑49 (MQ=255)
                                                                                                                                               caATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCCTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAATGACAACCGCGATCTGCCGGGTATGGAgg       <  1:363944/149‑1 (MQ=255)
                                                                                                                                                aaTACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAATGACAACCGCGATCTGCCGGGTATGGAggg      >  4:99398/1‑149 (MQ=255)
                                                                                                                                                    cGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAATGACAACCGCGATCTGCCGGGTATGGAGGGCAGt  <  5:10857/149‑1 (MQ=255)
                                                                                                                                                    cGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAATGACAACCGCGATCTGCCGGGTATGGAGGGCAGt  <  3:326370/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTAAGGCCAAAGGGAAGTATGAGCCTTACCCGTACAGCAAAGGAGCCTGCCATGCCCCGCTCAACCCTCCTGTCGATTGCCCTGCTCGCCTCGCCGCTGGCCTTTGCTGTCGGCACAGGTCCGACCTACCCCGACGACCCGCACAATACGGCGCCACAACCCGGCGTGGACAGCACCCTCAACCCGATCGAGAAGCCCTTGCCGCGGGCCACCGACCCGCGCATACAGGGCAACGAACCGCAAAGCCCGCCAGCCAAGCAGAATGACAACCGCGATCTGCCGGGTATGGAGGGCAGT  >  NC_002947/3869123‑3869419

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: