Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,103,779 C→A 9.5% F204L (TTC→TTA) ‡ PP_5353 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,103,7790CA 9.5% 187.0 / 9.4 63F204L (TTC→TTA) ‡PP_5353hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base A (3/3);  total (33/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 = 1.46e-01

CCGCACCGGTGATCAACGCGCCGGCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCG  >  NC_002947/6103631‑6103927
                                                                                                                                                    |                                                                                                                                                    
ccGCACCGGTGATCAACGCGCCGGCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTc                                                                                                                                                      <  3:212678/149‑1 (MQ=255)
  gCACCGGTGATCAACGCGCCGGCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCt                                                                                                                                                    <  4:235461/149‑1 (MQ=255)
  gCACCGGTGATCAACGCGCCGGCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCt                                                                                                                                                    >  1:285346/1‑149 (MQ=255)
   caacGGTGATCACCGCGCCGTCCTGGTTGATCGAGCAGAGCTTCGTGGCGGCAGCCCGTAACCAGTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGTTCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTc                                                                                                                                                   <  2:288625/146‑1 (MQ=255)
       ggTGATCACCGCGCCGTCCGGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGAAACCACTTCCAGTCGCTTCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAagc                                                                                                                                               <  5:320251/149‑1 (MQ=255)
       ggTGATCAACGCGCCGGCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAagc                                                                                                                                               <  7:81884/149‑1 (MQ=255)
       ggTGATCAACGCGCCGGCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAagc                                                                                                                                               <  2:160291/149‑1 (MQ=255)
          gATCAACGCGCCGGCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCATGAGTTCCTCAagcgtg                                                                                                                                            >  8:103789/1‑149 (MQ=255)
            tCAACGCGCCGGCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAagcgtgag                                                                                                                                          >  2:120358/1‑149 (MQ=255)
            tCAACGCGCCGGCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGGAACTCAagcattt                                                                                                                                           <  4:182426/148‑5 (MQ=255)
              aaCGCGCCGGCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGGAACTCAagcatttgc                                                                                                                                         <  4:116512/148‑7 (MQ=255)
                 gcgcCGTCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGc                                                                                                                                     >  2:193894/1‑149 (MQ=255)
                  cgcCGGCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCt                                                                                                                                    >  5:59322/1‑149 (MQ=255)
                     cGGCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAg                                                                                                                                 <  4:57256/149‑1 (MQ=255)
                                aaCGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATgg                                                                                                                      >  2:54683/1‑149 (MQ=255)
                                aaCGAGCAGAGCTTCCTGGCGGCAGCCCATAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATgg                                                                                                                      >  2:68659/1‑149 (MQ=255)
                                 aCGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGc                                                                                                                     >  2:152766/1‑149 (MQ=255)
                                  cGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCt                                                                                                                    >  1:5560/1‑149 (MQ=255)
                                     gCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTc                                                                                                                  >  5:80410/1‑148 (MQ=255)
                                     gCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCCCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAAAGCTTGAGTTCCTCAAGCGTGAGGGTGCTGAGCTGGGTGATGGCTTc                                                                                                                  >  5:95167/1‑148 (MQ=255)
                                        gagCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGt                                                                                                              <  6:238626/149‑1 (MQ=255)
                                           cTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTcta                                                                                                           >  5:119556/1‑149 (MQ=255)
                                              ccTGGCTGCATCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGTTCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTctacc                                                                                                         <  8:52390/148‑1 (MQ=255)
                                               cTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTctaccta                                                                                                       >  3:313281/1‑149 (MQ=255)
                                                  gcggcAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTctacctac                                                                                                      <  4:313281/147‑1 (MQ=255)
                                                           cGTAACCACTTCCAGTCGCTGCAGCAGCACTCGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCt                                                                                           <  1:152766/149‑1 (MQ=255)
                                                            gTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTg                                                                                          <  1:193894/149‑1 (MQ=255)
                                                                 cACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGt                                                                                     >  3:155440/1‑149 (MQ=255)
                                                                 cACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGt                                                                                     >  7:57883/1‑149 (MQ=255)
                                                                  aCTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTg                                                                                    >  3:135577/1‑149 (MQ=255)
                                                                  aCTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTAGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTg                                                                                    >  2:276559/1‑149 (MQ=255)
                                                                    ttCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGATGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGc                                                                                   >  2:126121/1‑148 (MQ=255)
                                                                       cAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAAc                                                                                <  1:120358/148‑1 (MQ=255)
                                                                              tGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATc                                                                                                   <  5:301247/122‑1 (MQ=255)
                                                                              tGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATc                                                                                                   >  6:301247/1‑122 (MQ=255)
                                                                                 agcagcaCTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTAGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTc                                                                     <  1:276559/149‑1 (MQ=255)
                                                                                   cagcaCTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGa                                                                   <  1:68659/149‑1 (MQ=255)
                                                                                      caCTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGa                                                                <  6:119556/149‑1 (MQ=255)
                                                                                                 cAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGAcc                                                     >  7:172909/1‑149 (MQ=255)
                                                                                                 cAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGAcc                                                     >  8:154493/1‑149 (MQ=255)
                                                                                                 cAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGGCGATGACCGCGCCGAcc                                                     >  8:306085/1‑149 (MQ=255)
                                                                                                  aaTGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCg                                                    >  4:114287/1‑149 (MQ=255)
                                                                                                            aCAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGc                                           >  3:251769/1‑148 (MQ=255)
                                                                                                             cAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTa                                         >  7:189427/1‑149 (MQ=255)
                                                                                                              aaGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTAc                                        >  3:236265/1‑149 (MQ=255)
                                                                                                                gATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTAccc                                      >  5:27745/1‑149 (MQ=255)
                                                                                                                           tCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCg                                                                    <  6:163066/108‑1 (MQ=255)
                                                                                                                           tCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCg                                                                    >  5:163066/1‑108 (MQ=255)
                                                                                                                            cGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGAAAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGt                           <  8:174509/148‑1 (MQ=255)
                                                                                                                             gTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTg                          <  1:54683/148‑1 (MQ=255)
                                                                                                                              ttACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGc                        <  6:59322/149‑1 (MQ=255)
                                                                                                                                aCCCCGCACATGCTTGAGTGCCTCCAGCGTGAGCGGGCGGAGCTGCGGGACGGCTGCCAGGCTCCCTACCTCGCTAACGTGGAAGTGCATCTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCAACCCTGACCTTCCGGGCCgt                      <  6:95167/149‑1 (MQ=255)
                                                                                                                                 ccccGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCgt                      <  7:154493/148‑1 (MQ=255)
                                                                                                                                  cccGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGCTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCgtgt                    >  6:296948/1‑149 (MQ=255)
                                                                                                                                      cAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCgaa                <  2:5560/149‑1 (MQ=255)
                                                                                                                                       tcacgCTTGAGGAACTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCg                  >  2:283812/5‑146 (MQ=255)
                                                                                                                                           gCTTGAGGAACTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCgaa                >  1:91793/1‑144 (MQ=255)
                                                                                                                                           gCTTGAGGAACTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCgaa                <  2:91793/144‑1 (MQ=255)
                                                                                                                                           gCTTGAGGAACTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAtt           >  7:43901/1‑149 (MQ=255)
                                                                                                                                               gAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTcgcg       >  3:143098/1‑149 (MQ=255)
                                                                                                                                               gAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTcgc        <  8:189427/148‑1 (MQ=255)
                                                                                                                                               gAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTcg         <  4:236265/147‑1 (MQ=255)
                                                                                                                                                aGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCgcgc      <  7:11129/149‑1 (MQ=255)
                                                                                                                                                 gTTCCTCAAGCGTGAGCGTGCTGAGCTGGGGGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCgcgc      <  7:103789/148‑1 (MQ=255)
                                                                                                                                                    ccTCAAGCGTGAGCGTGGTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTc   <  1:126121/148‑1 (MQ=255)
                                                                                                                                                    ccTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCg  >  4:207557/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGCACCGGTGATCAACGCGCCGGCCTGGTTCAACGAGCAGAGCTTCCTGGCGGCAGCCCGTAACCACTTCCAGTCGCTGCAGCAGCACTGGGATGCCAATGAAATGGACAAGATCAGCGAGTTCGTTACCCCGCAAATGCTTGAGTTCCTCAAGCGTGAGCGTGCTGAGCTGGGTGATGGCTTCCAGTCTACCTACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCG  >  NC_002947/6103631‑6103927

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