Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,080,204 A→T 5.5% Q31L (CAG→CTG)  PP_5332 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,080,2040AT 5.5% 272.6 / 5.2 87Q31L (CAG→CTG) PP_5332hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (42/40);  new base T (3/2);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AAACTCATTGCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAG  >  NC_002947/6080056‑6080352
                                                                                                                                                    |                                                                                                                                                    
aaaCTCATTGCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGGATCAGCAGTGACCGTGCCGAACTCCTGATCGCCcagcag                                                                                                                                                     <  2:71407/150‑1 (MQ=255)
   cTCATTGCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGgc                                                                                                                                                  <  2:1165421/150‑1 (MQ=255)
    tCATTGCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGgcg                                                                                                                                                 >  2:406199/1‑150 (MQ=255)
     cATTGCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACCCCGGATCGCCCAGCAGGgcgc                                                                                                                                                >  1:1317945/1‑150 (MQ=255)
         gCATCTGCCCGGGTCTAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAg                                                                                                                                            <  2:1142725/150‑1 (MQ=255)
         gCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAg                                                                                                                                            <  1:130355/150‑1 (MQ=255)
            tCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTc                                                                                                                                         <  2:305602/150‑1 (MQ=255)
            tCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTc                                                                                                                                         >  2:453849/1‑150 (MQ=255)
              tGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCt                                                                                                                                       >  2:1253919/1‑150 (MQ=255)
                cccGGGTCCAAGGCCCACTGGGAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCa                                                                                                                                     <  2:519043/150‑1 (MQ=255)
                cccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCa                                                                                                                                     <  2:518941/150‑1 (MQ=255)
                     gTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTAc                                                                                                                                <  1:1451242/150‑1 (MQ=255)
                         aaGGCCCACTGGCAATCACTGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCt                                                                                                                                       <  2:267092/139‑1 (MQ=255)
                         aaGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCt                                                                                                                                       >  1:267092/1‑139 (MQ=255)
                              ccACTGGCAATCAATGAGAAGAACGCAATGGCTACCGGTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTa                                                                                                                       >  1:950659/1‑150 (MQ=255)
                                aCTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAgg                                                                                                                     >  2:1409699/1‑150 (MQ=255)
                                    gCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAgg                                                                                                                     <  1:1507488/146‑1 (MQ=255)
                                    gCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAgg                                                                                                                     >  2:1507488/1‑146 (MQ=255)
                                      aatcaatGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAAcc                                                                                                               >  1:1224684/1‑150 (MQ=255)
                                          aatGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAAccgccg                                                                                                           >  2:461314/1‑150 (MQ=255)
                                           atGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACcgccgc                                                                                                          <  1:453849/150‑1 (MQ=255)
                                            tGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCgct                                                                                                         >  1:1333619/1‑150 (MQ=255)
                                               gaagaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCgctgct                                                                                                      >  1:473373/1‑150 (MQ=255)
                                                aagaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAgg                                                                                                                         >  1:865535/1‑130 (MQ=255)
                                                aagaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAgg                                                                                                                         <  2:865535/130‑1 (MQ=255)
                                                  gaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGACGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCATCAGGGCGCACAGGTCCTCAGCTAc                                                                                                                                >  2:1162493/1‑121 (MQ=255)
                                                  gaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGACGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTAc                                                                                                                                <  1:1162493/121‑1 (MQ=255)
                                                    aCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGc                                                                                                 <  1:271749/150‑1 (MQ=255)
                                                     cGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGTGCGCCCTGCTGGgcgatcag                                                                                                                                            >  1:489994/1‑101 (MQ=255)
                                                     cGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGTGCGCCCTGCTGGgcgatcag                                                                                                                                            <  2:489994/106‑6 (MQ=255)
                                                     cGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAAATGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCt                                                                                                <  2:1100694/150‑1 (MQ=255)
                                                       cAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGa                                                                                              <  1:280592/150‑1 (MQ=255)
                                                              tACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGc                                                                                            <  1:297702/145‑1 (MQ=255)
                                                              tACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGc                                                                                            >  2:297702/1‑145 (MQ=255)
                                                               aCCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAgg                                                                                                                         <  1:279691/115‑1 (MQ=255)
                                                               aCCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAgg                                                                                                                         >  2:279691/1‑115 (MQ=255)
                                                                 cTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAAc                                                                                                                >  1:774622/1‑122 (MQ=255)
                                                                 cTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAAc                                                                                                                <  2:774622/122‑1 (MQ=255)
                                                                  gttCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCATCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGcc                                                                                   >  1:377615/2‑150 (MQ=255)
                                                                         tCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGc                                                                                                                    >  2:1199575/1‑110 (MQ=255)
                                                                         tCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGc                                                                                                                    <  1:1199575/110‑1 (MQ=255)
                                                                                 gAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGc                                                                    >  2:1320231/1‑150 (MQ=255)
                                                                                   gCAACATGGCGAACTCAACTGCTGGCGTATCCGCTTCCACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCaa                                                                  <  2:950659/150‑1 (MQ=255)
                                                                                     aaCATGGCGAACTCAACTGCTGGCGTAGCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGt                                                                <  2:1646495/150‑1 (MQ=255)
                                                                                      aCATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCCGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCa                                                                        <  2:1196220/141‑1 (MQ=255)
                                                                                      aCATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCa                                                                        >  1:1196220/1‑141 (MQ=255)
                                                                                        aTGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCt                                                                                                <  1:1296065/115‑1 (MQ=255)
                                                                                        aTGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCt                                                                                                >  2:1296065/1‑115 (MQ=255)
                                                                                         tGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGt                                                            >  1:450101/1‑150 (MQ=255)
                                                                                         tGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGt                                                            >  2:1602278/1‑150 (MQ=255)
                                                                                                 tCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCg                                                    >  1:606714/1‑150 (MQ=255)
                                                                                                 tCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCg                                                    >  2:1281582/1‑150 (MQ=255)
                                                                                                 tCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCg                                                    >  2:1289892/1‑150 (MQ=255)
                                                                                                  cAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCgg                                                   <  1:1289892/150‑1 (MQ=255)
                                                                                                   aaCTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAgg                                                                                                                     >  1:684211/1‑83 (MQ=255)
                                                                                                   aaCTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAgg                                                                                                                     <  2:684211/83‑1 (MQ=255)
                                                                                                   aaCTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAgg                                                                                     >  2:1263049/1‑115 (MQ=255)
                                                                                                   aaCTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAgg                                                                                     <  1:1263049/115‑1 (MQ=255)
                                                                                                                tCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCt                                     >  1:316387/1‑150 (MQ=255)
                                                                                                                tCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGACAATGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCt                                     >  1:690740/1‑150 (MQ=255)
                                                                                                                   gcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCa                                                                                                                                     >  1:1457790/1‑51 (MQ=255)
                                                                                                                   gcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCa                                                                                                                                     <  2:1457790/51‑1 (MQ=255)
                                                                                                                   gcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGctgctg                                                                                                     >  1:1380457/1‑83 (MQ=255)
                                                                                                                   gcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGctgctg                                                                                                     <  2:1380457/83‑1 (MQ=255)
                                                                                                                   gcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGc                                  <  2:450101/150‑1 (MQ=255)
                                                                                                                          ccGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGtt                           <  1:922628/150‑1 (MQ=255)
                                                                                                                          ccGTGCCGAACCCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGtt                           <  2:1317945/150‑1 (MQ=255)
                                                                                                                           cGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGttt                          <  2:1453432/150‑1 (MQ=255)
                                                                                                                                 gAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAAc                    <  1:1320231/150‑1 (MQ=255)
                                                                                                                                  aaCTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGCACCGCCGCTGCTGTGGCTGAGCGAc                                                                                         <  2:1380166/80‑1 (MQ=255)
                                                                                                                                  aaCTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGAc                                                                                         >  1:1380166/1‑80 (MQ=255)
                                                                                                                                     tCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGc                >  1:131689/1‑150 (MQ=255)
                                                                                                                                     tCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGc                >  2:1112906/1‑150 (MQ=255)
                                                                                                                                      ccTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCa               <  2:316387/150‑1 (MQ=255)
                                                                                                                                       cTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAg              >  1:1245474/1‑150 (MQ=255)
                                                                                                                                       atgtgcGCCCTGCTCGGCGCGCAGTTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGATGAGCGACCAGGCCCTCTTCCGACGGGCCAAGCGGAAACGCGCTGGACTACCGGCGTGCTGAGCGTGGTTTGGCAACCTGCAg              >  2:100536/6‑150 (MQ=255)
                                                                                                                                        tGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAgc             >  2:192620/1‑150 (MQ=255)
                                                                                                                                        tGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAgc             >  1:1279042/1‑150 (MQ=255)
                                                                                                                                        tGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAgc             >  2:578580/1‑150 (MQ=255)
                                                                                                                                         gtgcGCCCTGCTGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCa                                                                        <  2:1631276/87‑1 (MQ=255)
                                                                                                                                         gtgcGCCCTGCTGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCa                                                                        >  1:1631276/4‑90 (MQ=255)
                                                                                                                                          aTCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTcc                                                                            <  2:686227/85‑1 (MQ=255)
                                                                                                                                          aTCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTcc                                                                            >  1:686227/1‑85 (MQ=255)
                                                                                                                                             gCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAc        <  2:1333619/150‑1 (MQ=255)
                                                                                                                                             gCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAc        <  2:131689/150‑1 (MQ=255)
                                                                                                                                             gCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAc        <  2:377615/150‑1 (MQ=255)
                                                                                                                                              cccAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAcc       <  2:1224684/150‑1 (MQ=255)
                                                                                                                                               ccggcagGGCGCACAGGTCCTCAGCTACCAGCGACGATGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCCCGGGGGCAATTCGGTAGGCGCCGGGGAACCGGTTTGCTGGCCATGGTTTGCCAACCTCCAGCGAAAccc      >  1:1365607/4‑150 (MQ=255)
                                                                                                                                               ccAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAccc      >  2:510774/1‑150 (MQ=255)
                                                                                                                                                cagcagGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCt     >  2:83914/1‑150 (MQ=255)
                                                                                                                                                  gcagGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCa   <  1:406199/150‑1 (MQ=255)
                                                                                                                                                   cagGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAg  <  1:1230370/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAACTCATTGCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAG  >  NC_002947/6080056‑6080352

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