Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,688,380 G→T 8.2% D214E (GAC→GAA pilJ ← twitching motility protein PilJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,688,3800GT 8.2% 224.5 / 4.6 72D214E (GAC→GAApilJtwitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base G (32/34);  new base T (2/4);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTA  >  NC_002947/5688234‑5688529
                                                                                                                                                   |                                                                                                                                                     
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                                                                                                            <  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                                                                                            <  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                                                                                                                                <  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)
                                                                            ggATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCt                                                                        <  1:1388620/150‑1 (MQ=255)
                                                                             gATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTg                                                                       <  1:1368482/150‑1 (MQ=255)
                                                                             gATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTg                                                                       <  1:995741/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:879110/150‑1 (MQ=255)
                                                                                                          ccAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCg                                          >  2:1260796/1‑150 (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                                            >  1:723383/1‑128 (MQ=255)
                                                                                                                              ccGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGc                                            <  2:723383/128‑1 (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)
                                                                                                                                                   gTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGACCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTa  >  1:954196/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGGCGTCGGCCGCCTGGGCGGCGCTGTCGTCACCGGCCAGCACGGTGTTGACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTA  >  NC_002947/5688234‑5688529

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: