Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 710,187 A→T 90.3% S392C (AGT→TGT)  pgm → phosphoglucomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273710,1870AT90.3% 435.5 / 31.9 144S392C (AGT→TGT) pgmphosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base A (8/6);  new base T (60/70);  total (68/76)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 150 total aligned reads displayed.

CCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTT  >  CP009273/710044‑710333
                                                                                                                                               |                                                                                                                                                  
ccGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGc                                                                                                                                                <  2:2853129/148‑1 (MQ=255)
  gTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAg                                                                                                                                              <  2:2476237/148‑1 (MQ=255)
  gTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAg                                                                                                                                              >  2:267504/1‑148 (MQ=255)
    cGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAggg                                                                                                                                            >  2:1477443/1‑148 (MQ=255)
    cGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAggg                                                                                                                                            >  2:1139763/1‑148 (MQ=255)
      gTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCATGGGc                                                                                                                                          <  2:1873488/148‑1 (MQ=255)
         aaaCGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGcttc                                                                                                                                       <  1:1278008/148‑1 (MQ=255)
           aCGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGcttcct                                                                                                                                     >  1:1366515/1‑148 (MQ=255)
           aCGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGcttcct                                                                                                                                     >  1:3013714/1‑148 (MQ=255)
             gCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGcttccttc                                                                                                                                   <  1:1082113/148‑1 (MQ=255)
             gCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGcttccttc                                                                                                                                   <  1:755866/148‑1 (MQ=255)
                 gTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGc                                                                                                                               >  2:1061882/1‑148 (MQ=255)
                   ttCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGt                                                                                                                             <  1:1383530/148‑1 (MQ=255)
                   ttCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGt                                                                                                                             <  2:1867244/148‑1 (MQ=255)
                    tCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGGAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAGGGGCTTCCTTCCTGCGtt                                                                                                                            >  1:1709123/1‑148 (MQ=255)
                    tCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGtt                                                                                                                            >  1:2066080/1‑148 (MQ=255)
                    tCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGtt                                                                                                                            <  2:405265/148‑1 (MQ=255)
                        cTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGa                                                                                                                        >  2:2043735/1‑148 (MQ=255)
                        cTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGa                                                                                                                        >  1:486345/1‑148 (MQ=255)
                            gatgatCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGc                                                                                                                    <  1:2628612/148‑1 (MQ=255)
                             atgatCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCa                                                                                                                   >  1:2457452/1‑148 (MQ=255)
                                   gACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGt                                                                                                             >  1:405525/1‑148 (MQ=255)
                                    aCCATGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTg                                                                                                            <  2:154855/148‑1 (MQ=255)
                                     ccGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTgg                                                                                                           >  1:450886/1‑148 (MQ=255)
                                         gtgGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCa                                                                                                       >  2:3267047/1‑148 (MQ=255)
                                             tCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGa                                                                                                   <  1:2553514/148‑1 (MQ=255)
                                             tCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGa                                                                                                   <  1:2768538/148‑1 (MQ=255)
                                                aCGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACaa                                                                                                <  1:2488536/148‑1 (MQ=255)
                                                aCGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACaa                                                                                                <  1:2488562/148‑1 (MQ=255)
                                                aCGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACaa                                                                                                <  2:535374/148‑1 (MQ=255)
                                                 cGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACaaa                                                                                               >  2:1004583/1‑148 (MQ=255)
                                                 cGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTACTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACaaa                                                                                               >  2:1929227/1‑148 (MQ=255)
                                                 cGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACaaa                                                                                               >  2:1690113/1‑148 (MQ=255)
                                                    cTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGAc                                                                                            <  2:1572639/148‑1 (MQ=255)
                                                      tGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACgg                                                                                          >  1:2810092/1‑148 (MQ=255)
                                                         gCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcat                                                                                       <  2:2530724/148‑1 (MQ=255)
                                                         gCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcat                                                                                       <  1:1566554/148‑1 (MQ=255)
                                                          ccGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcatc                                                                                      <  1:3024606/148‑1 (MQ=255)
                                                          ccGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcatc                                                                                      <  1:1277437/148‑1 (MQ=255)
                                                             tAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcatcatc                                                                                   <  1:320994/148‑1 (MQ=255)
                                                             tAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGcatcatc                                                                                   <  1:3013292/148‑1 (MQ=255)
                                                                 cTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCAtgt                                                                               <  1:664151/148‑1 (MQ=255)
                                                                   ggTAGAAGTCCCGGTAGGTTTCACATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATgtgt                                                                             <  2:273003/148‑1 (MQ=255)
                                                                   ggTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATgtgt                                                                             <  2:805982/148‑1 (MQ=255)
                                                                    gTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGCTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTc                                                                            <  2:957750/148‑1 (MQ=255)
                                                                      aGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTctg                                                                          >  1:197104/1‑148 (MQ=255)
                                                                      aGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTctg                                                                          >  2:570885/1‑148 (MQ=255)
                                                                      aGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGCCCACCGACAAAGACGGCATCATCATGTGTctg                                                                          >  1:197146/1‑148 (MQ=255)
                                                                       gAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTctgc                                                                         >  2:3037686/1‑148 (MQ=255)
                                                                       gAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTctgc                                                                         >  2:3037726/1‑148 (MQ=255)
                                                                       gAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTctgc                                                                         >  2:1583098/1‑148 (MQ=255)
                                                                          gTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTgg                                                                      >  2:2383637/1‑148 (MQ=255)
                                                                           tCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTggc                                                                     <  1:3293698/148‑1 (MQ=255)
                                                                           tCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTAGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTggc                                                                     >  2:820636/1‑148 (MQ=255)
                                                                           tCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTAGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTggc                                                                     >  2:820644/1‑148 (MQ=255)
                                                                           tCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGACAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTggc                                                                     >  2:171089/1‑148 (MQ=255)
                                                                                 tAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGa                                                               >  2:237078/1‑148 (MQ=255)
                                                                                   ggTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGaaa                                                             >  1:363397/1‑148 (MQ=255)
                                                                                     tttCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATc                                                           >  2:22025/1‑148 (MQ=255)
                                                                                      ttCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCa                                                          >  2:2466648/1‑148 (MQ=255)
                                                                                      ttCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCa                                                          <  2:2162123/148‑1 (MQ=255)
                                                                                      ttCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCa                                                          >  2:931449/1‑148 (MQ=255)
                                                                                      ttCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCa                                                          >  2:932755/1‑148 (MQ=255)
                                                                                       tCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCAc                                                         <  1:2369010/148‑1 (MQ=255)
                                                                                         aaaTGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCg                                                       <  2:950093/148‑1 (MQ=255)
                                                                                         aaaTGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCg                                                       <  2:2241194/148‑1 (MQ=255)
                                                                                         aaaTGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCg                                                       <  2:3181347/148‑1 (MQ=255)
                                                                                           aTGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCt                                                     >  2:1652086/1‑148 (MQ=255)
                                                                                            tGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTg                                                    <  2:1310862/148‑1 (MQ=255)
                                                                                              gTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGAGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTc                                                  <  1:1540292/148‑1 (MQ=255)
                                                                                               tttGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCa                                                 >  2:1238083/1‑148 (MQ=255)
                                                                                                  gTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGTGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCg                                              <  2:2843739/148‑1 (MQ=255)
                                                                                                    cGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGt                                            <  1:1202910/148‑1 (MQ=255)
                                                                                                     gATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTa                                           <  1:2231466/148‑1 (MQ=255)
                                                                                                      aTGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTaa                                          >  1:2963872/1‑148 (MQ=255)
                                                                                                      aTGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTaa                                          >  1:2740164/1‑148 (MQ=255)
                                                                                                      aTGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTaa                                          >  2:2548957/1‑148 (MQ=255)
                                                                                                        ggTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCCACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGa                                        >  2:410720/1‑148 (MQ=255)
                                                                                                           cTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAAcc                                     >  1:2199049/1‑148 (MQ=255)
                                                                                                             gTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCg                                   >  2:1073078/1‑148 (MQ=255)
                                                                                                              ttCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGc                                  >  2:674305/1‑148 (MQ=255)
                                                                                                              ttCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGc                                  >  2:674267/1‑148 (MQ=255)
                                                                                                              ttCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGAGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGc                                  <  2:3099225/148‑1 (MQ=255)
                                                                                                                cGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAg                                >  2:1171520/1‑148 (MQ=255)
                                                                                                                   cGGCAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGaa                             >  1:2296742/1‑148 (MQ=255)
                                                                                                                     gCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAaca                           <  1:2818850/148‑1 (MQ=255)
                                                                                                                      cAGCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAacac                          <  1:2043735/148‑1 (MQ=255)
                                                                                                                        gCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTa                        >  1:749845/1‑148 (MQ=255)
                                                                                                                        gCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTa                        <  2:1366515/148‑1 (MQ=255)
                                                                                                                        gCTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTa                        >  2:1408559/1‑148 (MQ=255)
                                                                                                                         cTTCGGCTTTGGCGGCGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTAc                       >  2:366347/1‑148 (MQ=255)
                                                                                                                                tttGGCGGAGAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAAc                <  1:1732757/148‑1 (MQ=255)
                                                                                                                                       gcgAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCaaa         <  1:671888/148‑1 (MQ=255)
                                                                                                                                       gcgAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCaaa         <  2:3121803/148‑1 (MQ=255)
                                                                                                                                       gcgAAGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCaaa         <  2:3121785/148‑1 (MQ=255)
                                                                                                                                          aaGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACg      <  2:2892390/148‑1 (MQ=255)
                                                                                                                                          aaGAGTGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACg      <  2:2892017/148‑1 (MQ=255)
                                                                                                                                            gagtgTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCt    <  1:2171789/148‑1 (MQ=255)
                                                                                                                                            gagtgTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCt    >  2:2201552/1‑148 (MQ=255)
                                                                                                                                              gagTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCttt  <  2:486345/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGTTTGTCGATGGTCTGTTCGACGGCAGCTTCGGCTTTGGCGGCGAAGAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTT  >  CP009273/710044‑710333

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: