Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,940,829 A→T 9.1% intergenic (+7/‑27) metZ → / → metW tRNA‑Met/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,940,8290AT 9.1% 498.0 / 2.3 171intergenic (+7/‑27)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (50/101);  new base T (8/11);  total (59/112)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

CGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGA  >  CP009273/2940683‑2940976
                                                                                                                                                  |                                                                                                                                                    
cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAtt                                                                                                                                                    <  2:2618006/149‑1 (MQ=255)
 gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                                                                                   <  2:3789527/149‑1 (MQ=255)
    aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                                <  2:2062754/149‑1 (MQ=255)
        gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAg                                                                                                                                            <  2:2181424/149‑1 (MQ=255)
          gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTa                                                                                                                                         <  1:1670351/149‑1 (MQ=255)
           tAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaa                                                                                                                                        <  2:3987680/149‑1 (MQ=255)
           tAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTGATGAATTaa                                                                                                                                        >  2:2383078/1‑149 (MQ=255)
             aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAg                                                                                                                                      >  1:1246446/1‑149 (MQ=255)
             aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAg                                                                                                                                      <  1:1028248/149‑1 (MQ=255)
              aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGc                                                                                                                                     >  1:1180637/1‑149 (MQ=255)
              aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaagc                                                                                                                                    <  2:337810/149‑4 (MQ=255)
                aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAg                                                                                                                                   <  2:811982/149‑1 (MQ=255)
                aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAg                                                                                                                                   <  1:3328575/149‑1 (MQ=255)
                aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTGATGAAGTAAAGCAg                                                                                                                                   >  2:3491630/1‑149 (MQ=255)
                   aGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTAc                                                                                                                                <  2:3095981/149‑1 (MQ=255)
                   aGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTAc                                                                                                                                >  2:4536899/1‑149 (MQ=255)
                   aGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTAc                                                                                                                                <  2:5188943/149‑1 (MQ=255)
                     gTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTGATGAAGTAAAGCAGTACgg                                                                                                                              >  1:4163762/1‑149 (MQ=255)
                      ttGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGt                                                                                                                             <  1:2239855/149‑1 (MQ=255)
                        gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGa                                                                                                                           <  2:4734679/149‑1 (MQ=255)
                            gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGAcgcg                                                                                                                       >  1:3823413/1‑149 (MQ=255)
                              aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCggg                                                                                                                     <  2:2862684/149‑1 (MQ=255)
                              aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTAGAGGAATAAAATCAGTACGTTTACGCggg                                                                                                                     >  1:558618/1‑149 (MQ=18)
                               aaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATTTTGATCAAGAAAAGCAGTACGGTGACGCgggg                                                                                                                    >  1:2566603/1‑149 (MQ=255)
                                aCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGt                                                                                                                   <  1:797688/149‑1 (MQ=255)
                                 cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTg                                                                                                                  <  1:4017607/149‑1 (MQ=255)
                                 cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTg                                                                                                                  <  1:1398389/149‑1 (MQ=255)
                                  gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTgg                                                                                                                 <  2:4002876/149‑1 (MQ=255)
                                  gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTgg                                                                                                                 <  1:406893/149‑1 (MQ=255)
                                  gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTgg                                                                                                                 <  1:406076/149‑1 (MQ=255)
                                  gagAGCTGAGTTTAGTGCTCATCCAGGGACGGTGGGTTGAGCATCCTGGTATATCGTCTTGCTATTAACCCGAATGTCGTAGGTTCAAATCCGGCACCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGAGGGTTgg                                                                                                                 <  2:4617511/147‑1 (MQ=255)
                                   cgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGa                                                                                                                <  2:3868858/149‑1 (MQ=255)
                                         gAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGAAAAGCAGTACGGTGACGCGGGGTGGAGCAGcc                                                                                                          >  1:3601531/1‑149 (MQ=255)
                                          aGTATAGTGCGCATCCACGGACGCGGGGTGGAGCATCCTGGTAGCTCGTCTGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTACAGCAGTACGGTGACGCGGGGTGGAGCAGACt                                                                                                         <  1:3937108/149‑1 (MQ=255)
                                           gTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATAGGCGGCAGAAGAGCAGTACGGTGACGCGGGGGGGAGCAGCCtc                                                                                                        >  1:3514529/1‑148 (MQ=17)
                                             ataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAATAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGt                                                                                                      >  1:2372484/1‑149 (MQ=255)
                                             agaGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGt                                                                                                      <  2:1066616/147‑1 (MQ=255)
                                              taGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATTTGAAGGACGAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTa                                                                                                     >  1:214052/1‑149 (MQ=14)
                                               aGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAATAAAAGCAGTACGGTGACGCGGTGTGGAGCAGCCTGGTAg                                                                                                    >  2:1870556/1‑149 (MQ=255)
                                               aGTGCGCAGCCCCGGACGGGGGGTGGATCATGCGGGGATGTCGTCGGGCTAATTACCCGTAGGTCGTCGGTTCAAATCCCGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTa                                                                                                     <  2:490307/149‑1 (MQ=11)
                                                gTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGc                                                                                                   <  2:767564/149‑1 (MQ=255)
                                                 tGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCt                                                                                                  <  2:2865846/149‑1 (MQ=255)
                                                 tGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCt                                                                                                  <  2:1238259/149‑1 (MQ=255)
                                                     cATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgtc                                                                                              >  1:2723437/1‑149 (MQ=255)
                                                     cATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCGCCTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAGGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgtc                                                                                              >  1:1287173/1‑149 (MQ=255)
                                                      aTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCACAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgtcg                                                                                             <  1:4261610/149‑1 (MQ=255)
                                                            ggACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCa                                                                                       <  1:228568/149‑1 (MQ=255)
                                                            ggACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCa                                                                                       <  2:306982/149‑1 (MQ=255)
                                                            ggACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCa                                                                                       <  2:306932/149‑1 (MQ=255)
                                                             gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCAt                                                                                      <  1:4004819/149‑1 (MQ=40)
                                                                gcgGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAc                                                                                   <  2:2860966/149‑1 (MQ=40)
                                                                      tGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAg                                                                             <  2:2301398/149‑1 (MQ=40)
                                                                      tGGAGCAGACTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAg                                                                             <  2:4911788/149‑1 (MQ=255)
                                                                         agcagcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtc                                                                          <  1:2799933/149‑1 (MQ=40)
                                                                         agcagcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtc                                                                          <  1:466688/149‑1 (MQ=40)
                                                                         agcagcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtc                                                                          <  1:2924102/149‑1 (MQ=40)
                                                                               cTGGTAGCTCGTCGGGCTCATAAGCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATTTGGATGAAGTCACGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGt                                                                    <  1:3097926/149‑1 (MQ=11)
                                                                                     gCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTGTGGAATAAAAGCAGTACGGTGACGCGGGGTGGAGCATCCGGGTAGCTCGTCGGGCGCAAATCCAGAGGGTGTTGGGTTctatt                                                              >  2:1649419/1‑145 (MQ=11)
                                                                                       tcgtcgGGCTCATAACCCGAAGGTCGTCGGTTCATATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAACGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATcc                                                            <  1:835806/149‑1 (MQ=255)
                                                                                         gtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAATAAAAGCAGTACGGTGAAGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGACGTCGGTTCAAATCCgg                                                          >  2:4807901/1‑149 (MQ=255)
                                                                                           cgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGcc                                                        <  2:4246891/149‑1 (MQ=38)
                                                                                             ggCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGcccc                                                      <  1:191086/149‑1 (MQ=40)
                                                                                              gCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtcg                                                                         >  1:115089/1‑129 (MQ=255)
                                                                                               cTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGCGCAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGGCGACGGTACACGTCCGGCCCCCg                                                    >  2:2639377/1‑149 (MQ=14)
                                                                                                 cATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGCCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGcc                                                  >  1:4669012/1‑148 (MQ=255)
                                                                                                  aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCaa                                                 >  1:1543836/1‑149 (MQ=40)
                                                                                                   tAACCCGAAGGTCGTCGGTTCTAATCCGGACCCCGCAACCAAATATATTTTGTTGTAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAc                                                <  1:4968707/149‑1 (MQ=11)
                                                                                                    aaCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTGGAGGAACAAAAGCAGTAAGGTGAAGCGGGGTGGAGCAGCCTGGTAGCTCGTCGCGCGCATAAACCGAAGGTCGGCGGTTAAAAACCGCCCCCCCCAAcc                                               >  1:4558453/1‑149 (MQ=11)
                                                                                                     aCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCa                                              >  1:1614259/1‑149 (MQ=40)
                                                                                                      cccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCaa                                             <  1:2782218/149‑1 (MQ=40)
                                                                                                      cccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAATTTTGATGAGGAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCTTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCCACCaa                                             >  1:3462468/1‑149 (MQ=11)
                                                                                                       ccGAAGGTCGTCTGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTCAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAt                                            <  2:1213957/149‑1 (MQ=255)
                                                                                                            ggTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTAGATGAAGAAAAGCATTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGACGTCGGGACAAATCCGGCCCCCGCAACCAAACAAAt                                       >  1:1257337/1‑149 (MQ=14)
                                                                                                              tcgtctGTTCAAATCCGGCCCCCGCAACCAATTACCAATTGTTGAAGTCACCGAGTCCGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                     <  2:4163762/149‑1 (MQ=14)
                                                                                                              tcgtcgGTTCAACTCCGGCCCCCGCAACCACTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                     <  1:2087527/149‑1 (MQ=255)
                                                                                                                  cgGTTCCAATCCGCCCCCAGCAACCAAATACATTTTGATGAAGTCTTGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgatg                                 <  2:671760/149‑1 (MQ=14)
                                                                                                                   ggTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgatga                                <  2:2144999/149‑1 (MQ=255)
                                                                                                                      tCAAATCCGGCCCCCGCAACCAATTACATTTTGAAGAAGTAAAGCAGTCCGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGt                             <  2:2851478/149‑1 (MQ=255)
                                                                                                                      tCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAGATGATGAAGt                             >  1:4388327/1‑149 (MQ=255)
                                                                                                                           tCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGaaactg                        >  2:443857/1‑146 (MQ=255)
                                                                                                                           acaggCACCCGCAACGAATTCAATTTTGATGTAGTAATGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAg                        <  2:1180637/146‑1 (MQ=17)
                                                                                                                              ggCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAg                     <  2:3969341/149‑1 (MQ=255)
                                                                                                                              ggCCCCCGCAACCAATTAAAATGTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGGGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAg                     <  1:3764535/149‑1 (MQ=255)
                                                                                                                                cccccGCAACCAATTACAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTa                   <  2:1504399/149‑1 (MQ=255)
                                                                                                                                 ccccGCAACCAATTCTAATTTGATGAAGTAAGGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTAc                  <  2:328699/149‑1 (MQ=255)
                                                                                                                                 ccccGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTAc                  <  1:1904564/149‑1 (MQ=255)
                                                                                                                                 ccccGCAACCAAATACAAATTGCTGAAGTAATGCAGTCCGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTAc                  <  2:1246539/149‑1 (MQ=18)
                                                                                                                                  cccGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACg                 <  1:1558498/149‑1 (MQ=255)
                                                                                                                                   ccGCAACCAATTATAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACgg                <  2:3386708/149‑1 (MQ=255)
                                                                                                                                   ccGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTAGAGGAAGTAAAAGCAGTACgg                >  2:2200228/1‑149 (MQ=255)
                                                                                                                                     gCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTg              <  1:4275688/149‑1 (MQ=255)
                                                                                                                                        aCCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTCAAGAAGTAAAAGCAGTACGGTGAcg           >  2:2356102/1‑149 (MQ=255)
                                                                                                                                        aCCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAATAAGAGGAAGAAAAAGAACTACAGTCAcc           >  1:2466092/1‑148 (MQ=14)
                                                                                                                                         ccAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAATTTCAACAAGAAAATGCTGTAGGGTGAcgc          >  2:614614/1‑149 (MQ=14)
                                                                                                                                           aaTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgg        <  1:3055843/149‑1 (MQ=255)
                                                                                                                                           aaTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgg        <  2:3916857/149‑1 (MQ=255)
                                                                                                                                            aTTCAATTTTGATGAAGTATCGGATTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCggg       <  2:507269/149‑1 (MQ=17)
                                                                                                                                            aTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGAGGAAGTAATAGCAGTAAGGTGACGAggg       >  2:129808/1‑149 (MQ=255)
                                                                                                                                                aaaTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTgg   >  1:3542475/1‑149 (MQ=255)
                                                                                                                                                 aaTTTGATGAAGTCACGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGa  <  1:3491630/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGA  >  CP009273/2940683‑2940976

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