Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,263,401 Δ1 bp 6.2% coding (1189/1635 nt) PP_2861 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,263,4010T. 6.2% 352.9 / 16.8 112coding (1189/1635 nt)PP_2861methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (55/50);  new base . (4/3);  total (59/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 112 total aligned reads displayed.

TTGACCGCAGCGTCCGACTGACTCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGCGGTCGATCTGCGC  >  NC_002947/3263253‑3263548
                                                                                                                                                      |                                                                                                                                                   
ttGACCGCAGCGTCCGACTGACTCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATg                                                                                                                                                    >  2:1275854/1‑150 (MQ=255)
   aCCGCAGCGTCCGACTGACTCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATggcg                                                                                                                                                 >  1:1875590/1‑150 (MQ=255)
     cGCAGCGTCCGACTGACTCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCACCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATggcggc                                                                                                                                               >  2:1074503/1‑150 (MQ=255)
       cAGCGTCCGACTGACTCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCGCGGCCTTGCTCGCCAGCCCGTGCCGCCATCGAGGCGGCac                                                                                                                                             <  2:398547/150‑2 (MQ=38)
       cAGCGTCCGACTGACTCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCATCGAGGCGGCac                                                                                                                                             <  1:2190281/150‑2 (MQ=255)
       cAGCGTCCGACTGACTCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCAt                                                                                                                                             >  1:1747624/1‑150 (MQ=255)
        aGCGTCCGACTGACTCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCAtt                                                                                                                                            <  1:896866/150‑1 (MQ=255)
        aGCGTCCGACTGACTCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCAtg                                                                                                                                            <  1:202049/150‑2 (MQ=255)
              cGACTGACTCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAgcgc                                                                                                                                      >  2:1841202/1‑150 (MQ=255)
               gactgactCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGcc                                                                                                                                     >  1:935726/1‑150 (MQ=255)
               gactgactCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGcc                                                                                                                                     >  1:1820691/1‑150 (MQ=255)
                 ctgactCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCcag                                                                                                                                   <  2:1747624/150‑1 (MQ=255)
                   gactCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCcagca                                                                                                                                 <  2:1820691/150‑1 (MQ=255)
                    actCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCATGCCGCCATCGAGGCGGCacgg                                                                                                                                           <  1:863391/139‑4 (MQ=38)
                    actCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCATGCCGCCATCGAGGCGGCacgg                                                                                                                                           >  2:863391/1‑136 (MQ=38)
                     ctctGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATTGCGGCATTCGGCGCCAGCAgg                                                                                                                               <  2:398884/150‑1 (MQ=255)
                        tGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTg                                                                                                                            >  1:1973493/1‑150 (MQ=255)
                          cAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGt                                                                                                                          <  1:1694088/150‑1 (MQ=255)
                              tGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGTCCCTCGTCGCTGGCCACGCTAACCGGAGAACACGCATGCGCAACCCGTACCGCCTCGATGGCGGCATTCAGCGCAAGCGGGTGGGTCTGt                                                                                                                      >  1:782706/1‑150 (MQ=255)
                                 tCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCg                                                                                                                   >  2:200982/1‑150 (MQ=255)
                                 tCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCg                                                                                                                   >  2:1581547/1‑150 (MQ=255)
                                   gATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATggcgg                                                                                                                                                <  1:1726626/119‑1 (MQ=255)
                                   gATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATggcgg                                                                                                                                                >  2:1726626/1‑119 (MQ=255)
                                      cATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGTCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAAt                                                                                                              <  2:304335/150‑1 (MQ=255)
                                          gTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGaa                                                                                                          <  1:556466/150‑1 (MQ=255)
                                           tATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAAc                                                                                                         <  1:1074503/150‑1 (MQ=255)
                                           tATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAAc                                                                                                         >  1:1479598/1‑150 (MQ=255)
                                                aTCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATc                                                                                                    <  2:1973493/150‑1 (MQ=255)
                                                aTCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGACGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATc                                                                                                    <  1:53238/150‑1 (MQ=255)
                                                  ctcGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTGGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCAc                                                                                                  <  1:1692654/150‑1 (MQ=255)
                                                     ggCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTc                                                                                               >  1:1491481/1‑150 (MQ=255)
                                                       ccGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTg                                                                                                                            <  2:1788977/119‑1 (MQ=255)
                                                       ccGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTg                                                                                                                            >  1:1788977/1‑119 (MQ=255)
                                                       ccGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCa                                                                                             <  2:731241/150‑1 (MQ=255)
                                                         gTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCg                                                                                                                   >  2:605772/1‑126 (MQ=255)
                                                         gTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCg                                                                                                                   <  1:605772/126‑1 (MQ=255)
                                                          tCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCa                                                                                          >  1:447058/1‑150 (MQ=255)
                                                             aTTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATc                                                                                                    <  2:1737024/137‑1 (MQ=255)
                                                             aTTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATc                                                                                                    >  1:1737024/1‑137 (MQ=255)
                                                                 tgCTGTACGCTGAGCCAACAGGCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATAGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACgccgcc                                                                                   >  2:106686/1‑150 (MQ=255)
                                                                  gCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGtt                                                                                                                             >  2:451475/1‑107 (MQ=255)
                                                                  gCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGtt                                                                                                                             <  1:451475/107‑1 (MQ=255)
                                                                   cTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGa                                                                                 >  1:2104588/1‑150 (MQ=255)
                                                                    tGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGtt                                                                                                                     >  2:1671806/1‑113 (MQ=255)
                                                                    tGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGAt                                                                                >  2:2176576/1‑150 (MQ=255)
                                                                    tGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGAt                                                                                >  2:1629029/1‑150 (MQ=255)
                                                                        cGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATgcgc                                                                            <  1:517307/150‑1 (MQ=255)
                                                                         gCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGcgcg                                                                           >  1:1205273/1‑150 (MQ=255)
                                                                          cTGAGCCAACAGTCGCACCTCGTCGGCGAGCACTGCAAATCTCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGt                                                                          <  2:444769/150‑1 (MQ=255)
                                                                           tGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTg                                                                         <  2:1875590/150‑1 (MQ=255)
                                                                               ccAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCtgt                                                                     >  1:616328/1‑150 (MQ=255)
                                                                               ccAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCtgt                                                                     <  1:1275854/150‑1 (MQ=255)
                                                                                cAACAGTCGCACCTCGTCGGCGACCACCGCAAAGCCCCGGCCTTGCTCGCCAGGCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGATCGAATCACGTCCAGCACGCCGCCGATGCGCGTGCtgtg                                                                    <  2:782706/150‑1 (MQ=255)
                                                                                    aGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCTAGGCGCGTGCTGTggc                                                                  <  2:1192256/148‑1 (MQ=255)
                                                                                    aGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTggc                                                                  >  1:1192256/1‑148 (MQ=255)
                                                                                        gCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAAt                                                                                                     >  1:461206/1‑109 (MQ=255)
                                                                                        gCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCa                                                            >  2:628697/1‑150 (MQ=255)
                                                                                          aCCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCCTGCTGTGGCCGGCCAgg                                                          >  2:2162540/1‑150 (MQ=255)
                                                                                              cgtcgGCGACCACGGCCTCCCCCCGGCCTTGTTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGAGGCGCGTGCTGTGGCCGGCCAGGTcgc                                                      <  1:736180/150‑1 (MQ=37)
                                                                                                    cGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCAGCTCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTa                                                >  1:23545/1‑150 (MQ=255)
                                                                                                     gACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTAc                                               >  1:2362429/1‑150 (MQ=255)
                                                                                                        cACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTg                                            >  1:402767/1‑150 (MQ=255)
                                                                                                            gCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCtgtg                                                                    >  1:1599080/1‑122 (MQ=255)
                                                                                                            gCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCtgtg                                                                    <  2:1599080/122‑1 (MQ=255)
                                                                                                            gCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTAc                                               <  1:113923/143‑1 (MQ=255)
                                                                                                              aaaGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGc                                      <  2:2362429/150‑1 (MQ=255)
                                                                                                                 gCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGACAGGTGGCGGATTACCTGTACGGCCtg                                   <  2:284226/150‑2 (MQ=255)
                                                                                                                   cccGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGGCTGTTCGGCAATGGAACTAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGGTTACCTGTACGGCCTcgt                                 >  1:1101687/1‑148 (MQ=255)
                                                                                                                    ccGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCCTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATc                                <  1:2162540/150‑1 (MQ=255)
                                                                                                                     cGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTggcc                                                                 >  1:623694/1‑116 (MQ=255)
                                                                                                                     cGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTggcc                                                                 <  2:623694/116‑1 (MQ=255)
                                                                                                                       gCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGAt                             >  2:2326667/1‑150 (MQ=255)
                                                                                                                          ttGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGt                          >  1:1113657/1‑150 (MQ=255)
                                                                                                                          ttGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTAAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGGCGCGGATTACCTGTACGGCCTCATCGATGGt                          >  1:48686/1‑150 (MQ=255)
                                                                                                                            gCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCa                        <  2:402767/150‑1 (MQ=255)
                                                                                                                              tCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGc                      >  1:1752299/1‑150 (MQ=255)
                                                                                                                              tCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGCCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGc                      >  1:221953/1‑150 (MQ=255)
                                                                                                                              tCGCCAGCCCGTGCCGCCTCGATGGCAGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCGCGCCGCCGATGCGCGTGCTATGGCCGGCCAGGACGCGGATTACCGGTACGGCCTCATCGATGGTCAGc                      >  2:1934698/1‑150 (MQ=39)
                                                                                                                                gCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATgg                           >  2:462179/1‑143 (MQ=255)
                                                                                                                                gCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGa                    >  2:2149840/1‑150 (MQ=255)
                                                                                                                                    aGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAg                 >  1:2043308/1‑150 (MQ=255)
                                                                                                                                     gCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCa                        >  2:2218450/1‑142 (MQ=255)
                                                                                                                                      tgaatGCCGCCATCGAGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGc               >  2:409736/5‑150 (MQ=255)
                                                                                                                                       ccGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGcc                                     <  1:1319631/127‑1 (MQ=255)
                                                                                                                                       ccGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGcc                                     >  2:1319631/1‑127 (MQ=255)
                                                                                                                                       ccGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGCg              >  1:1638270/1‑150 (MQ=255)
                                                                                                                                       ccGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGCCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGCg              <  1:624477/150‑1 (MQ=255)
                                                                                                                                       ccGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCAGCCCAGGGCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGCg              >  1:1143498/1‑150 (MQ=255)
                                                                                                                                        cGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGCgg             <  1:200982/150‑1 (MQ=255)
                                                                                                                                        aatGCCGCCATCGAGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGCgg             >  1:1668934/3‑150 (MQ=255)
                                                                                                                                         gTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTcgc                                                      <  2:228831/108‑1 (MQ=255)
                                                                                                                                         gTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTcgc                                                      >  1:228831/1‑108 (MQ=255)
                                                                                                                                         gTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGCGGt            <  1:106686/150‑1 (MQ=255)
                                                                                                                                         gTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGCGGt            >  2:75814/1‑150 (MQ=255)
                                                                                                                                         gTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGCGGt            <  1:1581547/150‑1 (MQ=255)
                                                                                                                                         gTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGCGGt            <  1:2176576/150‑1 (MQ=255)
                                                                                                                                          tGCCGCCATCGAGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCAAGCACGCCGCCGATGCGCGTGCTGTGTCCGGCCAGGTCGCGGATTACCTGTCCGGCCTCATAGGTGGTCAGGGACAGGCGGTc           >  1:825397/1‑150 (MQ=255)
                                                                                                                                                  tCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGCGGTCGATCTgc    <  1:102009/150‑1 (MQ=255)
                                                                                                                                                  tCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAACACGCCCCCGATCCGCCTACTTTGGCCGGCTCCAGCGCGTATTACCTGTACGGCCTCAGCGATGGTCAGCGACAGGCGGTCGATCTgc    <  2:1479598/150‑1 (MQ=37)
                                                                                                                                                    gATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGCGGTCGATCTgcgc  <  2:1113657/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TTGACCGCAGCGTCCGACTGACTCTGCAGGTGCTCGATCATGGTATGAATCTCGGCCGTCGATTGTGCTGTACGCTGAGCCAACAGTCGCACCTCGTCGGCGACCACGGCAAAGCCCCGGCCTTGCTCGCCAGCCCGTGCCGCCTCGATGGCGGCATTCAGCGCCAGCAGGTTGGTCTGTTCGGCAATGGAACGAATCACGTCCAGCACGCCGCCGATGCGCGTGCTGTGGCCGGCCAGGTCGCGGATTACCTGTACGGCCTCATCGATGGTCAGCGACAGGCGGTCGATCTGCGC  >  NC_002947/3263253‑3263548

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: