Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,143,685 A→T 5.3% Q141L (CAG→CTG)  PP_1899 → ExbD/TolR family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,143,6850AT 5.3% 300.1 / 7.2 95Q141L (CAG→CTG) PP_1899ExbD/TolR family protein
Reads supporting (aligned to +/- strand):  ref base A (45/45);  new base T (2/3);  total (47/48)
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.97e-01

TAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAGGGGT  >  NC_002947/2143538‑2143823
                                                                                                                                                     |                                                                                                                                          
tAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAgg                                                                                                                                          >  1:458060/1‑150 (MQ=255)
    gcgcGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCTGCCCCGGTTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGTGCC‑CTGGGcc                                                                                                                                       <  1:735935/150‑2 (MQ=255)
      gcgAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCATCAAGGCGTGCC‑CTGGGccgc                                                                                                                                     >  1:662470/1‑147 (MQ=255)
           tCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGa                                                                                                                               <  2:318731/150‑1 (MQ=255)
            cGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGAt                                                                                                                              <  1:2182349/150‑1 (MQ=255)
             gggCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATg                                                                                                                             <  1:1616682/150‑1 (MQ=255)
             gggCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATg                                                                                                                             >  1:1617263/1‑150 (MQ=255)
             gggCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCATCAAGGCGTGCC‑CTGGGC‑CGCCtcgacg                                                                                                                             <  2:662470/150‑6 (MQ=255)
              ggcggcGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATgg                                                                                                                            >  2:1864960/1‑150 (MQ=255)
               gcggcgACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGc                                                                                                                           >  1:1349921/1‑150 (MQ=255)
                  gcgACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCttt                                                                                                                        <  2:1056024/150‑1 (MQ=255)
                  gcgACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCttt                                                                                                                        <  1:391586/150‑1 (MQ=255)
                     acaacaAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGc                                                                                                                     <  2:153753/150‑1 (MQ=255)
                     acaacaAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGc                                                                                                                     <  2:577330/150‑1 (MQ=255)
                      caacaaGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGcc                                                                                                                    <  2:596415/150‑1 (MQ=255)
                       aacaaGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCtt                                                                                                                                 >  1:983590/1‑136 (MQ=255)
                       aacaaGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCtt                                                                                                                                 <  2:983590/136‑1 (MQ=255)
                          aaGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATc                                                                                                                >  2:1901197/1‑150 (MQ=255)
                            gTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGc                                                                                                              >  2:134093/1‑150 (MQ=255)
                             ttGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGTG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCt                                                                                                             >  2:419259/1‑150 (MQ=255)
                              tGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCtt                                                                                                            <  1:1129444/150‑1 (MQ=255)
                               gCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTg                                                                                                           <  1:1864960/150‑1 (MQ=255)
                                  gCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTc                                                                                                        >  1:1970077/1‑150 (MQ=255)
                                       ccATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTGAGCCAGATGTGCGT‑ATCAGCGCTGACTCG‑GCCTCACGGCGCGCCTTGATGGCTTTCGCCGCTCGCTTGCTCGCGGc                                                                                                   >  1:274124/1‑150 (MQ=255)
                                       ccATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGc                                                                                                   >  1:1051320/1‑150 (MQ=255)
                                        cATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGcc                                                                                                  >  2:830658/1‑150 (MQ=255)
                                         aTCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCt                                                                                                 <  2:1617263/150‑1 (MQ=255)
                                          tCAGCGCCGACGGCAAGACCCCGCATCCGGCTGTGTGCACCGGAAGGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTg                                                                                                <  1:521024/150‑1 (MQ=255)
                                             gcgcCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTa                                                                                             <  1:1901197/150‑1 (MQ=255)
                                                    cGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCggg                                                                                      <  1:1628866/150‑1 (MQ=255)
                                                     ggCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGc                                                                                     >  2:1118398/1‑150 (MQ=255)
                                                       cAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCAt                                                                                   >  2:1563445/1‑150 (MQ=255)
                                                        aaGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATc                                                                                  >  2:1400611/1‑150 (MQ=255)
                                                           aCCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATg                                                                                                                             <  2:1529850/104‑1 (MQ=255)
                                                           aCCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATg                                                                                                                             >  1:1529850/1‑104 (MQ=255)
                                                              ccGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCggcgc                                                                            >  2:398084/1‑150 (MQ=255)
                                                               cGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCggcgct                                                                           <  1:1047564/150‑1 (MQ=255)
                                                               cGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCggcgct                                                                           >  2:2233354/1‑150 (MQ=255)
                                                                     aGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCggcgctggcgct                                                                     >  2:2155701/1‑150 (MQ=255)
                                                                          gtgGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCtgctgc                                                                >  2:1767912/1‑150 (MQ=255)
                                                                             gTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTc                                                             <  2:342178/150‑1 (MQ=255)
                                                                              tCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGTGCC‑CTGGGc                                                                                                                                        >  1:2191630/1‑75 (MQ=255)
                                                                              tCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGTGCC‑CTGGGc                                                                                                                                        <  2:2191630/75‑1 (MQ=255)
                                                                              tCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTcc                                                            <  1:1956379/150‑1 (MQ=255)
                                                                              tCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATAAATACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTcc                                                            >  1:1333688/1‑150 (MQ=255)
                                                                                acagCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGGCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGc                                                          >  1:1760707/4‑150 (MQ=255)
                                                                                    cAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGTG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGa                                                      <  1:419259/150‑1 (MQ=255)
                                                                                    cAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTg                                                                                                <  1:1382843/108‑1 (MQ=255)
                                                                                    cAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTg                                                                                                >  2:1382843/1‑108 (MQ=255)
                                                                                     aaTGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAg                                                     <  1:398084/150‑1 (MQ=255)
                                                                                       tGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGc                                                   <  1:541198/150‑1 (MQ=255)
                                                                                         gATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCt                                                 <  2:782665/150‑1 (MQ=255)
                                                                                          aTGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCtt                                                >  1:1217834/1‑150 (MQ=255)
                                                                                          aTGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCtt                                                <  1:134093/150‑1 (MQ=255)
                                                                                          aTGCTGCCGGCAAGCTCGGTGTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGCCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCtt                                                <  2:130152/150‑1 (MQ=255)
                                                                                              tgcCGACAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTAt                                            <  1:1433437/150‑1 (MQ=255)
                                                                                                   gCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATcgccg                                       >  2:625622/1‑150 (MQ=255)
                                                                                                     aaGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCcgcg                                     >  1:18879/1‑150 (MQ=255)
                                                                                                      aGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGt                                    <  1:1790995/150‑1 (MQ=255)
                                                                                                       gCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTg                                   <  2:1970077/150‑1 (MQ=255)
                                                                                                        cTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCAAGCCTTGCTGGCTTTCGCCGATCGCGTGCTCGCTGCCTTGGATGCCGGGCTTCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTgg                                  >  2:851813/1‑150 (MQ=255)
                                                                                                         tCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGt                                 >  1:761023/1‑150 (MQ=255)
                                                                                                         tCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGt                                 >  1:1135247/1‑150 (MQ=255)
                                                                                                          cGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTc                                <  1:2110940/150‑1 (MQ=255)
                                                                                                          cGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGAGCATAAAATGCTTGTGGTCCTTCCTCCTCTGGAGGCCTTGTATCGCCGCGTGGTc                                <  1:407208/150‑1 (MQ=255)
                                                                                                              ttCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCAcgc                            <  1:128129/150‑1 (MQ=255)
                                                                                                                 aGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCggcgctggcgct                                                                     <  1:1000527/106‑1 (MQ=255)
                                                                                                                 aGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCggcgctggcgct                                                                     >  2:1000527/1‑106 (MQ=255)
                                                                                                                  gCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCaa                        >  1:214076/1‑150 (MQ=255)
                                                                                                                  gCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTAGATCGGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCaa                        >  1:1874040/1‑150 (MQ=255)
                                                                                                                  gCCACTTGCGCAGGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCaa                        <  2:1333688/150‑1 (MQ=255)
                                                                                                                   ccAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAg                       >  2:1030606/1‑150 (MQ=255)
                                                                                                                    cAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGc                      <  1:1544441/150‑1 (MQ=255)
                                                                                                                     aGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCg                                                                                                       >  1:2124680/1‑68 (MQ=255)
                                                                                                                     aGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCg                                                                                                       <  2:2124680/68‑1 (MQ=255)
                                                                                                                      gTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCgg                    <  2:1473088/150‑1 (MQ=255)
                                                                                                                          cgcATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGAGCAAGCGGGCAc                >  2:2246574/1‑150 (MQ=255)
                                                                                                                           gcATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACg               >  2:1424382/1‑150 (MQ=255)
                                                                                                                           gcATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACg               >  2:419504/1‑150 (MQ=255)
                                                                                                                            cATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGt              >  1:118468/1‑150 (MQ=255)
                                                                                                                             aTGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGa                                                      >  2:2003199/1‑109 (MQ=255)
                                                                                                                             aTGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGa                                                      <  1:2003199/109‑1 (MQ=255)
                                                                                                                                 ccaccGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCt         <  2:1020022/150‑1 (MQ=255)
                                                                                                                                  caccGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTa        >  2:561830/1‑150 (MQ=255)
                                                                                                                                    ccGTCGAGGCG‑GCC‑CAGGGCACGGCTGGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAg      <  2:135071/150‑1 (MQ=255)
                                                                                                                                      atCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAggg    >  1:1864665/2‑150 (MQ=255)
                                                                                                                                       tCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGt                                    <  2:120578/117‑1 (MQ=255)
                                                                                                                                       tCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGt                                    >  1:120578/1‑117 (MQ=255)
                                                                                                                                       tCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCg                     >  1:2095002/1‑132 (MQ=255)
                                                                                                                                       tCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCg                     <  2:2095002/132‑1 (MQ=255)
                                                                                                                                       tCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAgggg   <  1:830658/150‑1 (MQ=255)
                                                                                                                                       tCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGACGGTGCGGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAgggg   >  1:1352074/1‑150 (MQ=255)
                                                                                                                                        cGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAGGGGt  <  1:1030606/150‑1 (MQ=255)
                                                                                                                                                  c‑cAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAg                                                     >  2:34121/1‑90 (MQ=255)
                                                                                                                                                  c‑cAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAg                                                     <  1:34121/90‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
TAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCG‑GCC‑CAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAGGGGT  >  NC_002947/2143538‑2143823

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: ‑