Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,933,674 A→C 6.5% 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,6740AC 6.5% 184.7 / 4.5 62L27R (CTC→CGC) minDMinC‑MinD‑MinE system ATPase
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base C (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

CTTGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATAGCTAAATAAAGG  >  NC_002947/1933526‑1933820
                                                                                                                                                     |                                                                                                                                                   
cTTGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGa                                                                                                                                                     <  6:153555/149‑1 (MQ=255)
 ttGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCG‑TCCACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCGTGTGGCCGCGCAGTGCGAg                                                                                                                                                    >  4:234307/1‑149 (MQ=255)
    aTCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGcc                                                                                                                                                 <  5:199870/149‑1 (MQ=255)
    aTCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGc                                                                                                                                                  <  8:296952/148‑1 (MQ=255)
     tCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCg                                                                                                                                                >  7:119313/1‑149 (MQ=255)
     tCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCg                                                                                                                                                >  3:102548/1‑149 (MQ=255)
       aGGGCCTGTTGCAGGTGGGCTTCGCCGGTGACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCgg                                                                                                                                               <  6:195264/148‑1 (MQ=255)
                      ttGGCTTCGCCGTTGACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCt                                                                                                                              <  6:321570/149‑1 (MQ=255)
                      ttGGCTTCGCCGTTGACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCt                                                                                                                              <  6:218306/149‑1 (MQ=255)
                          cTTCGCCGTTGACCACGTTGACGAAGTCGTTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTgg                                                                                                                            >  1:195236/1‑148 (MQ=255)
                          cTTCGCCGTTGACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTa                                                                                                                          <  7:80531/149‑1 (MQ=255)
                               ccGTTGACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGt                                                                                                                     >  8:7699/1‑149 (MQ=255)
                                  ttGACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCtt                                                                                                                  <  3:248626/149‑1 (MQ=255)
                                    gACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTg                                                                                                                 <  4:102548/148‑1 (MQ=255)
                                    gACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTg                                                                                                                 <  8:119313/148‑1 (MQ=255)
                                     aCCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGcc                                                                                                               >  2:203587/1‑149 (MQ=255)
                                          gTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCacac                                                                                                          >  2:154441/1‑149 (MQ=255)
                                                gAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGcc                                                                                                               <  7:257757/138‑1 (MQ=255)
                                                gAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGcc                                                                                                               >  8:257757/1‑138 (MQ=255)
                                                gAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCt                                                                                                    >  6:56071/1‑149 (MQ=255)
                                                gAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGCGGCCTCGCACTGCGCGGccacaa                                                                                                                                             <  3:210569/109‑5 (MQ=255)
                                                gAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGCGGCCTCGCACTGCGCGGccacaa                                                                                                                                             >  4:210569/1‑105 (MQ=255)
                                                      g‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCg                                                                                               <  5:56071/148‑1 (MQ=255)
                                                          caccacGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAg                                                                                           >  7:13554/1‑149 (MQ=255)
                                                             cacGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGgcgg‑cg                                                                                                                                >  7:268833/1‑109 (MQ=255)
                                                             cacGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGgcgg‑cg                                                                                                                                <  8:268833/109‑1 (MQ=255)
                                                                             cccATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGaa                                                                               >  2:228342/1‑142 (MQ=255)
                                                                             cccATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGc                                                                        >  6:192742/1‑149 (MQ=255)
                                                                             cccATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGc                                                                        <  6:190046/149‑1 (MQ=255)
                                                                             cccATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGCTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGaa                                                                               <  1:228342/142‑1 (MQ=255)
                                                                               cATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGcc                                                                       >  8:110327/1‑148 (MQ=255)
                                                                                aTGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGaa                                                                                            >  6:62797/1‑126 (MQ=255)
                                                                                aTGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGaa                                                                                            <  5:62797/126‑1 (MQ=255)
                                                                                aTGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGcc                                                                       >  6:3467/1‑147 (MQ=255)
                                                                                    tCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCAAACACCCCCCTTGCCGGAAGTAAACACGAGAATCTTGGCCAAGgtg                                                                 >  8:262228/1‑149 (MQ=255)
                                                                                       ggTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGttt                                                              <  1:251808/149‑1 (MQ=255)
                                                                                         tCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGT‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACAGCGAGAATCTTGGCCAAGGTGTTTCa                                                            >  4:166556/1‑149 (MQ=255)
                                                                                            aGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCAccc                                                         <  2:195236/149‑1 (MQ=255)
                                                                                                   gCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGAAACCACCAGAATCTTGGCAAAGGTGTTTCACCCCTaaaaa                                                  >  4:181003/1‑149 (MQ=255)
                                                                                                     aGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATa                                                >  7:105118/1‑149 (MQ=255)
                                                                                                         ccACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCGAg                                                                    <  4:289061/125‑1 (MQ=255)
                                                                                                         ccACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCGAg                                                                    >  3:289061/1‑125 (MQ=255)
                                                                                                         ccACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCg                                            <  1:154441/149‑1 (MQ=255)
                                                                                                          cACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCTAGGCATTTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCgg                                           <  7:262228/149‑1 (MQ=255)
                                                                                                            cGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGa                                         <  7:7699/149‑1 (MQ=255)
                                                                                                                         aaTGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTg                             >  4:216541/1‑148 (MQ=255)
                                                                                                                           tGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCtcc                            >  1:72643/1‑146 (MQ=255)
                                                                                                                           tGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCa                           >  7:85290/1‑148 (MQ=255)
                                                                                                                           tGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCa                           >  1:57451/1‑148 (MQ=255)
                                                                                                                             aCGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCaaa                         >  1:327561/1‑148 (MQ=255)
                                                                                                                               ggTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCaaaaaa                      <  2:19132/149‑1 (MQ=255)
                                                                                                                               ggTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCaaaaaa                      <  3:181003/149‑1 (MQ=255)
                                                                                                                                    tgtgGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGa                  >  8:274695/1‑148 (MQ=255)
                                                                                                                                     gtgGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGaaa                >  2:43253/1‑149 (MQ=255)
                                                                                                                                      cggCCTCGCACTGCGCGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGaaa                >  5:21489/2‑148 (MQ=255)
                                                                                                                                        gCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCCCACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAagag             >  7:293775/1‑149 (MQ=255)
                                                                                                                                         ccGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCAcccc                                                        <  8:41439/105‑1 (MQ=255)
                                                                                                                                             gcACTGCGCGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATAGCtaa         <  6:21489/148‑1 (MQ=255)
                                                                                                                                               actGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCaaaaaa                      >  4:253873/3‑133 (MQ=255)
                                                                                                                                               actGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCaaaaaa                      <  3:253873/131‑1 (MQ=255)
                                                                                                                                                 tGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATAGCtaaataaa    <  2:57451/149‑1 (MQ=255)
                                                                                                                                                   cGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATAGCTAAATAAAg   <  2:327561/148‑1 (MQ=255)
                                                                                                                                                    gAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATAGCTAAATAAAgg  >  6:248140/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTTGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCG‑TACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGG‑CGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATAGCTAAATAAAGG  >  NC_002947/1933526‑1933820

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