Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,143,680:1 +T 5.4% coding (417/435 nt) 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,6801.T 5.4% 379.2 / 14.3 93coding (417/435 nt)PP_1899ExbD/TolR family protein
Reads supporting (aligned to +/- strand):  ref base . (44/44);  new base T (2/3);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

GAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAGGGGT  >  NC_002947/2143531‑2143823
                                                                                                                                                      |                                                                                                                                                
gAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAggcg                                                                                                                                                   <  2:1520167/150‑1 (MQ=255)
       tAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAgg                                                                                                                                          >  1:458060/1‑150 (MQ=255)
           gcgcGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCTGCCCCGGTTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGTGCCCTGGGcc                                                                                                                                       <  1:735935/150‑2 (MQ=255)
             gcgAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCATCAAGGCGTGCCCTGGGccgc                                                                                                                                     >  1:662470/1‑147 (MQ=255)
                  tCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGa                                                                                                                               <  2:318731/150‑1 (MQ=255)
                   cGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGAt                                                                                                                              <  1:2182349/150‑1 (MQ=255)
                    gggCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATg                                                                                                                             <  1:1616682/150‑1 (MQ=255)
                    gggCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATg                                                                                                                             >  1:1617263/1‑150 (MQ=255)
                    gggCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCATCAAGGCGTGCCCTGGGCCGCCtcgacg                                                                                                                             <  2:662470/150‑6 (MQ=255)
                     ggcggcGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATgg                                                                                                                            >  2:1864960/1‑150 (MQ=255)
                      gcggcgACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGc                                                                                                                           >  1:1349921/1‑150 (MQ=255)
                         gcgACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCttt                                                                                                                        <  1:391586/150‑1 (MQ=255)
                         gcgACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCttt                                                                                                                        <  2:1056024/150‑1 (MQ=255)
                            acaacaAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGc                                                                                                                     <  2:577330/150‑1 (MQ=255)
                            acaacaAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGc                                                                                                                     <  2:153753/150‑1 (MQ=255)
                             caacaaGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGcc                                                                                                                    <  2:596415/150‑1 (MQ=255)
                              aacaaGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCtt                                                                                                                                 <  2:983590/136‑1 (MQ=255)
                              aacaaGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCtt                                                                                                                                 >  1:983590/1‑136 (MQ=255)
                                 aaGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATc                                                                                                                >  2:1901197/1‑150 (MQ=255)
                                   gTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGc                                                                                                              >  2:134093/1‑150 (MQ=255)
                                    ttGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGTGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCt                                                                                                             >  2:419259/1‑150 (MQ=255)
                                     tGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCtt                                                                                                            <  1:1129444/150‑1 (MQ=255)
                                      gCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTg                                                                                                           <  1:1864960/150‑1 (MQ=255)
                                         gCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTc                                                                                                        >  1:1970077/1‑150 (MQ=255)
                                              ccATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTGAGCCAGATGTGCGTATCAGCGCTGACTCGGCCTCACGGCGCGCCTTGATGGCTTTCGCCGCTCGCTTGCTCGCGGc                                                                                                   >  1:274124/1‑150 (MQ=255)
                                              ccATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGc                                                                                                   >  1:1051320/1‑150 (MQ=255)
                                               cATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGcc                                                                                                  >  2:830658/1‑150 (MQ=255)
                                                aTCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCt                                                                                                 <  2:1617263/150‑1 (MQ=255)
                                                 tCAGCGCCGACGGCAAGACCCCGCATCCGGCTGTGTGCACCGGAAGGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTg                                                                                                <  1:521024/150‑1 (MQ=255)
                                                    gcgcCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTa                                                                                             <  1:1901197/150‑1 (MQ=255)
                                                           cGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCggg                                                                                      <  1:1628866/150‑1 (MQ=255)
                                                            ggCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGc                                                                                     >  2:1118398/1‑150 (MQ=255)
                                                              cAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCAt                                                                                   >  2:1563445/1‑150 (MQ=255)
                                                               aaGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATc                                                                                  >  2:1400611/1‑150 (MQ=255)
                                                                  aCCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATg                                                                                                                             >  1:1529850/1‑104 (MQ=255)
                                                                  aCCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATg                                                                                                                             <  2:1529850/104‑1 (MQ=255)
                                                                     ccGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCggcgc                                                                            >  2:398084/1‑150 (MQ=255)
                                                                      cGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCggcgct                                                                           <  1:1047564/150‑1 (MQ=255)
                                                                      cGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCggcgct                                                                           >  2:2233354/1‑150 (MQ=255)
                                                                            aGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCggcgctggcgct                                                                     >  2:2155701/1‑150 (MQ=255)
                                                                                 gtgGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCtgctgc                                                                >  2:1767912/1‑150 (MQ=255)
                                                                                    gTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTc                                                             <  2:342178/150‑1 (MQ=255)
                                                                                     tCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGTGCCCTGGGc                                                                                                                                        >  1:2191630/1‑75 (MQ=255)
                                                                                     tCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGTGCCCTGGGc                                                                                                                                        <  2:2191630/75‑1 (MQ=255)
                                                                                     tCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTcc                                                            <  1:1956379/150‑1 (MQ=255)
                                                                                     tCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATAAATACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTcc                                                            >  1:1333688/1‑150 (MQ=255)
                                                                                       acagCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGGCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGc                                                          >  1:1760707/4‑150 (MQ=255)
                                                                                           cAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGTGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGa                                                      <  1:419259/150‑1 (MQ=255)
                                                                                           cAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTg                                                                                                >  2:1382843/1‑108 (MQ=255)
                                                                                           cAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTg                                                                                                <  1:1382843/108‑1 (MQ=255)
                                                                                            aaTGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAg                                                     <  1:398084/150‑1 (MQ=255)
                                                                                              tGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGc                                                   <  1:541198/150‑1 (MQ=255)
                                                                                                gATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCt                                                 <  2:782665/150‑1 (MQ=255)
                                                                                                 aTGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCtt                                                >  1:1217834/1‑150 (MQ=255)
                                                                                                 aTGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCtt                                                <  1:134093/150‑1 (MQ=255)
                                                                                                 aTGCTGCCGGCAAGCTCGGTGTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGCCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCtt                                                <  2:130152/150‑1 (MQ=255)
                                                                                                     tgcCGACAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTAt                                            <  1:1433437/150‑1 (MQ=255)
                                                                                                          gCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATcgccg                                       >  2:625622/1‑150 (MQ=255)
                                                                                                            aaGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCcgcg                                     >  1:18879/1‑150 (MQ=255)
                                                                                                             aGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGt                                    <  1:1790995/150‑1 (MQ=255)
                                                                                                              gCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTg                                   <  2:1970077/150‑1 (MQ=255)
                                                                                                               cTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCAAGCCTTGCTGGCTTTCGCCGATCGCGTGCTCGCTGCCTTGGATGCCGGGCTTCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTgg                                  >  2:851813/1‑150 (MQ=255)
                                                                                                                tCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGt                                 >  1:761023/1‑150 (MQ=255)
                                                                                                                tCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGt                                 >  1:1135247/1‑150 (MQ=255)
                                                                                                                 cGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTc                                <  1:2110940/150‑1 (MQ=255)
                                                                                                                 cGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGAGCATAAAATGCTTGTGGTCCTTCCTCCTCTGGAGGCCTTGTATCGCCGCGTGGTc                                <  1:407208/150‑1 (MQ=255)
                                                                                                                     ttCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCAcgc                            <  1:128129/150‑1 (MQ=255)
                                                                                                                        aGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCggcgctggcgct                                                                     >  2:1000527/1‑106 (MQ=255)
                                                                                                                        aGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCggcgctggcgct                                                                     <  1:1000527/106‑1 (MQ=255)
                                                                                                                         gCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCaa                        >  1:214076/1‑150 (MQ=255)
                                                                                                                         gCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTAGATCGGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCaa                        >  1:1874040/1‑150 (MQ=255)
                                                                                                                         gCCACTTGCGCAGGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCaa                        <  2:1333688/150‑1 (MQ=255)
                                                                                                                          ccAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAg                       >  2:1030606/1‑150 (MQ=255)
                                                                                                                           cAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGc                      <  1:1544441/150‑1 (MQ=255)
                                                                                                                            aGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCg                                                                                                       >  1:2124680/1‑68 (MQ=255)
                                                                                                                            aGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCg                                                                                                       <  2:2124680/68‑1 (MQ=255)
                                                                                                                             gTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCgg                    <  2:1473088/150‑1 (MQ=255)
                                                                                                                                 cgcATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGAGCAAGCGGGCAc                >  2:2246574/1‑150 (MQ=255)
                                                                                                                                  gcATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACg               >  2:419504/1‑150 (MQ=255)
                                                                                                                                  gcATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACg               >  2:1424382/1‑150 (MQ=255)
                                                                                                                                   cATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGt              >  1:118468/1‑150 (MQ=255)
                                                                                                                                    aTGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGa                                                      <  1:2003199/109‑1 (MQ=255)
                                                                                                                                    aTGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGa                                                      >  2:2003199/1‑109 (MQ=255)
                                                                                                                                        ccaccGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCt         <  2:1020022/150‑1 (MQ=255)
                                                                                                                                         caccGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTa        >  2:561830/1‑150 (MQ=255)
                                                                                                                                           ccGTCGAGGCGGCCCAGGGCACGGCTGGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAg      <  2:135071/150‑1 (MQ=255)
                                                                                                                                             atCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAggg    >  1:1864665/2‑150 (MQ=255)
                                                                                                                                              tCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGt                                    <  2:120578/117‑1 (MQ=255)
                                                                                                                                              tCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGt                                    >  1:120578/1‑117 (MQ=255)
                                                                                                                                              tCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCg                     <  2:2095002/132‑1 (MQ=255)
                                                                                                                                              tCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCg                     >  1:2095002/1‑132 (MQ=255)
                                                                                                                                              tCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAgggg   <  1:830658/150‑1 (MQ=255)
                                                                                                                                              tCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGACGGTGCGGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAgggg   >  1:1352074/1‑150 (MQ=255)
                                                                                                                                               cGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAGGGGt  <  1:1030606/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAGGCGGCCCAGGGCACGCCTTGATGGCTTTCGCCGATCGCTTGCTCGCTGCCTGGTATGCCGGGCATCCGGCGCTGGCGCTGCTGCGTCCGCTGGAGGCCTTGTATCGCCGCGTGGTCACGCGCAAGCGGGCACGTTTTCTAAGGGGT  >  NC_002947/2143531‑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: