Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 799,838 A→T 8.1% L213Q (CTG→CAG)  ispB ← octaprenyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947799,8380AT 8.1% 265.6 / 12.0 87L213Q (CTG→CAG) ispBoctaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (40/39);  new base T (4/3);  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 = 8.46e-01

CACCCTTCTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCC  >  NC_002947/799689‑799984
                                                                                                                                                     |                                                                                                                                                      
cACCCTTCTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCcagca                                                                                                                                                        >  1:2398355/1‑150 (MQ=255)
  cccTTCTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAgg                                                                                                                                                      >  2:249215/1‑150 (MQ=255)
   ccTTCTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGgt                                                                                                                                                     >  1:979232/1‑150 (MQ=255)
       cTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGTGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGgtcgtc                                                                                                                                                 <  1:3794/150‑1 (MQ=255)
               ccTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCTGGTCGACGACCTGCt                                                                                                                                         <  1:29669/150‑1 (MQ=255)
                cTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCCGGGGCGTCGCCGGTTTGCCCACGGCCAGATCGTCACCGACGATCTTGCCCAGGGTCTGCGAGTCCCCCTTGTAGTCCAGCAGGTCGTCGACCAGCTg                                                                                                                                        >  2:371608/1‑150 (MQ=255)
                cTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTg                                                                                                                                        <  2:979232/150‑1 (MQ=255)
                   gCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGACCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGaa                                                                                                                                     <  1:400468/150‑1 (MQ=255)
                    cGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAg                                                                                                                                    <  1:115009/150‑1 (MQ=255)
                     gTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAgg                                                                                                                                   <  1:1814076/150‑1 (MQ=255)
                     gTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCAGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAgg                                                                                                                                   >  2:2087298/1‑150 (MQ=255)
                      taacAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGc                                                                                                                                  >  1:2407444/4‑150 (MQ=255)
                         cAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCAc                                                                                                                               >  1:1235700/1‑150 (MQ=255)
                                   tGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCCGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATgg                                                                                                                     >  1:1823240/1‑150 (MQ=255)
                                     cTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTc                                                                                                                   <  1:677328/150‑1 (MQ=255)
                                     cTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTc                                                                                                                   <  1:1866828/150‑1 (MQ=255)
                                     cTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTc                                                                                                                   <  2:2268544/150‑1 (MQ=255)
                                       cGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGc                                                                                                                 <  1:249215/150‑1 (MQ=255)
                                       cGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGc                                                                                                                 <  1:1809039/150‑1 (MQ=255)
                                            gTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAg                                                                                                            <  1:652753/150‑1 (MQ=255)
                                            gTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAg                                                                                                            >  1:918771/1‑150 (MQ=255)
                                              gCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGt                                                                                                          <  1:850395/150‑1 (MQ=255)
                                                cTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCTGgtcgac                                                                                                                                                 <  1:1582242/109‑3 (MQ=255)
                                                cTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCTGgtcgac                                                                                                                                                 >  2:1582242/1‑107 (MQ=255)
                                                cTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTAc                                                                                                        >  1:1385788/1‑150 (MQ=255)
                                                    gcgcATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAgg                                                                                                                                   >  1:1911957/1‑119 (MQ=255)
                                                    gcgcATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAgg                                                                                                                                   <  2:1911957/119‑1 (MQ=255)
                                                          ggTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCCCCCGTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTc                                                                                              <  2:1994794/150‑1 (MQ=255)
                                                             gtAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAggg                                                                                                  <  2:587331/143‑1 (MQ=255)
                                                             gtAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAggg                                                                                                  >  1:587331/1‑143 (MQ=255)
                                                              tAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑gcgc                                                                                        <  2:843419/150‑1 (MQ=255)
                                                                gATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATg                                                                                                                      <  1:1207417/120‑1 (MQ=255)
                                                                gATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATg                                                                                                                      >  2:1207417/1‑120 (MQ=255)
                                                                  tCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCgctgct                                                                                    >  1:1630442/1‑150 (MQ=255)
                                                                        gCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCgg                                                                             >  2:1885247/1‑150 (MQ=255)
                                                                          aGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTg                                                                           <  2:1823240/150‑1 (MQ=255)
                                                                            ggTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAgg                                                                                                                                   >  1:888545/1‑95 (MQ=255)
                                                                            ggTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAgg                                                                                                                                   <  2:888545/95‑1 (MQ=255)
                                                                                  tttGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGc                                                                   >  1:1240017/1‑150 (MQ=255)
                                                                                       ccTCGGCCAGGTCGTCAGCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAgcg                                                              >  2:1238434/1‑150 (MQ=255)
                                                                                        cTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGTTACGCAGTTCGTG‑‑TCGCTGCTCGTCGGTGGACCGGGCCCgcgc                                                             >  2:1600704/1‑150 (MQ=255)
                                                                                        cTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAgcgc                                                             >  2:1371620/1‑150 (MQ=255)
                                                                                        cTCGGCCAGGTCGTCACCGACGTTCTTACCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAgcgc                                                             >  1:2275442/1‑150 (MQ=255)
                                                                                         tCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAg                                                                                                            <  2:1887328/105‑1 (MQ=255)
                                                                                         tCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAg                                                                                                            >  1:1887328/1‑105 (MQ=255)
                                                                                             ccAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCgcag                                                        >  1:2196610/1‑150 (MQ=255)
                                                                                                   cgtcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgt                                                 <  1:1885247/150‑1 (MQ=255)
                                                                                                    gtcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgtg                                                >  1:794225/1‑150 (MQ=255)
                                                                                                    gtcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgtg                                                >  2:1112894/1‑150 (MQ=255)
                                                                                                     tcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTgg                                               <  2:794225/150‑1 (MQ=255)
                                                                                                       aCCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCAATGTGGGt                                             >  2:904781/1‑150 (MQ=255)
                                                                                                         cGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgt                                                 >  1:1032964/1‑144 (MQ=255)
                                                                                                         cGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgt                                                 <  2:1032964/144‑1 (MQ=255)
                                                                                                         cGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAg                                           >  1:1146734/1‑150 (MQ=255)
                                                                                                         cGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAg                                           >  1:2306309/1‑150 (MQ=255)
                                                                                                          gACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGa                                          >  2:793527/1‑150 (MQ=255)
                                                                                                                cTTGCTCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGGCCCTAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                    <  1:162741/150‑1 (MQ=255)
                                                                                                                cTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                    <  2:1385788/150‑1 (MQ=255)
                                                                                                                   gCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGaa                                 <  2:1235700/150‑1 (MQ=255)
                                                                                                                     ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCtt                                     <  1:1968137/144‑1 (MQ=255)
                                                                                                                     ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCtt                                     >  2:1968137/1‑144 (MQ=255)
                                                                                                                     ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACa                               >  1:829175/1‑150 (MQ=255)
                                                                                                                     ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACa                               <  1:518994/150‑1 (MQ=255)
                                                                                                                     ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGCTGGccca                                                                      >  2:1591028/1‑110 (MQ=255)
                                                                                                                     ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGCTGGccca                                                                      <  1:1591028/112‑3 (MQ=255)
                                                                                                                       aGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGc                             <  2:2398355/150‑1 (MQ=255)
                                                                                                                          gTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATc                          <  1:1603236/150‑1 (MQ=255)
                                                                                                                           tCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCAAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCg                         >  1:1230680/1‑150 (MQ=255)
                                                                                                                            cTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGc                        <  2:2407444/150‑1 (MQ=255)
                                                                                                                             tGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgtg                                                >  2:515461/1‑125 (MQ=255)
                                                                                                                             tGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgtg                                                <  1:515461/125‑1 (MQ=255)
                                                                                                                                gAGTCACCCTTGTAGTCCAGCAGGTCGTCGAACAGCTGGAAGGCCACGCCCAGATGGTCGGAAACGGTACGCAGGACATCATGGGCTACTGCATCACTGCGCTG‑‑‑CAGCGCCGCAGCACTGTGGGTAGAAGGTTCGGACGGCATCGCCGTc                    >  1:1151734/1‑150 (MQ=255)
                                                                                                                                  gTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTa                                            >  1:164888/1‑124 (MQ=255)
                                                                                                                                  gTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTa                                            <  2:164888/124‑1 (MQ=255)
                                                                                                                                     aCCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGa                                                                                                                        <  2:975489/49‑1 (MQ=255)
                                                                                                                                     aCCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGa                                                                                                                        >  1:975489/1‑49 (MQ=255)
                                                                                                                                     aCCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGcc               >  2:2362232/1‑150 (MQ=255)
                                                                                                                                       ccTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCCGAGCCGCAGGCTGACGTGTAGAAGCTTCTAACAGCAGCGCCGTCTTGCcgc             >  2:1833267/1‑150 (MQ=255)
                                                                                                                                       ccTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTcgccgc             >  1:2196270/1‑150 (MQ=255)
                                                                                                                                         ttGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGa           <  1:793527/150‑1 (MQ=255)
                                                                                                                                           gTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCtt                  >  2:1156272/1‑141 (MQ=255)
                                                                                                                                           gTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCtt                  <  1:1156272/141‑1 (MQ=255)
                                                                                                                                           gTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAAt         <  2:2196610/150‑1 (MQ=255)
                                                                                                                                             cttCCAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGa       >  1:1907513/3‑150 (MQ=255)
                                                                                                                                             cttCCAGCTGGTCGACGACCTGCTGGAAGGACACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGa       >  1:1635436/3‑150 (MQ=255)
                                                                                                                                             aGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGa       <  2:1630442/150‑1 (MQ=255)
                                                                                                                                               tCCAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTgg                                               >  1:215827/1‑108 (MQ=255)
                                                                                                                                               tCCAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTgg                                               <  2:215827/108‑1 (MQ=255)
                                                                                                                                                 cagctgGTCGACGACCTGCTGGACGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTgg                                                                     <  1:336269/85‑1 (MQ=38)
                                                                                                                                                 cagctgGTCGACGACCTGCTGGACGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTgg                                                                     >  2:336269/1‑85 (MQ=38)
                                                                                                                                                 cagctgGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATg        >  1:1674396/1‑145 (MQ=255)
                                                                                                                                                 cagctgGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATg        <  2:1674396/145‑1 (MQ=255)
                                                                                                                                                 cagctgGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTc   >  1:1668847/1‑150 (MQ=255)
                                                                                                                                                 cagcagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGTATCGCCGTCTTGCCGCGAATGACCTc   <  2:1146734/150‑1 (MQ=255)
                                                                                                                                                 cagcagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTc   >  2:1752577/1‑150 (MQ=255)
                                                                                                                                                  agctgGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGAcc     >  2:225295/1‑147 (MQ=255)
                                                                                                                                                  agctgGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTcc  >  2:2373106/1‑150 (MQ=255)
                                                                                                                                                  agcagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                    <  2:2348989/116‑1 (MQ=255)
                                                                                                                                                  agcagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCATGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                    >  1:2348989/1‑116 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CACCCTTCTGGATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTC‑‑GCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCC  >  NC_002947/799689‑799984

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: