Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 351,749 C→T 6.3% 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 6.3% 286.8 / 5.3 95Q79* (CAG→TAG) ‡hisAphosphoribosylformimino‑5‑aminoimidazole carboxamide ribonucleotide isomerase
Reads supporting (aligned to +/- strand):  ref base C (44/45);  new base T (3/3);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

CGACGACCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCATCGTCG  >  NC_002947/351601‑351896
                                                                                                                                                     |                                                                                                                                                   
cgacgaCCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGGCAAGCGCTACCCGAACCTGCCGATcc                                                                                                                                                     <  1:232860/149‑1 (MQ=255)
cgacgaCCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATcc                                                                                                                                                     <  1:164063/149‑1 (MQ=255)
    gacCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTGCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGACCCTGCCGCTACAgct                                                                                                                                                 >  2:426519/1‑147 (MQ=255)
         ggTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggcg                                                                                                                                            >  7:129869/1‑149 (MQ=255)
          gTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCTGGATcggcag                                                                                                                                           <  3:427265/149‑3 (MQ=255)
          gTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggcgg                                                                                                                                           <  1:13972/149‑1 (MQ=255)
           tGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggcggc                                                                                                                                          <  6:283799/149‑1 (MQ=255)
             aGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggcggcg                                                                                                                                         <  5:303761/148‑1 (MQ=255)
               cATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTa                                                                                                                                      >  7:46739/1‑149 (MQ=255)
                aTGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCTGGATcggcag                                                                                                                                           >  6:201694/1‑141 (MQ=255)
                aTGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCTGGATcggcag                                                                                                                                           <  5:201694/143‑3 (MQ=255)
                aTGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTAt                                                                                                                                     >  1:21259/1‑149 (MQ=255)
                 tGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATc                                                                                                                                    >  2:284249/1‑149 (MQ=255)
                  ggCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATcc                                                                                                                                   >  5:15813/1‑149 (MQ=255)
                   gccgccAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGCTCCAGATcggcgg                                                                                                                                           >  3:280799/1‑140 (MQ=255)
                   gccgccAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATcc                                                                                                                                   <  5:138620/148‑1 (MQ=255)
                   gccgccAAGTGGGGCGAGGGGGTCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggcgg                                                                                                                                           <  4:280799/140‑1 (MQ=255)
                    ccgccAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGt                                                                                                                                 >  6:165358/1‑149 (MQ=255)
                    ccgccAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCAGGTCGACCTCAACGGCGCCTTCGAAGCCCAGCCGGTCAATGGTGAAGAGGTAACCGCCAACGCCAACCGCAACCCGAACCTGCCGATCCAGATAGGCGGCGCTATCCg                                                                                                                                  >  1:355834/1‑148 (MQ=255)
                        cAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGc                                                                                                                             <  1:8861/149‑1 (MQ=255)
                        cAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGACGGTCAATCGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGCTCGGCGGCGGTATCCGTTCGc                                                                                                                             >  6:129399/1‑149 (MQ=255)
                        cAAGTGGGTCGAGGGTGGCTGCCGCCGCATGCACCTGGTCGACCTCAACGGCGCATTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGc                                                                                                                             >  8:270057/1‑149 (MQ=255)
                         aaGTGGGTCGAGGGTGGCTGCCGCCGCCTGCTCCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCt                                                                                                                            <  5:419574/149‑1 (MQ=255)
                          aGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTg                                                                                                                           <  8:164320/149‑1 (MQ=255)
                               gTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGaaa                                                                                                                       <  1:284249/149‑1 (MQ=255)
                                tCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAAc                                                                                                                      >  1:187742/1‑148 (MQ=255)
                                tCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAAc                                                                                                                      >  3:208450/1‑148 (MQ=255)
                                   aGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGTGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATc                                                                                                                  <  3:197522/149‑1 (MQ=255)
                                   aGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATc                                                                                                                  <  6:252476/149‑1 (MQ=255)
                                        ggCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAaca                                                                                                             <  8:81808/149‑1 (MQ=255)
                                           tGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGCTCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAacacaa                                                                                                          >  8:277805/1‑147 (MQ=255)
                                             ccgccgccTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACg                                                                                                        <  2:21259/149‑1 (MQ=255)
                                             ccgccgccTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGATGGTTTCCAGCGAacggat                                                                                                          <  1:325721/147‑5 (MQ=255)
                                                ccgccTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTc                                                                                                      >  4:219241/1‑148 (MQ=255)
                                                   ccTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAgg                                                                                                  <  3:95916/149‑1 (MQ=255)
                                                     tGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggc                                                                                                                                             >  8:241788/1‑104 (MQ=255)
                                                     tGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggc                                                                                                                                             <  7:241788/104‑1 (MQ=255)
                                                     tGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCt                                                                                                >  5:383667/1‑149 (MQ=255)
                                                        aCCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGc                                                                                             <  7:167726/149‑1 (MQ=255)
                                                          cTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGt                                                                                           >  5:99901/1‑149 (MQ=255)
                                                           tGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGt                                                                                           >  8:150258/1‑148 (MQ=255)
                                                             gTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAg                                                                                        <  6:15813/149‑1 (MQ=255)
                                                                 aCCTCAACGGCGCCTTCGAAGGCCAGCCGGTAAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTAc                                                                                    >  6:105499/1‑149 (MQ=255)
                                                                      aaCGGCGCCTTCGAAGCCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCAAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGAATCCGTTAGCTGGAAACCATCGAACACTACGACAAGGCTGGCGTCAGCAACGtca                                                                                >  2:165388/1‑146 (MQ=255)
                                                                        cGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGatca                                                                             >  4:107498/1‑149 (MQ=255)
                                                                        cGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGatca                                                                             >  6:195/1‑149 (MQ=255)
                                                                          gcgcCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTTTCCGTTCGCTGGAAACCATCGATCACTACGTCAAGGCTGGCGTCAGCTACGTGatcatc                                                                           <  7:277805/149‑1 (MQ=255)
                                                                          gcgcCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAg                                                                                        >  1:375140/1‑136 (MQ=255)
                                                                          gcgcCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAg                                                                                        <  2:375140/136‑1 (MQ=255)
                                                                           cgcCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCg                                                                          <  2:187742/149‑1 (MQ=255)
                                                                            gcCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCgg                                                                         >  2:277845/1‑149 (MQ=255)
                                                                             ccTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGc                                                                        <  5:27047/149‑1 (MQ=255)
                                                                              cTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCa                                                                       >  1:246281/1‑149 (MQ=255)
                                                                              cTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCa                                                                       >  5:353536/1‑149 (MQ=255)
                                                                               ttCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCAc                                                                      >  7:366141/1‑149 (MQ=255)
                                                                                tCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAAc                                                                                                                      >  3:176334/1‑100 (MQ=255)
                                                                                tCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAAc                                                                                                                      <  4:176334/100‑1 (MQ=255)
                                                                                  gAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCaa                                                                   <  3:107498/149‑1 (MQ=255)
                                                                                      gCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGcc                                                               >  3:126707/1‑149 (MQ=255)
                                                                                        cAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGTTCATCGGCACCAAGGCCGt                                                             <  8:46739/149‑1 (MQ=255)
                                                                                         aGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTc                                                            <  5:232667/149‑1 (MQ=255)
                                                                                            cGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCaa                                                          <  7:165875/148‑1 (MQ=255)
                                                                                                          gTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTAc                                                                                    <  2:403513/108‑1 (MQ=255)
                                                                                                          gTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTAc                                                                                    >  1:403513/1‑108 (MQ=255)
                                                                                                          gTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTtcgt                                           >  1:283770/1‑149 (MQ=255)
                                                                                                           tGGTTACCGCCATCTCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACTTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTtcgtc                                          <  7:270057/149‑1 (MQ=255)
                                                                                                           tGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTtcgtc                                          <  7:150258/149‑1 (MQ=255)
                                                                                                             gTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTtcgtcg                                         <  5:168496/148‑1 (MQ=255)
                                                                                                              ttACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGcc                                                               >  6:155482/1‑125 (MQ=255)
                                                                                                              ttACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGcc                                                               <  5:155482/125‑1 (MQ=255)
                                                                                                              ttACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCt                                       >  1:211391/1‑149 (MQ=255)
                                                                                                              ttACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCt                                       >  6:420534/1‑149 (MQ=255)
                                                                                                              ttACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCt                                       <  4:208450/149‑1 (MQ=255)
                                                                                                                  cGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAg                                   <  3:219241/149‑1 (MQ=255)
                                                                                                                   gCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCTCCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGc                                  <  1:76888/149‑1 (MQ=255)
                                                                                                                    ccATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCAc                                                                      >  3:285295/1‑112 (MQ=255)
                                                                                                                    ccATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCAc                                                                      <  4:285295/112‑1 (MQ=255)
                                                                                                                     cATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGt                                >  8:152274/1‑149 (MQ=255)
                                                                                                                     cATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGt                                >  8:229513/1‑149 (MQ=255)
                                                                                                                      aTCGCCAAGCGCTACCCGAACGTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTGGTCGCTGAAGCGTg                               <  1:165388/149‑1 (MQ=255)
                                                                                                                      aTCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTg                               <  7:152274/149‑1 (MQ=255)
                                                                                                                         gCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCa                             >  3:262080/1‑148 (MQ=255)
                                                                                                                          ccAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCaaa                           <  8:129869/149‑1 (MQ=255)
                                                                                                                           cAAGCGCTACCCGAACCTGCCGATCCAGCTCGGCGTCTGTATCCGTTCTCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGTCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAg                          <  5:129399/149‑1 (MQ=255)
                                                                                                                            aaGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGc                         <  4:126707/149‑1 (MQ=255)
                                                                                                                                  tACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTccc                   >  2:31010/1‑149 (MQ=255)
                                                                                                                                         cgccgCCGATCTGGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAg            >  1:411673/5‑149 (MQ=255)
                                                                                                                                         cgccgCCGATCTGGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAg            <  2:411673/145‑1 (MQ=255)
                                                                                                                                         cccTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAg            <  7:105918/148‑1 (MQ=255)
                                                                                                                                          ccTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGTACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAgg           >  8:247962/1‑149 (MQ=255)
                                                                                                                                           ccgCCGATCTGGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGt          >  6:125552/3‑149 (MQ=255)
                                                                                                                                           cTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGt          >  5:412910/1‑149 (MQ=255)
                                                                                                                                             gCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCa        >  8:78942/1‑149 (MQ=255)
                                                                                                                                               ccGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCAt       <  5:439102/149‑1 (MQ=255)
                                                                                                                                                 gATCCAGATCGGCGGCGTTATCCGTTCGCTGGAAACCATCGAACACTACGTCCAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCAtcg     <  5:105499/149‑1 (MQ=255)
                                                                                                                                                 gATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCAtcg     <  5:420534/149‑1 (MQ=255)
                                                                                                                                                 gATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCAtcg     >  4:70072/1‑149 (MQ=255)
                                                                                                                                                     cAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCAtcgtcg  <  5:195/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGACGACCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCGCCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCATCGTCG  >  NC_002947/351601‑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: