Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 351,749 C→T 7.4% Q79* (CAG→TAG) ‡ 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,7490CT 7.4% 240.1 / 6.4 81Q79* (CAG→TAG) ‡hisAphosphoribosylformimino‑5‑aminoimidazole carboxamide ribonucleotide isomerase
Reads supporting (aligned to +/- strand):  ref base C (43/32);  new base T (3/3);  total (46/35)
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.52e-01

ACCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCATCGTCG  >  NC_002947/351606‑351896
                                                                                                                                               |                                                                                                                                                       
acCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATc                                                                                                                                                    <  2:392101/149‑1 (MQ=255)
  ccGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGACCCTGCCGATCCAGATcgg                                                                                                                                                  >  5:71803/1‑149 (MQ=255)
   cGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggc                                                                                                                                                 <  6:84628/149‑1 (MQ=255)
        aGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggcggcg                                                                                                                                             <  1:266848/148‑1 (MQ=255)
          cATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTa                                                                                                                                          >  7:72902/1‑149 (MQ=255)
             ggCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCTGGATcggcag                                                                                                                                               <  8:344540/141‑3 (MQ=255)
             ggCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCTGGATcggcag                                                                                                                                               >  7:344540/1‑139 (MQ=255)
                cgccAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGtt                                                                                                                                 >  6:158273/1‑149 (MQ=255)
                 gccAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTc                                                                                                                                <  8:50788/149‑1 (MQ=255)
                         ggTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGaaa                                                                                                                        >  3:381334/1‑149 (MQ=255)
                           tCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAAc                                                                                                                       >  4:390262/1‑148 (MQ=255)
                                ggTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCg                                                                                                                  >  8:166217/1‑148 (MQ=255)
                                    gCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAaca                                                                                                              <  3:316364/148‑1 (MQ=255)
                                         cgccgccTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTAGCTGGAAACCATCGAACACTACGt                                                                                                        <  7:257104/149‑1 (MQ=255)
                                            cgccTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCaa                                                                                                     <  1:365681/149‑1 (MQ=255)
                                             gccTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAg                                                                                                    >  1:10737/1‑149 (MQ=255)
                                              ccTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGAGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAgg                                                                                                   >  7:96078/1‑149 (MQ=255)
                                                tGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCt                                                                                                 >  8:42162/1‑149 (MQ=255)
                                                 gCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTg                                                                                                <  4:381334/149‑1 (MQ=255)
                                                   aCCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGc                                                                                              <  7:166217/149‑1 (MQ=255)
                                                     cTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGt                                                                                            >  4:120209/1‑149 (MQ=255)
                                                      tGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATc                                                                                                                   <  5:251412/125‑1 (MQ=255)
                                                      tGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATc                                                                                                                   >  6:251412/1‑125 (MQ=255)
                                                           gACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCt                                                                                       <  7:190479/148‑1 (MQ=255)
                                                              cTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGt                                                                                   <  7:42162/149‑1 (MQ=255)
                                                              cTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCCATCCAGATCGGCGGCGGTAGGTCGCGTTAGCATGAAACCATACAACACAACGTAAAAGCTGGCCTCatcaa                                                                                      >  1:131696/1‑145 (MQ=255)
                                                                     gcgcCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGatcatc                                                                            <  8:72902/149‑1 (MQ=255)
                                                                        ccTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGc                                                                         <  2:10737/149‑1 (MQ=255)
                                                                         cTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGtt                                                                                                                                 <  3:396860/92‑1 (MQ=255)
                                                                         cTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGtt                                                                                                                                 >  4:396860/1‑92 (MQ=255)
                                                                                   cAGCCGGTCAATGGTGAAGTTGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGatcatc                                                                            <  6:403496/135‑1 (MQ=255)
                                                                                   cAGCCGGTCAATGGTGAAGTTGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGatcatc                                                                            >  5:403496/1‑135 (MQ=255)
                                                                                   cAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCg                                                                           >  8:194259/1‑136 (MQ=255)
                                                                                   cAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCg                                                                           <  7:194259/136‑1 (MQ=255)
                                                                                   cAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCg                                                               >  2:405408/1‑148 (MQ=255)
                                                                                   cAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGt                                                              >  2:411843/1‑149 (MQ=255)
                                                                                     gCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTAc                                                                                                          <  7:194326/103‑1 (MQ=255)
                                                                                     gCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTAc                                                                                                          >  8:194326/1‑103 (MQ=255)
                                                                                      ccGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCaa                                                           <  1:411843/149‑1 (MQ=255)
                                                                                      ccGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCaa                                                           <  1:405408/149‑1 (MQ=255)
                                                                                            aaTGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAg                                                                   <  3:168267/135‑1 (MQ=255)
                                                                                            aaTGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAg                                                                   >  4:168267/1‑135 (MQ=255)
                                                                                            aaTGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGcc                                                    >  7:158625/1‑149 (MQ=255)
                                                                                                gTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAg                                                >  4:85250/1‑149 (MQ=255)
                                                                                                 tGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGt                                               <  3:120209/149‑1 (MQ=255)
                                                                                                     gTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTtcgt                                           <  5:158273/149‑1 (MQ=255)
                                                                                                         ttACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCt                                       >  5:248484/1‑149 (MQ=255)
                                                                                                         ttACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCt                                       >  3:55526/1‑149 (MQ=255)
                                                                                                               ccATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCg                                 >  4:47544/1‑149 (MQ=255)
                                                                                                               ccATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCg                                 >  3:315386/1‑149 (MQ=255)
                                                                                                                    gCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCaa                            <  4:315386/149‑1 (MQ=255)
                                                                                                                    gCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCa                             >  5:281690/1‑148 (MQ=255)
                                                                                                                       aaGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCCGAGTTCGTCGCTGAAGCGTGCAAAg                          <  1:4723/149‑1 (MQ=255)
                                                                                                                        aGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGcc                        <  8:96078/149‑1 (MQ=255)
                                                                                                                           gcTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCCTGCAAAGCCTTc                     >  8:27381/1‑149 (MQ=255)
                                                                                                                             tACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTccc                   >  6:30252/1‑149 (MQ=255)
                                                                                                                             tACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTccc                   >  2:279743/1‑149 (MQ=255)
                                                                                                                             tACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTccc                   >  2:307731/1‑149 (MQ=255)
                                                                                                                              aCCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGt                                >  1:145959/1‑135 (MQ=255)
                                                                                                                              aCCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGt                                <  2:145959/135‑1 (MQ=255)
                                                                                                                              aCCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTCCAAAGCATTccct                  >  4:76284/1‑148 (MQ=255)
                                                                                                                               cccGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCgg                 >  4:60299/1‑149 (MQ=255)
                                                                                                                               cccGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTAGATGAACCGTGCAAAGCCTTCCAgg                 >  8:367287/1‑149 (MQ=255)
                                                                                                                                ccGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCggg                <  2:131696/149‑1 (MQ=255)
                                                                                                                                ccGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCggg                >  3:323284/1‑149 (MQ=255)
                                                                                                                                  gAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGc               >  8:123591/1‑148 (MQ=255)
                                                                                                                                    aCCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAg            >  6:414732/1‑149 (MQ=255)
                                                                                                                                     ccTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAgg           <  7:367287/149‑1 (MQ=255)
                                                                                                                                     ccTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAgg           <  3:60299/149‑1 (MQ=255)
                                                                                                                                      ccgCCGATCTGGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGcc                        >  1:179635/3‑135 (MQ=255)
                                                                                                                                      ccgCCGATCTGGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGcc                        <  2:179635/133‑1 (MQ=255)
                                                                                                                                      ccgCCGATCTGGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAgg           >  3:87216/3‑148 (MQ=255)
                                                                                                                                      ccgCCGATCTGGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAgg           <  4:87216/146‑1 (MQ=255)
                                                                                                                                      cTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTtcgtcg                                         <  8:17371/118‑1 (MQ=255)
                                                                                                                                      cTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTtcgtcg                                         >  7:17371/1‑118 (MQ=255)
                                                                                                                                      cTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGt          >  5:192457/1‑149 (MQ=255)
                                                                                                                                        gCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCa        >  7:365856/1‑149 (MQ=255)
                                                                                                                                        gCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCa        >  4:265124/1‑149 (MQ=255)
                                                                                                                                        gCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCa        >  4:27436/1‑149 (MQ=255)
                                                                                                                                            gtCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGc                              <  8:275609/122‑1 (MQ=255)
                                                                                                                                            gtCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGc                              >  7:275609/2‑123 (MQ=255)
                                                                                                                                              ccAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCAtcgtcg  >  4:141792/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                       
ACCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTA‑‑TCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCATCGTCG  >  NC_002947/351606‑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: