Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,296,145 C→A 5.7% E129* (GAG→TAG) ‡ PP_2023 ← glutathione S‑transferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,296,1450CA 5.7% 284.1 / 3.8 88E129* (GAG→TAG) ‡PP_2023glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base C (34/49);  new base A (2/3);  total (36/52)
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.95e-01

GGATAAGCCGCCTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCAGCCGCAGGGTTGGCCTGTAGCCTC  >  NC_002947/2295997‑2296290
                                                                                                                                                    |                                                                                                                                                  
ggATAAGCCGCCTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGGCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCg                                                                                                                                                   <  2:1341604/150‑1 (MQ=255)
 gATAAGCCGCCTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCgg                                                                                                                                                  <  1:1334209/150‑1 (MQ=255)
     aGCCGCCTCGGCGAACGTCTTGGTTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCAAGCGCCAACATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCctg                                                                                                                                             <  1:119600/150‑1 (MQ=255)
     aGCCGCCTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCAGGCCGAGGCctgg                                                                                                                                            <  2:1933885/150‑2 (MQ=255)
       ccgccTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCctgct                                                                                                                                           >  1:1929628/1‑150 (MQ=255)
       ccgccTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCctgct                                                                                                                                           >  2:896223/1‑150 (MQ=255)
         gccTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTgg                                                                                                                                         >  2:744226/1‑150 (MQ=255)
         gccTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTgg                                                                                                                                         >  2:563843/1‑150 (MQ=255)
          ccTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGc                                                                                                                                        <  1:483494/150‑1 (MQ=255)
          ccTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGc                                                                                                                                        <  1:1734658/150‑1 (MQ=255)
          ccTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGCACACCTAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGc                                                                                                                                        >  1:1089535/1‑150 (MQ=255)
           cTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCg                                                                                                                                       <  1:1033716/150‑1 (MQ=255)
           cTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCg                                                                                                                                       <  2:1602582/150‑1 (MQ=255)
           cTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCCGCCTGCTGGCg                                                                                                                                       >  2:279285/1‑150 (MQ=255)
            tCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCgg                                                                                                                                      >  1:1166044/1‑150 (MQ=255)
            tCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCgg                                                                                                                                      >  1:1708683/1‑150 (MQ=255)
               gCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTaa                                                                                                                                   <  2:863365/150‑1 (MQ=255)
                     cACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTcg                                                                                                                             <  1:563843/150‑1 (MQ=255)
                     cACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTcg                                                                                                                             <  1:1133772/150‑1 (MQ=255)
                        tGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCgg                                                                                                                          <  1:199288/150‑1 (MQ=255)
                          cGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGa                                                                                                                        <  2:473699/150‑1 (MQ=255)
                                 aCCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTcc                                                                                                                 <  2:1137539/150‑1 (MQ=255)
                                   cccGGCAAACTGGGGCATCGGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCgg                                                                                                               <  2:878996/150‑1 (MQ=255)
                                   cccGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCg                                                                                                                                       <  1:1658067/126‑1 (MQ=255)
                                   cccGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCg                                                                                                                                       >  2:1658067/1‑126 (MQ=255)
                                   cccGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCgg                                                                                                               >  1:1905367/1‑150 (MQ=255)
                                     cGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGt                                                                                                             >  1:1036063/1‑150 (MQ=255)
                                     cGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGt                                                                                                             >  1:245900/1‑150 (MQ=255)
                                     cGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGt                                                                                                             >  1:45260/1‑150 (MQ=255)
                                      ggCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCg                                                                                                                >  1:1012480/1‑146 (MQ=255)
                                      ggCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCg                                                                                                                <  2:1012480/146‑1 (MQ=255)
                                       gCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTaa                                                                                                           <  2:1402532/150‑1 (MQ=255)
                                           aCTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTTGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCTGCGAACCATTCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGtt                                                                                                       <  2:958360/150‑1 (MQ=255)
                                             tGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCg                                                                                                     <  2:1708683/150‑1 (MQ=255)
                                               gcgcATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTcc                                                                                                                 >  1:1251521/1‑136 (MQ=255)
                                               gcgcATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTcc                                                                                                                 <  2:1251521/136‑1 (MQ=255)
                                                 gcATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAAt                                                                                                                                  <  2:1686767/117‑1 (MQ=255)
                                                 gcATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAAt                                                                                                                                  >  1:1686767/1‑117 (MQ=255)
                                                 gcATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCAGGCCGAGGCCTGGCTGGCg                                                                                                                                       >  1:281940/1‑112 (MQ=38)
                                                 gcATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCAGGCCGAGGCCTGGCTGGCg                                                                                                                                       <  2:281940/112‑1 (MQ=38)
                                                  cATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCACGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATa                                                                                                <  2:1036063/150‑1 (MQ=255)
                                                  cATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATa                                                                                                >  2:1152999/1‑150 (MQ=255)
                                                    tCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACt                                                                                              <  2:1905367/150‑1 (MQ=255)
                                                       ggggCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCgg                                                                                                                          <  1:1068324/119‑1 (MQ=255)
                                                       ggggCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCgg                                                                                                                          >  2:1068324/1‑119 (MQ=255)
                                                        gggCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTa                                                                                          <  1:744226/150‑1 (MQ=255)
                                                          gcagcaACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACa                                                                                        >  1:1255362/1‑150 (MQ=255)
                                                           cagcaACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAg                                                                                       <  2:1929628/150‑1 (MQ=255)
                                                             gcaCCGCGGCATCGGCGATGCTGGGGGGGTTGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCGAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGt                                                                                     <  2:1871428/150‑1 (MQ=255)
                                                                       aTCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTgcgc                                                                           >  1:1730920/1‑150 (MQ=255)
                                                                       aTCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTgcgc                                                                           >  1:281443/1‑150 (MQ=255)
                                                                         cGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCt                                                                         <  2:616047/150‑1 (MQ=255)
                                                                          ggCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCtt                                                                        >  2:576520/1‑150 (MQ=255)
                                                                           gCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTg                                                                       <  2:1089535/150‑1 (MQ=255)
                                                                           gCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTg                                                                       >  2:1665485/1‑150 (MQ=255)
                                                                           gCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCAAGTgcgaattg                                                                       >  2:1665496/1‑145 (MQ=255)
                                                                              aTGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAAc                                                                    <  1:896223/150‑1 (MQ=255)
                                                                                gCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACgt                                                                  <  1:1270446/150‑1 (MQ=255)
                                                                                  tCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCg                                                                                                                >  2:1170278/1‑102 (MQ=255)
                                                                                  tCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCg                                                                                                                <  1:1170278/102‑1 (MQ=255)
                                                                                   cGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACgtggt                                                               <  2:1892692/150‑1 (MQ=255)
                                                                                    gggTGGTCGGCCAGCAGGTAAGCCCTCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCt                                                                         <  2:461213/139‑1 (MQ=255)
                                                                                    gggTGGTCGGCCAGCAGGTAAGCCCTCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCt                                                                         >  1:461213/1‑139 (MQ=255)
                                                                                    gggTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACgtggtg                                                              <  1:1665496/150‑1 (MQ=255)
                                                                                    gggTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACgtggtg                                                              <  1:1665485/150‑1 (MQ=255)
                                                                                               cagcagGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCAGGCCGAGGCCTGGCTGGCgga                                                                                                                                     >  2:833245/1‑67 (MQ=18)
                                                                                               cagcagGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCAGGCCGAGGCCTGGCTGGCgga                                                                                                                                     <  1:833245/68‑2 (MQ=18)
                                                                                                     gTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTAc                                                                                         >  2:1260386/1‑106 (MQ=255)
                                                                                                     gTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTAc                                                                                         <  1:1260386/106‑1 (MQ=255)
                                                                                                     gTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGt                                                                                     >  2:1074776/1‑110 (MQ=255)
                                                                                                     gTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGt                                                                                     <  1:1074776/110‑1 (MQ=255)
                                                                                                     gTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGgtgt                                                             <  1:836279/134‑1 (MQ=255)
                                                                                                     gTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGgtgt                                                             >  2:836279/1‑134 (MQ=255)
                                                                                                            cgccccgccAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCa                                      <  1:279285/150‑1 (MQ=255)
                                                                                                              ccccgccAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGc                                    <  2:1166044/150‑1 (MQ=255)
                                                                                                                ccgccAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCct                                  >  2:1210991/1‑150 (MQ=255)
                                                                                                                         aGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCAg                         >  2:1842754/1‑150 (MQ=255)
                                                                                                                           ccTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGt                                                                               >  1:160678/1‑94 (MQ=255)
                                                                                                                           ccTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGt                                                                               <  2:160678/94‑1 (MQ=255)
                                                                                                                             ttCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCa                          <  1:893119/145‑1 (MQ=255)
                                                                                                                             ttCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCa                          >  2:893119/1‑145 (MQ=255)
                                                                                                                                 aGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCAGCCGCAggg                 <  2:45260/150‑1 (MQ=255)
                                                                                                                                         cagccagGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCAGCCGCAGGGTTGGCCTg         <  2:245900/150‑1 (MQ=255)
                                                                                                                                           gccagGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCAGCCGCAGGGTTGGCCTGTa       <  2:281443/150‑1 (MQ=255)
                                                                                                                                                gCCTCGGCTTGCTGGCGGTAACGCTCGCTGTAATGCTCCGGGTAACATTCGTCATCCTTGTACAGGTTCACGTGCGCCTTGAACGCGTTGTAATCGCGAGCCGCCAGCCCTTCAGCCTGCTAAGCAGCCGCAAGTTTGGCTTgaacaaca  <  2:761709/150‑8 (MQ=255)
                                                                                                                                                gCCGCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCAGCCGCAGGGTTGGCCTGTAGCCTc  <  1:1486961/150‑1 (MQ=255)
                                                                                                                                                   tCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGt                                                                                           <  2:1361938/58‑1 (MQ=255)
                                                                                                                                                   tCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGt                                                                                           >  1:1361938/1‑58 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGATAAGCCGCCTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGGCCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCTGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCAGCCGCAGGGTTGGCCTGTAGCCTC  >  NC_002947/2295997‑2296290

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