Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,688,378 G→C 6.9% A215G (GCC→GGC)  pilJ ← twitching motility protein PilJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,688,3780GC 6.9% 223.7 / 6.3 72A215G (GCC→GGC) pilJtwitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base C (2/3);  total (34/38)
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.67e-01

CCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGA  >  NC_002947/5688229‑5688516
                                                                                                                                                      |                                                                                                                                          
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                                                                                                                           <  2:900528/141‑1 (MQ=255)
                          tCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGggcg                                                                                                                           >  1:900528/1‑141 (MQ=255)
                          tCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTgtcg                                                                                                                  >  2:5364/1‑150 (MQ=255)
                             gCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGGCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCa                                                                                                               >  2:125788/1‑150 (MQ=255)
                                             gcgcACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGgtgt                                                                                               >  1:638912/1‑150 (MQ=255)
                                             gcgcACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGgtgt                                                                                               <  2:542465/150‑1 (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                                                                               <  1:620416/150‑1 (MQ=255)
                                                             gTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGAt                                                                               <  2:849528/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                                                                                                                   >  2:1196030/1‑96 (MQ=255)
                                                                               cTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTgtc                                                                                                                   <  1:1196030/96‑1 (MQ=255)
                                                                               cTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGc                                                                      <  1:1301695/141‑1 (MQ=255)
                                                                               cTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAg                                                             >  1:993655/1‑150 (MQ=255)
                                                                               cTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGGGTTTACCGAGCCAAGGATGCGTTCGGc                                                                      >  2:1301695/1‑141 (MQ=255)
                                                                                 ggATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCt                                                           <  1:1388620/150‑1 (MQ=255)
                                                                                  gATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTg                                                          <  1:1368482/150‑1 (MQ=255)
                                                                                  gATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTg                                                          <  1:995741/150‑1 (MQ=255)
                                                                                       tCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGACGCCTTCGgccgt                                                                                                                                  <  1:317093/73‑2 (MQ=38)
                                                                                       tCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGACGCCTTCGgccgt                                                                                                                                  >  2:317093/1‑72 (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                                                                    >  1:1369364/1‑124 (MQ=255)
                                                                                                  ttGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCa                                                                    <  2:1369364/124‑1 (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                                                                                >  1:1084068/1‑101 (MQ=255)
                                                                                                             cTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGa                                                                                <  2:1084068/101‑1 (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:614019/150‑1 (MQ=255)
                                                                                                                     ccTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCAc                       <  1:1260796/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                                                                                              >  2:821071/1‑77 (MQ=255)
                                                                                                                       tGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGtt                                                                                              <  1:821071/77‑1 (MQ=255)
                                                                                                                           cAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGctct                 >  2:376519/1‑150 (MQ=255)
                                                                                                                                   ccGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGtt                                                                                              <  2:308185/65‑1 (MQ=255)
                                                                                                                                   ccGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGtt                                                                                              >  1:308185/1‑65 (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)
                                                                                                                                                      |                                                                                                                                          
CCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGA  >  NC_002947/5688229‑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: