Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,198,338 C→T 5.2% A214A (GCC→GCT gspE → protein secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,198,3380CT 5.2% 552.3 / 12.3 174A214A (GCC→GCTgspEprotein secretion protein
Reads supporting (aligned to +/- strand):  ref base C (76/89);  new base T (6/3);  total (82/92)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 183 total aligned reads displayed.

ATGGCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGCCT  >  NC_002947/1198189‑1198487
                                                                                                                                                     |                                                                                                                                                     
aTGGCCCGCCTGGACATTGCGGAAAAGCGACAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGcc                                                                                                                                                       >  1:2908867/1‑150 (MQ=255)
  ggCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCGTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTAACCACCAAGCCCCGCATACACCGCCAGCGGGTGGTGAATCGCCTGCTCGACAACAAGGCCAg                                                                                                                                                     >  2:2331045/1‑150 (MQ=255)
   gCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGc                                                                                                                                                    <  1:3902075/150‑1 (MQ=255)
     ccGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCct                                                                                                                                                  >  1:3599175/1‑150 (MQ=255)
     ccGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCct                                                                                                                                                  >  1:5144996/1‑150 (MQ=255)
      cGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCctg                                                                                                                                                 <  1:3784259/150‑1 (MQ=255)
        ccTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCctgct                                                                                                                                               <  2:194575/150‑1 (MQ=255)
         cTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCctgctg                                                                                                                                              <  1:1154063/150‑1 (MQ=255)
         cTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCctgctg                                                                                                                                              <  1:2116180/150‑1 (MQ=255)
          tGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCTGGCctgcttg                                                                                                                                             <  2:3827274/150‑3 (MQ=255)
          tGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCTGGCctgcttg                                                                                                                                             <  1:4093416/150‑3 (MQ=255)
           ggACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTAgc                                                                                                                                            <  2:3069749/150‑1 (MQ=255)
            gACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGgcg                                                                                                                                           <  2:4785150/150‑1 (MQ=255)
            gACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGgcg                                                                                                                                           >  1:1017882/1‑150 (MQ=255)
             aCATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCGTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCATGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGgcgc                                                                                                                                          <  2:4005251/150‑1 (MQ=255)
              cATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCt                                                                                                                                         <  1:3812672/150‑1 (MQ=255)
               aTTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCAGGCTGgcctgc                                                                                                                                        <  1:2900638/150‑5 (MQ=255)
                ttGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCct                                                                                                                                                  >  2:4584277/1‑139 (MQ=255)
                 tGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGGCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGa                                                                                                                                      >  1:5325293/1‑150 (MQ=255)
                 tGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGa                                                                                                                                      <  1:2011866/150‑1 (MQ=255)
                 tGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGa                                                                                                                                      >  1:5325615/1‑150 (MQ=255)
                    ggAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATaa                                                                                                                                   >  1:1046255/1‑150 (MQ=255)
                       caaGCGGCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGCGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCt                                                                                                                                <  2:320965/149‑1 (MQ=255)
                        aaGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTg                                                                                                                               <  2:2794158/150‑1 (MQ=255)
                         aGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTgg                                                                                                                              <  1:4902289/150‑1 (MQ=255)
                         aGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTgg                                                                                                                              <  2:4860464/150‑1 (MQ=255)
                             tCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCAt                                                                                                                          >  2:1793636/1‑150 (MQ=255)
                                gCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGcc                                                                                                                       >  2:4510574/1‑150 (MQ=255)
                                 ccccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGGGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCg                                                                                                                      >  1:2461948/1‑150 (MQ=255)
                                  cccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCgg                                                                                                                     >  1:1918050/1‑150 (MQ=255)
                                  cccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCgg                                                                                                                     >  1:4457092/1‑150 (MQ=255)
                                   ccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGc                                                                                                                    >  2:4386540/1‑150 (MQ=255)
                                      ggACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGc                                                                                                                 <  1:541907/150‑1 (MQ=255)
                                      ggACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGc                                                                                                                 <  2:4220194/150‑1 (MQ=255)
                                      ggACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGc                                                                                                                 <  1:3241949/150‑1 (MQ=255)
                                       gACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCg                                                                                                                <  1:4188219/150‑1 (MQ=255)
                                               cATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCa                                                                                                        >  2:3018955/1‑150 (MQ=255)
                                                aTCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCTCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTTCAg                                                                                                       <  2:3865512/150‑1 (MQ=255)
                                                aTCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAg                                                                                                       >  2:3977324/1‑150 (MQ=255)
                                                 tCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAgg                                                                                                      <  1:2012749/150‑1 (MQ=255)
                                                      ttGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTg                                                                                                 >  1:3455426/1‑150 (MQ=255)
                                                      ttGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTg                                                                                                 >  1:2363084/1‑150 (MQ=255)
                                                       tGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTgc                                                                                                <  1:995378/150‑1 (MQ=255)
                                                          gcgcgGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCa                                                                                             <  1:4005004/150‑1 (MQ=255)
                                                               gcggGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTCTTGCAGGGCCTGCGCAGTTGc                                                                                        >  1:998732/1‑150 (MQ=255)
                                                                  gggCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTg                                                                                     <  2:359794/150‑1 (MQ=255)
                                                                      gcgAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGgcgc                                                                                 >  2:3499479/1‑150 (MQ=255)
                                                                      gcgAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGgcgc                                                                                 >  1:1783426/1‑150 (MQ=255)
                                                                               aCGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACg                                                                        <  2:1168908/150‑1 (MQ=255)
                                                                                 gTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGc                                                                      <  1:4646345/150‑1 (MQ=255)
                                                                                   aCGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCAt                                                                    <  1:614922/150‑1 (MQ=255)
                                                                                    cGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCAtt                                                                   <  1:716566/150‑1 (MQ=255)
                                                                                    cGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCAtt                                                                   <  1:3018955/150‑1 (MQ=255)
                                                                                       gTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATtgtg                                                                <  1:288286/150‑1 (MQ=255)
                                                                                       gTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATtgtg                                                                <  2:1716839/150‑1 (MQ=255)
                                                                                       gTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATtgtg                                                                <  1:1195314/150‑1 (MQ=255)
                                                                                         ctcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCt                                                              <  2:4297725/150‑1 (MQ=255)
                                                                                         ctcGACCGTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTAGTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCt                                                              <  2:18381/150‑1 (MQ=255)
                                                                                          tcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTg                                                             >  1:4335919/1‑150 (MQ=255)
                                                                                               cTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCAc                                                        >  1:5130320/1‑150 (MQ=255)
                                                                                               cTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCAc                                                        <  1:1571459/150‑1 (MQ=255)
                                                                                               cTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCAc                                                        >  1:1445111/1‑150 (MQ=255)
                                                                                                 tGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCg                                                      <  2:1472930/150‑1 (MQ=255)
                                                                                                   cccGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCggg                                                    <  1:2974844/150‑1 (MQ=255)
                                                                                                     cGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGcc                                                  >  1:1548304/1‑150 (MQ=255)
                                                                                                        cATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGAc                                               <  1:2822286/150‑1 (MQ=255)
                                                                                                        cATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGAc                                               <  1:5017183/150‑1 (MQ=255)
                                                                                                           ccACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCgg                                            >  2:574569/1‑150 (MQ=255)
                                                                                                                  gAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGc                                     <  2:2908867/150‑1 (MQ=255)
                                                                                                                  gAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGc                                     <  2:4089/150‑1 (MQ=255)
                                                                                                                    gCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCaa                                   <  1:2981389/150‑1 (MQ=255)
                                                                                                                    gCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCaa                                   <  1:2253557/150‑1 (MQ=255)
                                                                                                                        gtggtgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGAcc                               >  2:4163085/1‑150 (MQ=255)
                                                                                                                           gtgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACg                            >  2:4684753/1‑150 (MQ=255)
                                                                                                                             gATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGAc                          >  1:1874056/1‑150 (MQ=255)
                                                                                                                                gcgcGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCtt                       >  1:2716755/1‑150 (MQ=255)
                                                                                                                                 cgcgTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTg                      <  1:4555896/150‑1 (MQ=255)
                                                                                                                                   cgTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTa                    <  2:1610206/150‑1 (MQ=255)
                                                                                                                                        tCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCa               >  1:3922676/1‑150 (MQ=255)
                                                                                                                                          gACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGc             >  2:3785034/1‑150 (MQ=255)
                                                                                                                                           gccaGCAGGCTGGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGcc            >  2:2596645/4‑150 (MQ=255)
                                                                                                                                           gccaGCAGGCTGGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAACGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGcc            >  2:617733/4‑150 (MQ=255)
                                                                                                                                            ccaGCAGGCTGGCCTGCTTGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCt           >  1:3144509/3‑150 (MQ=255)
                                                                                                                                            cAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCt           <  1:2598901/150‑1 (MQ=255)
                                                                                                                                             caGCAGGCTGGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTg          >  2:3148813/2‑150 (MQ=255)
                                                                                                                                             caGCAGGCTGGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTg          >  1:5018012/2‑150 (MQ=255)
                                                                                                                                             caGCAGGCTGGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTg          >  1:340978/2‑150 (MQ=255)
                                                                                                                                              aGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGa         <  2:5421502/150‑1 (MQ=255)
                                                                                                                                              aGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGa         >  1:2752285/1‑150 (MQ=255)
                                                                                                                                               gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGTGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGaa        >  2:5029477/1‑150 (MQ=255)
                                                                                                                                               gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCa               >  1:2505540/1‑143 (MQ=255)
                                                                                                                                               gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCa               <  2:2505540/143‑1 (MQ=255)
                                                                                                                                               gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGaa        >  1:3700390/1‑150 (MQ=255)
                                                                                                                                               gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGaa        <  2:2064663/150‑1 (MQ=255)
                                                                                                                                               gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGaa        >  1:2807174/1‑150 (MQ=255)
                                                                                                                                               gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGaa        >  1:5313445/1‑150 (MQ=255)
                                                                                                                                                   gCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGc    >  2:412185/1‑150 (MQ=255)
                                                                                                                                                   gCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGc    >  1:561678/1‑150 (MQ=255)
                                                                                                                                                     cAGCCTGCTGGCGCTGGATACCCTGGGCATGCCGGCAGCGGTATTGCGGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGCCt  <  2:2205933/150‑1 (MQ=255)
                                                                                                                                                     cAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGCCt  >  1:1655585/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGGCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGCCT  >  NC_002947/1198189‑1198487

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