Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 351,750 A→G 6.9% Q79R (CAG→CGG) ‡ hisA → phosphoribosylformimino‑5‑aminoimidazole carboxamide ribonucleotide isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947351,7500AG 6.9% 214.1 / 4.5 72Q79R (CAG→CGG) ‡hisAphosphoribosylformimino‑5‑aminoimidazole carboxamide ribonucleotide isomerase
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base G (2/3);  total (35/37)
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.86e-01

CGACCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCATCGTCG  >  NC_002947/351604‑351896
                                                                                                                                                   |                                                                                                                                                  
cgacCCGGAGAGCAGGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGGGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGa                                                                                                                                                  <  6:354672/149‑1 (MQ=255)
    ccGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCTGGATcgg                                                                                                                                              <  4:133018/149‑1 (MQ=255)
      ggTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCTGGATcggca                                                                                                                                            <  1:114403/149‑2 (MQ=255)
       gTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggcgg                                                                                                                                           >  6:10314/1‑149 (MQ=255)
                gccgccAAGTTGGGCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATcc                                                                                                                                   <  8:320772/148‑1 (MQ=255)
                 ccgccAAGTGGGTCGAGGGTGGCTGCCGCGGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGt                                                                                                                                 <  5:150185/149‑1 (MQ=255)
                 ccgccAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCCGCCGATCTGGATcggcag                                                                                                                                           >  6:97678/1‑137 (MQ=255)
                 ccgccAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCCGCCGATCTGGATcggcag                                                                                                                                           <  5:97678/139‑3 (MQ=255)
                          gggTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGCCGGCGGTATCCGTTCGCTGGa                                                                                                                         >  6:244033/1‑148 (MQ=255)
                           ggTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGaaa                                                                                                                       <  7:94388/149‑1 (MQ=255)
                             tCGAGGTTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAAc                                                                                                                      >  2:67146/1‑148 (MQ=255)
                              cGAGGGTGGCTGCCGCCGTCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCa                                                                                                                    <  8:341581/149‑1 (MQ=255)
                               gAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAATGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGc                                                                                                                             <  5:230708/139‑1 (MQ=255)
                               gAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAATGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGc                                                                                                                             >  6:230708/1‑139 (MQ=255)
                                aGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATc                                                                                                                  <  7:242453/149‑1 (MQ=255)
                                aGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATc                                                                                                                  <  2:280969/149‑1 (MQ=255)
                                 gggTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCg                                                                                                                 <  2:103682/149‑1 (MQ=255)
                                 gggTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCg                                                                                                                 >  3:60696/1‑149 (MQ=255)
                                      gCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCTAGCGCTACCCGAACCTGCCTGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAac                                                                                                              <  1:77979/148‑1 (MQ=255)
                                        tGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTa                                                                                                          >  2:94795/1‑149 (MQ=255)
                                          ccgccgccTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACg                                                                                                        <  1:67146/149‑1 (MQ=255)
                                            gccgccTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTc                                                                                                      <  5:41674/149‑1 (MQ=255)
                                            gccgccTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTc                                                                                                      <  3:346763/149‑1 (MQ=255)
                                                   gCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTg                                                                                               >  3:264277/1‑149 (MQ=255)
                                                      ccTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCg                                                                                            >  2:265123/1‑149 (MQ=255)
                                                       cTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGt                                                                                           <  5:10314/149‑1 (MQ=255)
                                                         ggTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGa                                                                                                                         >  2:105945/1‑117 (MQ=255)
                                                           tCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGc                                                                                       >  7:333338/1‑149 (MQ=255)
                                                               ccTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACg                                                                                   <  1:52470/149‑1 (MQ=255)
                                                               ccTCAACGGCGCCTTCGAAGGACAGCCGGTCAATGGTGAAGTGCATACCGACATCGCCAACCGCTACCCGAACCTCCCGATACACATCGGAGGCGTTATCCGTTCGCTGGAAACCATCGAACACTACGTCAACGCTGGCGTCAGCTACg                                                                                   >  4:342540/1‑149 (MQ=255)
                                                                  cAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGa                                                                                >  8:170322/1‑149 (MQ=255)
                                                                         gcCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCgg                                                                         >  8:293062/1‑149 (MQ=255)
                                                                               gAAGGCCAGCCGGTCTATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCTGCGGCGGTATCCTTTCTCTGGAAACCATCGAACACTACGTCAAGGCTGGGGTCAGCTACGTTATCATCGGCACCaa                                                                   <  3:342540/149‑1 (MQ=255)
                                                                                aaGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAGCCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAg                                                                  >  5:388755/1‑149 (MQ=255)
                                                                                aaGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAg                                                                  <  1:265123/149‑1 (MQ=255)
                                                                                    ccAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCg                                                              <  6:215055/149‑1 (MQ=255)
                                                                                    ccAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCg                                                              >  2:256142/1‑149 (MQ=255)
                                                                                     cAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGGCGt                                                             <  2:171337/149‑1 (MQ=255)
                                                                                     cAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACATGCCGATCCAGAACGCCGGCGGTAACCGTACGATGGAAACCATCGACCACTACTACAAGCAGGGCGACAGGTACGTGATCATCGGCACCAAGCCCg                                                              >  2:91085/1‑148 (MQ=255)
                                                                                      aGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCTAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTc                                                            <  1:94795/149‑1 (MQ=255)
                                                                                       cccGGTCAATGGTGAAGTTGGTACCGCCAACGCCAATCGCTAACCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCa                                                           <  1:91085/148‑1 (MQ=255)
                                                                                            tCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAagcag                                                      >  3:301946/1‑149 (MQ=255)
                                                                                              aaTGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGcc                                                    <  4:264277/149‑1 (MQ=255)
                                                                                                       gTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTtcgt                                           >  6:61198/1‑149 (MQ=255)
                                                                                                           ttACCGCCATCGCCAAGCACTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCt                                       >  7:331255/1‑149 (MQ=255)
                                                                                                               cGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAg                                   >  4:269582/1‑149 (MQ=255)
                                                                                                                      gCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGc                              >  7:37807/1‑147 (MQ=255)
                                                                                                                      gCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGc                              <  8:37807/147‑1 (MQ=255)
                                                                                                                      gCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCa                             >  6:378852/1‑148 (MQ=255)
                                                                                                                       ccAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCaaa                           <  5:61198/149‑1 (MQ=255)
                                                                                                                       ccAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCa                             >  3:371474/1‑147 (MQ=255)
                                                                                                                        cAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAg                          >  2:233446/1‑149 (MQ=255)
                                                                                                                        cAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAg                          >  6:234259/1‑149 (MQ=255)
                                                                                                                        cAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAg                          >  8:263288/1‑149 (MQ=255)
                                                                                                                          aGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGAACATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGcc                        >  2:221503/1‑149 (MQ=255)
                                                                                                                              cTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCt                                       >  8:21705/1‑130 (MQ=255)
                                                                                                                              cTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCt                                       <  7:21705/130‑1 (MQ=255)
                                                                                                                               tACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAagc                                                        >  1:211038/1‑112 (MQ=255)
                                                                                                                               tACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAagc                                                        <  2:211038/112‑1 (MQ=255)
                                                                                                                               tACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTccc                   >  7:111121/1‑149 (MQ=255)
                                                                                                                                 cccGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTccc                   <  4:301946/147‑1 (MQ=255)
                                                                                                                                  ccGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCggg                <  7:293062/149‑1 (MQ=255)
                                                                                                                                     aaCCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAACGCCTTcagggtc               >  6:265215/1‑141 (MQ=255)
                                                                                                                                      aCCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAg            <  8:111121/149‑1 (MQ=255)
                                                                                                                                       gccgCCGATCTGGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAgg           >  1:95367/4‑149 (MQ=255)
                                                                                                                                           ccGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTc                                                            >  2:369347/1‑96 (MQ=255)
                                                                                                                                           ccGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTc                                                            <  1:369347/96‑1 (MQ=255)
                                                                                                                                            cGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAGCACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCAtc      <  4:60696/149‑1 (MQ=255)
                                                                                                                                                tCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTtcgtc                                          >  6:105275/1‑110 (MQ=255)
                                                                                                                                                tCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTtcgtc                                          <  5:105275/110‑1 (MQ=255)
                                                                                                                                                tCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCAtcgtc   <  5:234259/149‑1 (MQ=255)
                                                                                                                                                  cAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCAtcgtcg  <  7:263288/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGACCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCATCGTCG  >  NC_002947/351604‑351896

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: