Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 351,753 A→G 6.2% 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,7530AG 6.2% 238.5 / 4.3 81Q79R (CAG→CGG) ‡hisAphosphoribosylformimino‑5‑aminoimidazole carboxamide ribonucleotide isomerase
Reads supporting (aligned to +/- strand):  ref base A (38/38);  new base G (2/3);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

GACGACCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCATCGTCGGC  >  NC_002947/351605‑351901
                                                                                                                                                      |                                                                                                                                                    
gacgacCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCCATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCTAACCTGCCGATCCAg                                                                                                                                                     >  1:2230234/1‑150 (MQ=255)
gacgacCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAg                                                                                                                                                     >  2:1234212/1‑150 (MQ=255)
 acgacCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGa                                                                                                                                                    <  1:805423/150‑1 (MQ=255)
   gacCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CGTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCACCACCAAGCCCTACCCGCACCTGCCGATCCAGATc                                                                                                                                                  >  2:1026687/1‑150 (MQ=255)
    acCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGCAATCAACGGAGAACACACGTCTGAACTCCAGTCACTCTAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcg                                                                                                                                                 >  1:531600/1‑150 (MQ=255)
       cGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggcg                                                                                                                                              <  1:1259299/150‑1 (MQ=255)
       cGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggcg                                                                                                                                              >  2:1819149/1‑150 (MQ=255)
       cGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggcg                                                                                                                                              >  2:1146606/1‑150 (MQ=255)
         gTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAGGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggcggc                                                                                                                                            >  2:1330316/1‑150 (MQ=255)
          tGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCTAACCTGCCGATCCAGATcggcggcg                                                                                                                                           <  2:2230234/150‑1 (MQ=255)
          tGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACATCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATcggcggcg                                                                                                                                           <  1:1538570/150‑1 (MQ=255)
             gCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACTCGAACCCGCCGATCTGGATcggcaggttcg                                                                                                                                        <  1:298408/150‑8 (MQ=255)
                tGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATcc                                                                                                                                     <  2:1732057/150‑1 (MQ=255)
                       cAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCt                                                                                                                              >  1:1900657/1‑150 (MQ=255)
                              gTCGAGGTTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCTGGATcggcag                                                                                                                                             <  2:457225/128‑3 (MQ=255)
                              gTCGAGGTTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCTGGATcggcag                                                                                                                                             >  1:457225/1‑126 (MQ=255)
                               tCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCT‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTTGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCa                                                                                                                      <  1:1234212/150‑1 (MQ=255)
                               tCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCa                                                                                                                      <  1:2274575/150‑1 (MQ=255)
                                    ggTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGaa                                                                                                                 <  1:439976/150‑1 (MQ=255)
                                          tGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTAc                                                                                                           <  1:1752948/150‑1 (MQ=255)
                                             cgccgccTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTc                                                                                                        >  1:1308334/1‑150 (MQ=255)
                                             cgccgccTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTc                                                                                                        >  1:2107491/1‑150 (MQ=255)
                                                cgccTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAg                                                                                                     <  2:234665/150‑1 (MQ=255)
                                                 gccTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAgg                                                                                                    <  1:1146606/150‑1 (MQ=255)
                                                     gCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTgg                                                                                                >  1:1155829/1‑150 (MQ=255)
                                                      cACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGa                                                                                                                                                    >  1:1091817/1‑97 (MQ=255)
                                                      cACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGa                                                                                                                                                    <  2:1091817/97‑1 (MQ=255)
                                                      cACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGc                                                                                               <  2:1308334/150‑1 (MQ=255)
                                                       aCCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTg                                                                                                                             >  1:1649137/1‑119 (MQ=255)
                                                       aCCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTg                                                                                                                             <  2:1649137/119‑1 (MQ=255)
                                                        ccTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGGAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGt                                                                                             <  1:954919/150‑1 (MQ=255)
                                                         cTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTc                                                                                            >  1:319961/1‑150 (MQ=255)
                                                           ggTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAg                                                                                          <  1:459631/150‑1 (MQ=255)
                                                            gTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGc                                                                                         <  1:1965524/150‑1 (MQ=255)
                                                              cGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTa                                                                                       <  2:753097/150‑1 (MQ=255)
                                                                aCCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACg                                                                                     <  2:156611/150‑1 (MQ=255)
                                                                aCCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCCGCCGATCTGGATcggcag                                                                                                                                             <  1:1303436/94‑3 (MQ=255)
                                                                aCCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCCGCCGATCTGGATcggcag                                                                                                                                             >  2:1303436/1‑92 (MQ=255)
                                                                 ccTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGt                                                                                    >  1:1414456/1‑150 (MQ=255)
                                                                  cTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCa                                                                                                                                                      <  1:592742/83‑1 (MQ=255)
                                                                  cTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCa                                                                                                                                                      >  2:592742/1‑83 (MQ=255)
                                                                  cTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTg                                                                                   >  1:1690302/1‑150 (MQ=255)
                                                                    cAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGCAACTGCAAGGCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGat                                                                                 <  1:869845/150‑1 (MQ=255)
                                                                        ggCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGattatc                                                                             >  2:2112914/1‑150 (MQ=255)
                                                                                   gAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATc                                                                                                                    >  2:1464160/1‑102 (MQ=255)
                                                                                   gAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATc                                                                                                                    <  1:1464160/102‑1 (MQ=255)
                                                                                      ggCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGcc                                                                 <  2:30814/150‑1 (MQ=255)
                                                                                          aGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCAAAACCTGCCGATCCAGATGGGCGGGAGTGAGCAATCGCTGGAAACCATCGAAGACTACGTCAAGGCTGGCGCCAGCGACGTGATAATTGGCACCAAGGCCGTCa                                                             >  1:1873226/1‑150 (MQ=255)
                                                                                                tCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAagcagc                                                       <  2:319961/150‑1 (MQ=255)
                                                                                                  aaTGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATcc                                                                                                                                     <  2:1746238/70‑1 (MQ=255)
                                                                                                  aaTGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATcc                                                                                                                                     >  1:1746238/1‑70 (MQ=255)
                                                                                                  aaTGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGaaa                                                                                                                         <  2:2108919/82‑1 (MQ=255)
                                                                                                  aaTGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGaaa                                                                                                                         >  1:2108919/1‑82 (MQ=255)
                                                                                                   aTGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCTGTTCGATGGAAACCATCGAACACTACGTCAAGGCTGCCGTCAGCTACGTGATCATCGGCACCAAGGCGGTCAAGCAGCCCg                                                    >  2:288666/1‑150 (MQ=255)
                                                                                                     ggTGAAGTGGTTACCGCCATAGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAg                                                  >  1:1975670/1‑150 (MQ=255)
                                                                                                      gTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAagc                                                          <  1:1926666/141‑1 (MQ=255)
                                                                                                      gTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAagc                                                          >  2:1926666/1‑141 (MQ=255)
                                                                                                          ggTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCCTACACCACTACGTCAAGGCTGGCGTCAGCCACGTGACCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTtcgt                                             <  1:1330316/149‑1 (MQ=255)
                                                                                                              gTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCt                                         >  2:1498087/1‑150 (MQ=255)
                                                                                                              gTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCt                                         >  1:1478945/1‑150 (MQ=255)
                                                                                                                  ccGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGGTACGTGATCCTCCCCCCAAAGGCCGTGAAGCAGCCCGAGTTCGTCGCTGCAg                                     >  2:774900/1‑150 (MQ=255)
                                                                                                                   cGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGcc                                                                 >  1:231724/1‑121 (MQ=255)
                                                                                                                   cGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGcc                                                                 <  2:231724/121‑1 (MQ=255)
                                                                                                                   cGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGc                                    <  2:531600/150‑1 (MQ=255)
                                                                                                                    gCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTg                                                                                                 >  2:968236/1‑88 (MQ=255)
                                                                                                                    gCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTg                                                                                                 <  1:968236/88‑1 (MQ=255)
                                                                                                                     ccATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCa                                                                         <  2:1534107/111‑1 (MQ=255)
                                                                                                                     ccATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCa                                                                         >  1:1534107/1‑111 (MQ=255)
                                                                                                                     ccATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGt                                  <  2:1867620/150‑1 (MQ=255)
                                                                                                                      cATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTg                                 >  1:497323/1‑150 (MQ=255)
                                                                                                                            cAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGc                           >  2:1975038/1‑150 (MQ=255)
                                                                                                                            cAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGc                           <  2:617748/150‑1 (MQ=255)
                                                                                                                              aGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGTCGTCAGCTACGTGATCCTCGGAAATTCAGCGCCCAATCAGCCCTAGTTCGTCGCTGAAGCGTGCAAAGCCt                         >  1:2270952/1‑150 (MQ=255)
                                                                                                                                  cTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTccc                     >  2:120911/1‑150 (MQ=255)
                                                                                                                                         aaCCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAg              <  1:1819149/150‑1 (MQ=255)
                                                                                                                                              gCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGatc                                                                                <  1:1994379/79‑1 (MQ=255)
                                                                                                                                              gCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGatc                                                                                >  2:1994379/1‑79 (MQ=255)
                                                                                                                                                cGATCCAGATCGGCGGCGGCATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCAtcg       >  1:72337/1‑150 (MQ=255)
                                                                                                                                                  aTCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTtc                                               >  1:57218/1‑108 (MQ=255)
                                                                                                                                                  aTCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTtc                                               <  2:57218/108‑1 (MQ=255)
                                                                                                                                                  aTCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCAtcgtc     <  2:575974/150‑1 (MQ=255)
                                                                                                                                                  aTCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCAtcgtc     >  1:2214256/1‑150 (MQ=255)
                                                                                                                                                     cAGATCGGCGGCGGTATCCGTTCGCTGCAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCATCGTCGGc  <  2:1783385/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GACGACCCGGTGAGCATGGCCGCCAAGTGGGTCGAGGGTGGCTGCCGCCGCCTGCACCTGGTCGACCTCAACGGCG‑‑CCTTCGAAGGCCAGCCGGTCAATGGTGAAGTGGTTACCGCCATCGCCAAGCGCTACCCGAACCTGCCGATCCAGATCGGCGGCGGTATCCGTTCGCTGGAAACCATCGAACACTACGTCAAGGCTGGCGTCAGCTACGTGATCATCGGCACCAAGGCCGTCAAGCAGCCCGAGTTCGTCGCTGAAGCGTGCAAAGCCTTCCCGGGCAAGGTCATCGTCGGC  >  NC_002947/351605‑351901

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: