Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,933,717 A→C 5.8% L27R (CTC→CGC)  minD ← MinC‑MinD‑MinE system ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,933,7170AC 5.8% 198.5 / 3.4 69L27R (CTC→CGC) minDMinC‑MinD‑MinE system ATPase
Reads supporting (aligned to +/- strand):  ref base A (35/30);  new base C (2/2);  total (37/32)
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.16e-01

CCTTGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATAGCTAAATAAAGGC  >  NC_002947/1933568‑1933864
                                                                                                                                                     |                                                                                                                                                   
ccTTGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGa                                                                                                                                                     >  2:1378812/1‑150 (MQ=255)
 cTTGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAg                                                                                                                                                    <  1:1424439/150‑1 (MQ=255)
  ttGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAgg                                                                                                                                                   >  1:902533/1‑150 (MQ=255)
   tGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGc                                                                                                                                                  <  1:597231/150‑1 (MQ=255)
      tCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCgg                                                                                                                                               >  1:1187771/1‑150 (MQ=255)
                ttGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATggc                                                                                                                                     >  2:73419/1‑150 (MQ=255)
                           cTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCTCGCACTGCGCGGccaca                                                                                                                                              >  2:48374/1‑127 (MQ=255)
                           cTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCTCGCACTGCGCGGccaca                                                                                                                                              <  1:48374/130‑4 (MQ=255)
                             tCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTCGCCGCTCGGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTg                                                                                                                        <  1:1549462/150‑1 (MQ=255)
                               gCCGTTGAAAACGTTGACGAATTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGt                                                                                                                      >  2:1302591/1‑150 (MQ=255)
                                       ccACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCg                                                                                                                                                >  1:144610/1‑116 (MQ=255)
                                       ccACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCg                                                                                                                                                <  2:144610/116‑1 (MQ=255)
                                        ctcGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCcac                                                                                                             >  1:1382958/3‑150 (MQ=255)
                                        cACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGc                                                                                                                 <  1:81062/146‑1 (MQ=255)
                                        cACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGc                                                                                                                 >  2:81062/1‑146 (MQ=255)
                                        cACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCcac                                                                                                             >  1:211716/1‑150 (MQ=255)
                                         aCGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCcaca                                                                                                            >  2:579792/1‑150 (MQ=255)
                                                cGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGccaa                                                                                                              >  1:1303884/1‑139 (MQ=255)
                                                cGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCca                                                                                                              <  2:1303884/141‑1 (MQ=255)
                                                 gAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCtt                                                                                                    <  1:1013440/150‑1 (MQ=255)
                                                 gAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCtt                                                                                                    >  2:1013440/1‑150 (MQ=255)
                                                   aGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGc                                                                                                  <  2:1382958/150‑1 (MQ=255)
                                                    gTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGcc                                                                                                 <  2:902533/150‑1 (MQ=255)
                                                        tACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGaa                                                                                             >  2:1242641/1‑150 (MQ=255)
                                                           accacGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTa                                                                                          >  1:928285/1‑150 (MQ=255)
                                                             cacGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATggcg                                                                                                                                    >  2:753371/1‑106 (MQ=255)
                                                             cacGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAg                                                                                            >  2:324125/1‑146 (MQ=255)
                                                             cacGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAg                                                                                            <  1:324125/146‑1 (MQ=255)
                                                             cacGCGGCGTTCGCAGCCCATGAGCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATggcg                                                                                                                                    <  1:753371/106‑1 (MQ=255)
                                                              acGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAg                                                                                                                                                    >  2:1414868/1‑89 (MQ=255)
                                                              acGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAg                                                                                                                                                    <  1:1414868/89‑1 (MQ=255)
                                                                      ttCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAAt                                                                                                                                        <  1:727727/93‑1 (MQ=255)
                                                                      ttCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAAt                                                                                                                                        >  2:727727/1‑93 (MQ=255)
                                                                          cAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTg                                                                           <  1:1302591/150‑1 (MQ=255)
                                                                           aGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTgg                                                                          <  1:1378812/150‑1 (MQ=255)
                                                                            gCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGc                                                                         <  2:928285/150‑1 (MQ=255)
                                                                                aTGATCAGGTCGAGGTTACGCAGGCCCACGTCGACGCCCAAAACATCGACAAACAGATCGGAAAGTTCGAAGCCTGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAg                                                                     >  1:444789/1‑150 (MQ=255)
                                                                                    tCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTgg                                                                                                                             >  1:1539275/1‑90 (MQ=255)
                                                                                    tCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTgg                                                                                                                             <  2:1539275/90‑1 (MQ=255)
                                                                                     cAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGtt                                                                >  2:34266/1‑150 (MQ=255)
                                                                                        gTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGTCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGc                                                                         <  1:1342234/138‑1 (MQ=255)
                                                                                        gTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGc                                                                         >  2:1342234/1‑138 (MQ=255)
                                                                                                 aCGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTaaaa                                                    >  2:369258/1‑150 (MQ=255)
                                                                                                       gCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTACCCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTc                                              <  1:1242641/150‑1 (MQ=255)
                                                                                                         ccACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCgg                                            >  1:265186/1‑150 (MQ=255)
                                                                                                              tCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGacac                                       <  1:73419/150‑1 (MQ=255)
                                                                                                               cGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACt                                      >  2:822952/1‑150 (MQ=255)
                                                                                                                 aaGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGc                                                                         >  1:422170/1‑113 (MQ=255)
                                                                                                                 aaGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGc                                                                         <  2:422170/113‑1 (MQ=255)
                                                                                                                    tCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATc                                                                              <  2:1225791/105‑1 (MQ=255)
                                                                                                                    tCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATc                                                                              >  1:1225791/1‑105 (MQ=255)
                                                                                                                     cGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTccc                                <  1:34266/150‑1 (MQ=255)
                                                                                                                      gACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCt                               <  2:444789/150‑1 (MQ=255)
                                                                                                                      gACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCt                               >  1:1214655/1‑150 (MQ=255)
                                                                                                                            gACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGccc                                                                                                               <  1:1125717/64‑1 (MQ=255)
                                                                                                                            gACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGccc                                                                                                               >  2:1125717/1‑64 (MQ=255)
                                                                                                                                    tgtgGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAATGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGaaa                 <  2:1187771/150‑1 (MQ=255)
                                                                                                                                     ccggCCTCGCACTGCGCGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAAt                >  1:637930/3‑150 (MQ=255)
                                                                                                                                        gCCGCGCAGTGCGAGGCCGGTACCAGTGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATAGc             >  1:481846/1‑150 (MQ=255)
                                                                                                                                         ccGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCggg                                           >  1:1082117/1‑119 (MQ=255)
                                                                                                                                         ccGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCggg                                           <  2:1082117/119‑1 (MQ=255)
                                                                                                                                          cgcgCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTgg                                                                          <  2:702049/87‑1 (MQ=255)
                                                                                                                                          cgcgCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTgg                                                                          <  2:701503/87‑1 (MQ=255)
                                                                                                                                          cgcgCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTgg                                                                          >  1:702049/1‑87 (MQ=255)
                                                                                                                                          cgcgCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTgg                                                                          >  1:701503/1‑87 (MQ=255)
                                                                                                                                           gcgcAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATAGCtaa          >  1:594231/1‑150 (MQ=255)
                                                                                                                                               actGCGCGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATAGCtaaataa      <  2:637930/148‑1 (MQ=255)
                                                                                                                                                 tGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCg                                             <  1:806683/109‑1 (MQ=255)
                                                                                                                                                 tGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCg                                             >  2:806683/1‑109 (MQ=255)
                                                                                                                                                   cGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACTAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATATCTAAATAAAGGc  >  2:667905/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCTTGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATAGCTAAATAAAGGC  >  NC_002947/1933568‑1933864

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