Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 795,308 G→C 7.0% A483P (GCG→CCG)  PP_0680 → ATP‑dependent protease‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947795,3080GC 7.0% 220.6 / 6.3 71A483P (GCG→CCG) PP_0680ATP‑dependent protease‑like protein
Reads supporting (aligned to +/- strand):  ref base G (33/33);  new base C (2/3);  total (35/36)
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.66e-01

ACGACCCGGACTTCCAGGAGATGTTCCGCGTGCTGGTGGACTTCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCGCCAGCTGGCCAGTGACGAAATGACCGACGCAG  >  NC_002947/795159‑795443
                                                                                                                                                     |                                                                                                                                       
aCGACCCGGACTTCCAGGAGATGTTCCGCGTGCTGGTGGACTTCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAg                                                                                                                                         >  1:716832/1‑150 (MQ=255)
     ccGGACTTCCAGGAGATGTTCCGCGTGCTGGTGGACTTCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGt                                                                                                                                    >  2:1149093/1‑150 (MQ=255)
       ggACTGCCAGGGGATGTGCCGCGTGCTGGTGGACTTCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCt                                                                                                                                  <  2:441021/150‑1 (MQ=255)
           ttCCAGGAGATGTTCCGCGTGCTGGTGGACTTCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGcc                                                                                                                              >  2:288220/1‑150 (MQ=255)
               aGGAGATGTTCCGCGTGCTGGTGGACTTCGATGAAGAAATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCt                                                                                                                          >  2:1144938/1‑150 (MQ=255)
                 gagaTGTTCCGCGTGCTGGTGGACTTCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTAc                                                                                                                        >  2:1282805/1‑150 (MQ=255)
                 gagaTGTTCCGCGTGCTGGTGGACTTCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCACCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTAc                                                                                                                        >  2:1164159/1‑150 (MQ=255)
                  agaTGTTCCGCGTGCTGGTGGACTTCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCACCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACa                                                                                                                       <  1:1164159/150‑1 (MQ=255)
                        tCCGCGTGCTGGTGGACTTCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTc                                                                                                                 <  1:119025/150‑1 (MQ=255)
                                 tggtggACTTCGATGAAGACATGCCGATGGTCGACGAGAACCTGGCGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAGCGAAGAGGGCATGGGGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACATTGCTCGTCTGGCCg                                                                                                        >  1:579729/1‑150 (MQ=255)
                                         ttCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCAc                                                                                                                            >  2:282465/1‑122 (MQ=255)
                                         ttCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCAc                                                                                                                            <  1:282465/122‑1 (MQ=255)
                                         ttCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTc                                                                                                                 >  2:2678/1‑133 (MQ=255)
                                         ttCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTc                                                                                                                 <  1:2678/133‑1 (MQ=255)
                                           cGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTGCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTc                                                                                              <  2:1226011/150‑1 (MQ=255)
                                            gATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCAc                                                                                                                            <  2:971141/119‑1 (MQ=255)
                                            gATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCAc                                                                                                                            >  1:971141/1‑119 (MQ=255)
                                                 aGACATGCCGATGTTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCCCCGCTGTAAGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCt                                                                                        <  1:88168/150‑1 (MQ=255)
                                                   aCATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCTCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGt                                                                                      <  2:867238/150‑1 (MQ=255)
                                                     aTGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGACTGTCa                                                                                    >  1:1077919/1‑150 (MQ=255)
                                                        ccGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAgcg                                                                                 >  1:504727/1‑150 (MQ=255)
                                                         cGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAgcgc                                                                                <  1:1231644/150‑1 (MQ=255)
                                                         cGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAgcgc                                                                                <  1:1149093/150‑1 (MQ=255)
                                                         cGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAgcgc                                                                                <  1:94474/150‑1 (MQ=255)
                                                           aTGATCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCAGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCGCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGcgcgg                                                                              >  2:1355987/1‑149 (MQ=255)
                                                               gcgacgaGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGGCGCGCCTGTCAGCGCGTATCg                                                                          >  2:673190/2‑150 (MQ=255)
                                                                 gactaGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGCTACCGCGTcact                                                                                                                                <  1:1386628/94‑4 (MQ=255)
                                                                 gacgaGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAAc                                                                                                   >  1:1004102/1‑123 (MQ=255)
                                                                 gacgaGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAAc                                                                                                   <  2:1004102/123‑1 (MQ=255)
                                                                 gacgaGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGCTACCGCGTcact                                                                                                                                >  2:1386628/1‑91 (MQ=255)
                                                                       aaCCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGAGACGCGCTACCGCGTcactg                                                                                                                               >  2:1299850/1‑85 (MQ=255)
                                                                       aaCCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGAGACGCGCTACCGCGTcactg                                                                                                                               <  1:1299850/89‑5 (MQ=255)
                                                                          cTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTc                                                               >  1:1006138/1‑150 (MQ=255)
                                                                          cTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTc                                                               >  1:751404/1‑150 (MQ=255)
                                                                           tGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTcc                                                              <  1:39563/150‑1 (MQ=255)
                                                                                cAGTTCGCCAAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTGTGGCCAACTACAGTGCTCGTCTGGCCGAAAACCAGTAGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTg                                                         <  2:623865/150‑1 (MQ=255)
                                                                                        ccAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGcc                                                                                         <  2:1353743/110‑1 (MQ=255)
                                                                                        ccAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGcc                                                                                         >  1:1353743/1‑110 (MQ=255)
                                                                                        ccAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGa                                                 >  2:369373/1‑150 (MQ=255)
                                                                                         cAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAg                                                <  2:716832/150‑1 (MQ=255)
                                                                                              gTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGa                                           <  1:288220/150‑1 (MQ=255)
                                                                                               ttGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGAt                                          >  1:212916/1‑150 (MQ=255)
                                                                                                  cgcacccgcaccAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTc                                       >  2:2696/1‑150 (MQ=255)
                                                                                                  cgcacccgcaccAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTc                                       <  1:2696/150‑1 (MQ=255)
                                                                                                   gcacccgcaccAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTcgc                                                                                            >  2:72971/1‑96 (MQ=255)
                                                                                                   gcacccgcaccAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTcgc                                                                                            <  1:72971/96‑1 (MQ=255)
                                                                                                      ccggcaccAGCGAAGAGGCCATGGCGCCGCTGACCAATGACTCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTTTTCCAGCTGTTCAGCGAGGCTGATTTCATTc                                   <  2:579729/147‑1 (MQ=255)
                                                                                                      cccgcaccAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTc                                   <  2:1069633/150‑1 (MQ=255)
                                                                                                      cccgcaccAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTc                                   >  1:510289/1‑150 (MQ=255)
                                                                                                       ccgcaccAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTa                                                                             >  1:1010954/1‑107 (MQ=255)
                                                                                                       ccgcaccAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTa                                                                             <  2:1010954/107‑1 (MQ=255)
                                                                                                       ccgcaccAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCt                                                          >  2:1219563/1‑126 (MQ=255)
                                                                                                       ccgcaccAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCt                                                          <  1:1219563/126‑1 (MQ=255)
                                                                                                       ccgcaccAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCg                                  >  2:1111413/1‑150 (MQ=255)
                                                                                                        cgcaccAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCGc                                 <  1:1282805/150‑1 (MQ=255)
                                                                                                             taaCGAAGAGGGCATGGCGCCGCTGACCAGTGACGCCCTCCGGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCGCCAGCt                            <  1:374565/149‑1 (MQ=255)
                                                                                                               aCGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCCGGGCATATCCGTAAAAGGCCACAGCTGGGCAGCAAGACTGATTTATTTAGCctgctgc                          >  2:374565/1‑149 (MQ=255)
                                                                                                                    gagGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCGCCAGCTGGCCAGt                     >  2:1019022/1‑150 (MQ=255)
                                                                                                                       ggCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTc                                   <  1:537697/133‑1 (MQ=255)
                                                                                                                       ggCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTc                                   >  2:537697/1‑133 (MQ=255)
                                                                                                                        gCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCGCCAGCTGGCCAGTGACg                 >  1:1426580/1‑150 (MQ=255)
                                                                                                                           tGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCGCCAGCTGGCCAGTGACGaaa              <  2:510289/150‑1 (MQ=255)
                                                                                                                               gccgcTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCGCCAGCTGGCCAGTGACGAAATGAc          <  2:504727/150‑1 (MQ=255)
                                                                                                                               gccgcTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCGCCAGCTGGCCAGTGACGAAATGAc          <  1:369373/150‑1 (MQ=255)
                                                                                                                                     gACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCGCCAGCTGGCCAGTGACGAAATGACCGACGc    >  1:561509/1‑150 (MQ=255)
                                                                                                                                       ccAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCGCCAGCTGGCCAGTGACGAAATGACCGACGCAg  <  2:1006138/150‑1 (MQ=255)
                                                                                                                                        cAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTc                                   >  2:746911/1‑116 (MQ=255)
                                                                                                                                        cAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTc                                   <  1:746911/116‑1 (MQ=255)
                                                                                                                                               gcgGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTg                                                         >  1:1402389/1‑87 (MQ=255)
                                                                                                                                               gcgGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTg                                                         <  2:1402389/87‑1 (MQ=255)
                                                                                                                                                 ggTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCGCCAGCTGGCCAGTGACGAAATGACCGAc      <  2:76261/136‑1 (MQ=255)
                                                                                                                                                 ggTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCGCCAGCTGGCCAGTGACGAAATGACCGAc      >  1:76261/1‑136 (MQ=255)
                                                                                                                                                     |                                                                                                                                       
ACGACCCGGACTTCCAGGAGATGTTCCGCGTGCTGGTGGACTTCGATGAAGACATGCCGATGGTCGACGAGAACCTGGAGCAGTTCGCCCAGCTGTTGCGCACCCGCACCAACGAAGAGGGCATGGCGCCGCTGACCAGTGACGCGGTAGCGCGTCTGGCCACCTACAGTGCTCGTCTGGCCGAAAACCAGTCGCGCCTGTCAGCGCGTATCGGTGACCTGTTCCAGCTGGTCAGCGAGGCTGATTTCATTCGCCAGCTGGCCAGTGACGAAATGACCGACGCAG  >  NC_002947/795159‑795443

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: