Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,688,376 A→C 6.6% F216V (TTC→GTC)  pilJ ← twitching motility protein PilJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,688,3760AC 6.6% 219.5 / 5.7 77F216V (TTC→GTC) pilJtwitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base A (33/38);  new base C (2/3);  total (35/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.40e-01

CACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGA  >  NC_002947/5688227‑5688516
                                                                                                                                                      |                                                                                                                                            
cACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGaa                                                                                                                                              >  2:1418576/1‑150 (MQ=255)
 aCCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTTGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTAACGGCCGaaa                                                                                                                                             >  1:213474/1‑149 (MQ=255)
 aCCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAg                                                                                                                                             <  1:557148/150‑1 (MQ=255)
 aCCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAg                                                                                                                                             <  1:556940/150‑1 (MQ=255)
 aCCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAg                                                                                                                                             <  1:30289/150‑1 (MQ=255)
  ccGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAgg                                                                                                                                            <  1:682981/150‑1 (MQ=255)
  ccGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAgg                                                                                                                                            >  2:167957/1‑150 (MQ=255)
       ccGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCg                                                                                                                                       <  1:1075171/150‑1 (MQ=255)
           cAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGACGCCTTCGgccg                                                                                                                                   <  1:575968/150‑1 (MQ=255)
             aCGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGgccgcc                                                                                                                                 >  2:1368482/1‑150 (MQ=255)
                         aacTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCGTAGACGCGGGTGACCTGGATTGTCGCATTGCCCTTGATCGTGCTTTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTg                                                                                                                     >  1:324328/1‑150 (MQ=255)
                           cTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTgtc                                                                                                                   <  2:1088013/150‑1 (MQ=255)
                            tCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGggcg                                                                                                                           >  1:900528/1‑141 (MQ=255)
                            tCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGggcg                                                                                                                           <  2:900528/141‑1 (MQ=255)
                            tCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTgtcg                                                                                                                  >  2:5364/1‑150 (MQ=255)
                               gCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGGCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCa                                                                                                               >  2:125788/1‑150 (MQ=255)
                                               gcgcACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGgtgt                                                                                               <  2:542465/150‑1 (MQ=255)
                                               gcgcACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGgtgt                                                                                               >  1:638912/1‑150 (MQ=255)
                                                      cgcgTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCGGTCGTCACCGGCCAGCACGGTGTTGACCGa                                                                                        <  2:92516/150‑1 (MQ=255)
                                                          tCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCa                                                                                    >  2:620416/1‑150 (MQ=255)
                                                           cGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCaa                                                                                   <  2:358007/150‑1 (MQ=255)
                                                               gtcgaggACGCCGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGAt                                                                               <  1:1021169/144‑1 (MQ=255)
                                                               gTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGAt                                                                               <  2:849528/150‑1 (MQ=255)
                                                               gTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGAt                                                                               <  1:620416/150‑1 (MQ=255)
                                                                   tCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGt                                                                           >  1:494820/1‑150 (MQ=255)
                                                                   tCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGt                                                                           >  1:1103936/1‑150 (MQ=255)
                                                                     gACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTc                                                                         >  1:76975/1‑150 (MQ=255)
                                                                               aCCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAAc                                                               <  1:5364/150‑1 (MQ=255)
                                                                                 cTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTgtc                                                                                                                   <  1:1196030/96‑1 (MQ=255)
                                                                                 cTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTgtc                                                                                                                   >  2:1196030/1‑96 (MQ=255)
                                                                                 cTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGc                                                                      <  1:1301695/141‑1 (MQ=255)
                                                                                 cTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAg                                                             >  1:993655/1‑150 (MQ=255)
                                                                                 cTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGGGTTTACCGAGCCAAGGATGCGTTCGGc                                                                      >  2:1301695/1‑141 (MQ=255)
                                                                                   ggATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGaa                                                                                                                                              >  2:981951/1‑67 (MQ=255)
                                                                                   ggATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGaa                                                                                                                                              <  1:981951/67‑1 (MQ=255)
                                                                                   ggATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCt                                                           <  1:1388620/150‑1 (MQ=255)
                                                                                    gATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTg                                                          <  1:995741/150‑1 (MQ=255)
                                                                                    gATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTg                                                          <  1:1368482/150‑1 (MQ=255)
                                                                                         tCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGACGCCTTCGgccgt                                                                                                                                  >  2:317093/1‑72 (MQ=38)
                                                                                         tCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGACGCCTTCGgccgt                                                                                                                                  <  1:317093/73‑2 (MQ=38)
                                                                                           gCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTgg                                                   >  2:434312/1‑150 (MQ=255)
                                                                                           gCATTGCCATTGATCATGCCCGCCAGCACCTGGCCAAAGAGACCGGCCTGCAATGCTCTTGGTTTGCCGCCCGGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTgg                                                   <  1:167957/150‑1 (MQ=255)
                                                                                            cATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTggg                                                  <  2:638912/150‑1 (MQ=255)
                                                                                             aTTGCCATTGATCCTGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAg                                                             <  1:930945/138‑1 (MQ=255)
                                                                                             aTTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAg                                                             >  2:930945/1‑138 (MQ=255)
                                                                                                  cATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCg                                            >  1:1215929/1‑150 (MQ=255)
                                                                                                    ttGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCa                                                                    <  2:1369364/124‑1 (MQ=255)
                                                                                                    ttGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCa                                                                    >  1:1369364/1‑124 (MQ=255)
                                                                                                    ttGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTAACCTGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCa                                          >  1:246539/1‑150 (MQ=255)
                                                                                                        tCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCt                                      <  1:1258013/150‑1 (MQ=255)
                                                                                                          aTGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGa                                    <  1:1297187/150‑1 (MQ=255)
                                                                                                          aTGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGa                                    >  1:504657/1‑150 (MQ=255)
                                                                                                           tGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGAt                                   <  1:1316682/150‑1 (MQ=255)
                                                                                                               cTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGa                                                                                <  2:1084068/101‑1 (MQ=255)
                                                                                                               cTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGa                                                                                >  1:1084068/1‑101 (MQ=255)
                                                                                                                tCCAGCACCTGGCCATAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGcc                              <  2:993655/150‑1 (MQ=255)
                                                                                                                tCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACTGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGcc                              >  1:879110/1‑150 (MQ=255)
                                                                                                                tCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGcc                              >  2:614019/1‑150 (MQ=255)
                                                                                                                 ccAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCg                             >  2:1260796/1‑150 (MQ=255)
                                                                                                                 ccAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCg                             <  2:879110/150‑1 (MQ=255)
                                                                                                                  cAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCgg                            <  1:434312/150‑1 (MQ=255)
                                                                                                                   aGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAATCTGATTGGCCggg                           <  2:324328/150‑1 (MQ=255)
                                                                                                                       ccTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCAc                       <  1:1260796/150‑1 (MQ=255)
                                                                                                                       ccTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCAc                       <  1:614019/150‑1 (MQ=255)
                                                                                                                        cTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCAcc                      <  1:1124439/150‑1 (MQ=255)
                                                                                                                        cTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCAcc                      >  2:231779/1‑150 (MQ=255)
                                                                                                                        atggaaaaaGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCAcc                      >  2:167504/7‑150 (MQ=255)
                                                                                                                         tGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGtt                                                                                              <  1:821071/77‑1 (MQ=255)
                                                                                                                         tGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGtt                                                                                              >  2:821071/1‑77 (MQ=255)
                                                                                                                             cAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGctct                 >  2:376519/1‑150 (MQ=255)
                                                                                                                                     ccGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGtt                                                                                              >  1:308185/1‑65 (MQ=255)
                                                                                                                                     ccGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGtt                                                                                              <  2:308185/65‑1 (MQ=255)
                                                                                                                                     ccGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGc                               <  2:723383/128‑1 (MQ=255)
                                                                                                                                     ccGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGc                               >  1:723383/1‑128 (MQ=255)
                                                                                                                                     ccGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGAACCGCTCTGCAGCAgg         >  2:389151/1‑150 (MQ=255)
                                                                                                                                      cGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGct                   <  2:698350/139‑1 (MQ=255)
                                                                                                                                      cGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGct                   >  1:698350/1‑139 (MQ=255)
                                                                                                                                           tCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGTATCTCg   >  2:163404/1‑150 (MQ=255)
                                                                                                                                             cACGGCCGTAGGCGTGGGCCGCGTGGGCGGCGCTGTCGTTACCGGCCAGCACGGTGTCGACCAAGGCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGa  <  2:268163/150‑1 (MQ=255)
                                                                                                                                              gcGGCCGACGCCTTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTgg                                                   >  2:653675/2‑100 (MQ=255)
                                                                                                                                              gcGGCCGACGCCTTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTgg                                                   <  1:653675/99‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                            
CACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGA  >  NC_002947/5688227‑5688516

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: