Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,375,260 G→T 100% D387Y (GAT→TAT)  degQ → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,375,2600GT99.5% 623.0 / ‑4.1 192D387Y (GAT→TAT) degQserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (88/103);  total (88/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 195 total aligned reads displayed.

CACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCA  >  CP009273/3375115‑3375403
                                                                                                                                                 |                                                                                                                                               
cACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAatat                                                                                                                                               <  2:1941557/148‑1 (MQ=255)
 aCGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACGAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCATCTAAAATATg                                                                                                                                              <  1:2508641/148‑1 (MQ=255)
 aCGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATg                                                                                                                                              <  1:468274/148‑1 (MQ=255)
   gAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATggc                                                                                                                                            >  1:1175255/1‑148 (MQ=255)
       gTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCa                                                                                                                                        <  1:1007978/148‑1 (MQ=255)
         gAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCaaa                                                                                                                                      <  2:3759607/148‑1 (MQ=255)
           aGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAgg                                                                                                                                    <  2:2322218/148‑1 (MQ=255)
            gCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGt                                                                                                                                   >  2:1205687/1‑148 (MQ=255)
               tGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTAtt                                                                                                                                >  1:260994/1‑148 (MQ=255)
               tGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTAtt                                                                                                                                <  1:3530068/148‑1 (MQ=255)
                 gCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTaa                                                                                                                              >  2:1290465/1‑148 (MQ=255)
                    tgctgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAAt                                                                                                                           <  2:2653581/148‑1 (MQ=255)
                       tgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGa                                                                                                                        <  1:3633261/148‑1 (MQ=255)
                       tgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGa                                                                                                                        <  2:2618049/148‑1 (MQ=255)
                       tgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGa                                                                                                                        <  2:866814/148‑1 (MQ=255)
                            aaCGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAg                                                                                                                   <  2:196441/148‑1 (MQ=255)
                             aCGGCAAACCACTGGAAGTAGAAGTGACGCTCGATAGCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACNCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGt                                                                                                                  >  1:1933310/1‑148 (MQ=255)
                             aCGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGt                                                                                                                  >  2:2123828/1‑148 (MQ=255)
                              cGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGtt                                                                                                                 >  2:3508213/1‑148 (MQ=255)
                                   aaCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCaa                                                                                                            <  2:2783355/148‑1 (MQ=255)
                                    aCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCACAGGTATTAAAATCGATGAAGTTGTCaaa                                                                                                           >  1:701617/1‑148 (MQ=255)
                                    aCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCaaa                                                                                                           >  1:520037/1‑148 (MQ=255)
                                    aCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCaaa                                                                                                           >  1:1666409/1‑148 (MQ=255)
                                      cACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAgg                                                                                                         <  2:1926156/148‑1 (MQ=255)
                                          ggAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGc                                                                                                     <  2:3029618/148‑1 (MQ=255)
                                              gTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAg                                                                                                 <  1:3230626/148‑1 (MQ=255)
                                              gTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAg                                                                                                 <  1:1413638/148‑1 (MQ=255)
                                               tAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgc                                                                                                >  2:2901137/1‑148 (MQ=255)
                                                 gAAGTGACGCTCGATACCAGCACCTCTACGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgctg                                                                                              <  2:2165512/148‑1 (MQ=255)
                                                  aaGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgctgc                                                                                             <  1:22471/148‑1 (MQ=255)
                                                    gTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTc                                                                                           <  2:1315762/148‑1 (MQ=255)
                                                       aCGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAgg                                                                                        >  2:2797741/1‑148 (MQ=255)
                                                       aCGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAgg                                                                                        <  2:804436/148‑1 (MQ=255)
                                                           tCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctgg                                                                                    <  2:1634041/148‑1 (MQ=255)
                                                             gATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctggct                                                                                  <  2:53028/148‑1 (MQ=255)
                                                              aTACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCtt                                                                                 >  1:1527001/1‑148 (MQ=255)
                                                              aTACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCtt                                                                                 >  1:3427092/1‑148 (MQ=255)
                                                                   aGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaaa                                                                            >  1:2315692/1‑148 (MQ=255)
                                                                    gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaaaa                                                                           <  1:3230217/148‑1 (MQ=255)
                                                                    gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaaaa                                                                           >  1:2653020/1‑148 (MQ=255)
                                                                    gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaaaa                                                                           <  1:347004/148‑1 (MQ=255)
                                                                     cACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaaaaa                                                                          >  1:1463881/1‑148 (MQ=255)
                                                                      aCCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAg                                                                         >  1:509010/1‑148 (MQ=255)
                                                                       ccTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGa                                                                        >  2:573373/1‑148 (MQ=255)
                                                                        ctattCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGAc                                                                       <  2:1810131/145‑1 (MQ=255)
                                                                         tctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACg                                                                      <  2:1241818/148‑1 (MQ=255)
                                                                            tcgtcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGAtg                                                                   <  2:1016622/148‑1 (MQ=255)
                                                                             cgtcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGAtgt                                                                  <  1:3392249/148‑1 (MQ=255)
                                                                               tcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGa                                                                <  1:79031/148‑1 (MQ=255)
                                                                                 ggCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAATGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATc                                                              <  2:2144096/148‑1 (MQ=255)
                                                                                   ccAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAt                                                            <  1:276834/148‑1 (MQ=255)
                                                                                   ccAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAt                                                            >  2:3800746/1‑148 (MQ=255)
                                                                                    cAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAtt                                                           >  1:3138903/1‑148 (MQ=255)
                                                                                    cAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAtt                                                           <  2:260994/148‑1 (MQ=255)
                                                                                     aGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTg                                                          <  2:2861736/148‑1 (MQ=255)
                                                                                       cgcTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATGAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGc                                                        >  1:3507106/1‑148 (MQ=255)
                                                                                          tGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTc                                                     >  1:2692000/1‑148 (MQ=255)
                                                                                          tGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTc                                                     >  1:222685/1‑148 (MQ=255)
                                                                                           gAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCa                                                    <  1:537499/148‑1 (MQ=255)
                                                                                           gAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCa                                                    >  1:2952407/1‑148 (MQ=255)
                                                                                           gAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCa                                                    >  1:2952832/1‑148 (MQ=255)
                                                                                           gAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCa                                                    >  1:3300796/1‑148 (MQ=255)
                                                                                           gAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCa                                                    >  1:1885719/1‑148 (MQ=255)
                                                                                            aaaTGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCaa                                                   >  2:1405391/1‑148 (MQ=255)
                                                                                            aaaTGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCAAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCaa                                                   >  2:1681164/1‑148 (MQ=255)
                                                                                              aTGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAAcc                                                 >  1:1043624/1‑148 (MQ=255)
                                                                                                 aTCACGCCAGCTCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACcgcg                                              <  1:2731757/148‑1 (MQ=255)
                                                                                                  tCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGa                                             <  2:1067246/148‑1 (MQ=255)
                                                                                                  tCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGa                                             <  2:3578813/148‑1 (MQ=255)
                                                                                                     cGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg                                          <  1:2923286/148‑1 (MQ=255)
                                                                                                     cGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg                                          <  1:37207/148‑1 (MQ=255)
                                                                                                       ccAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCggg                                        <  2:1962314/148‑1 (MQ=255)
                                                                                                       ccAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCggg                                        <  2:1773097/148‑1 (MQ=255)
                                                                                                       ccAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCggg                                        <  2:1480497/148‑1 (MQ=255)
                                                                                                        cAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGt                                       >  2:3732099/1‑148 (MQ=255)
                                                                                                               ggAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCg                                <  1:1767823/148‑1 (MQ=255)
                                                                                                                gAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGa                               <  1:2711440/148‑1 (MQ=255)
                                                                                                                gAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGa                               >  1:3576744/1‑148 (MQ=255)
                                                                                                                gAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGa                               >  1:1264579/1‑148 (MQ=255)
                                                                                                                gAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGa                               <  2:1322735/148‑1 (MQ=255)
                                                                                                                 aaGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGAt                              <  2:2315692/148‑1 (MQ=255)
                                                                                                                   ggTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTg                            <  1:1205742/148‑1 (MQ=255)
                                                                                                                   ggTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTg                            >  1:28172/1‑148 (MQ=255)
                                                                                                                   ggTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTg                            >  1:3098174/1‑148 (MQ=255)
                                                                                                                     tGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCt                          <  2:101322/148‑1 (MQ=255)
                                                                                                                        aaCGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGaa                       <  1:1546587/148‑1 (MQ=255)
                                                                                                                         aCGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGaaa                      <  1:3667167/148‑1 (MQ=255)
                                                                                                                         aCGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGaaa                      <  2:765779/148‑1 (MQ=255)
                                                                                                                            ttGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGc                   >  2:3428728/1‑148 (MQ=255)
                                                                                                                             tGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCg                  >  1:1343551/1‑148 (MQ=255)
                                                                                                                                gCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaa               <  1:402350/148‑1 (MQ=255)
                                                                                                                                   aTGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGt            >  1:2115525/1‑148 (MQ=255)
                                                                                                                                    tGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTg           <  1:3771285/148‑1 (MQ=255)
                                                                                                                                     ggTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGc          >  2:1901427/1‑148 (MQ=255)
                                                                                                                                      gTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCt         >  1:3346096/1‑148 (MQ=255)
                                                                                                                                      gTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCt         <  1:2375215/148‑1 (MQ=255)
                                                                                                                                         aGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggc      >  1:589768/1‑148 (MQ=255)
                                                                                                                                           cTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggcgg    >  2:210717/1‑148 (MQ=255)
                                                                                                                                             aaaaTATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCa  >  2:2583506/1‑148 (MQ=255)
                                                                                                                                             aaaaTATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCa  >  2:1808944/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCA  >  CP009273/3375115‑3375403

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: