Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,585,597 T→G 6.3% K942Q (AAG→CAG) ‡ treSB ← bifunctional trehalose synthase B/maltokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,585,5970TG 6.3% 1174.1 / 46.4 347K942Q (AAG→CAG) ‡treSBbifunctional trehalose synthase B/maltokinase
Reads supporting (aligned to +/- strand):  ref base T (161/162);  new base G (10/12);  total (173/174)
Fisher's exact test for biased strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 363 total aligned reads displayed.

TAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGC  >  NC_002947/4585449‑4585746
                                                                                                                                                    |                                                                                                                                                     
 aaGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtgcttgc                                                                                                                                                     <  1:2508236/150‑1 (MQ=255)
    ggCTGTGTTTTGCCCGGCGCTCTTCATGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCTTCACCCTGCACCAAGAGCACCTGGCACAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGTGCCTGCTTGCTCa                                                                                                                                                  <  1:966496/150‑1 (MQ=255)
    ggCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCa                                                                                                                                                  >  2:4305079/1‑150 (MQ=255)
     gCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAg                                                                                                                                                 <  2:456638/150‑1 (MQ=255)
      cTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAgg                                                                                                                                                <  1:1162466/150‑1 (MQ=255)
      cTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAgg                                                                                                                                                <  1:6423718/150‑1 (MQ=255)
      cTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAgg                                                                                                                                                <  1:3501322/150‑1 (MQ=255)
      cTGTGTTTGGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAgg                                                                                                                                                <  2:1537181/150‑1 (MQ=255)
       tgtgTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAggg                                                                                                                                               >  1:5575145/1‑150 (MQ=255)
           ttttGCCCGGCGCTCTTCAAGGGTCCCGGCGGGTTCAACTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGAGCAGCTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCAGCGAGCTCAGGGTGTc                                                                                                                                           >  2:5223869/1‑150 (MQ=255)
            tttGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCg                                                                                                                                          >  1:5965306/1‑150 (MQ=255)
              tGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAt                                                                                                                                        <  2:8244943/150‑1 (MQ=255)
                cccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAtgt                                                                                                                                      >  1:6827146/1‑150 (MQ=255)
                cccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAtgt                                                                                                                                      >  1:4402064/1‑150 (MQ=255)
                 ccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGTACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAtgtg                                                                                                                                     >  2:3892765/1‑150 (MQ=255)
                 ccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAtgtg                                                                                                                                     >  2:3152857/1‑150 (MQ=255)
                 ccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAtgtg                                                                                                                                     >  2:6308157/1‑150 (MQ=255)
                     cgcTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTgg                                                                                                                                 <  1:1007980/150‑1 (MQ=255)
                     cgcTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTgg                                                                                                                                 <  1:7094363/150‑1 (MQ=255)
                     cgcTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTgg                                                                                                                                 >  2:3351356/1‑150 (MQ=255)
                     cgcTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTgg                                                                                                                                 <  1:4960595/150‑1 (MQ=255)
                      gcTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTggg                                                                                                                                >  2:6630052/1‑150 (MQ=255)
                        tcttctAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGcc                                                                                                                              >  2:771981/1‑150 (MQ=255)
                        tctTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGcc                                                                                                                              >  2:268301/1‑150 (MQ=255)
                        tctTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGcc                                                                                                                              <  1:5356876/150‑1 (MQ=255)
                         ctTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCa                                                                                                                             >  1:48483/1‑150 (MQ=255)
                          ttCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAg                                                                                                                            <  1:3656351/150‑1 (MQ=255)
                          ttCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAg                                                                                                                            <  2:1424003/150‑1 (MQ=255)
                          ttCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAg                                                                                                                            >  2:238948/1‑150 (MQ=255)
                          ttCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAg                                                                                                                            >  2:8494783/1‑150 (MQ=255)
                          cacAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGCTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCACGCCCCCCTGAGCCTGCTTTCTCAGGGTGTCGATGTGCTGGGCCAg                                                                                                                            <  1:1749562/148‑1 (MQ=255)
                           tCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAgg                                                                                                                           >  1:731982/1‑150 (MQ=255)
                             aaGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGTCAGGcc                                                                                                                         >  1:3121029/1‑150 (MQ=255)
                             aaGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGcc                                                                                                                         >  2:1058635/1‑150 (MQ=255)
                             aaGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGcc                                                                                                                         <  1:4190598/150‑1 (MQ=255)
                             aaGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGcc                                                                                                                         <  2:3805430/150‑1 (MQ=255)
                               ggggCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGt                                                                                                                       <  1:5326731/150‑1 (MQ=255)
                               ggggCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGGGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGt                                                                                                                       <  2:7036077/150‑1 (MQ=255)
                                 ggCCGGGCGGGTTCACCTTCAAAATCGCTGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCa                                                                                                                     <  2:6077221/150‑1 (MQ=255)
                                 ggCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACTAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCa                                                                                                                     <  2:7919129/150‑1 (MQ=255)
                                    cgggcgggTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGc                                                                                                                  >  1:3382162/1‑150 (MQ=255)
                                       gcgggTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGt                                                                                                               <  1:4504147/150‑1 (MQ=255)
                                       gcgggTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGt                                                                                                               <  2:1497874/150‑1 (MQ=255)
                                        cgggTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGtt                                                                                                              <  1:3724499/150‑1 (MQ=255)
                                             tCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTtgctgc                                                                                                         <  1:8728036/150‑1 (MQ=255)
                                             tCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTtgctgc                                                                                                         <  1:6310765/150‑1 (MQ=255)
                                             tCACCTTCAAAATCGATGAGGTAGGCATCACCCGGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTtgctgc                                                                                                         <  1:1761466/150‑1 (MQ=255)
                                               aCCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAg                                                                                                       <  2:2793167/150‑1 (MQ=255)
                                                   tCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgtcg                                                                                                   >  1:7848711/1‑150 (MQ=255)
                                                     aaaaTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgtcgtc                                                                                                 <  2:1984688/150‑1 (MQ=255)
                                                      aaaTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGtcgtcg                                                                                                <  2:1696162/150‑1 (MQ=255)
                                                       aaTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGa                                                                                               >  1:5154098/1‑150 (MQ=255)
                                                         tCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCCGGTCGTCGAcc                                                                                             >  2:7026970/1‑150 (MQ=255)
                                                         tCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGAcc                                                                                             <  1:3287440/150‑1 (MQ=255)
                                                           gATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCaa                                                                                           >  1:8377101/1‑150 (MQ=255)
                                                           gATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCaa                                                                                           <  2:4750421/150‑1 (MQ=255)
                                                            aTGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAg                                                                                          <  2:3715404/150‑1 (MQ=255)
                                                            aTGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAg                                                                                          <  1:7575517/150‑1 (MQ=255)
                                                            aTGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAg                                                                                          <  2:3508233/150‑1 (MQ=255)
                                                                ggTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACt                                                                                      <  1:591841/150‑1 (MQ=255)
                                                                  tAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTgg                                                                                    <  1:6685508/150‑1 (MQ=255)
                                                                     gCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTa                                                                                 >  2:5578169/1‑150 (MQ=255)
                                                                     gCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTa                                                                                 >  1:6968789/1‑150 (MQ=255)
                                                                       aTCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACACATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATc                                                                               >  2:7807404/1‑150 (MQ=255)
                                                                          aCCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGTCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCt                                                                            >  2:2899312/1‑150 (MQ=255)
                                                                            ccTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGt                                                                          <  2:4211130/150‑1 (MQ=255)
                                                                             cTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTc                                                                         <  1:5673210/150‑1 (MQ=255)
                                                                                  ccacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgtg                                                                    >  1:7418414/1‑150 (MQ=255)
                                                                                   cacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgtgt                                                                   <  2:8378123/150‑1 (MQ=255)
                                                                                          aCCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGa                                                            >  1:7255947/1‑150 (MQ=255)
                                                                                           ccTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGAt                                                           >  1:1414623/1‑150 (MQ=255)
                                                                                               gCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTg                                                       <  2:5997193/150‑1 (MQ=255)
                                                                                                cccAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTgg                                                      <  2:2265173/150‑1 (MQ=255)
                                                                                                 ccAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTggg                                                     <  1:4032330/150‑1 (MQ=255)
                                                                                                   aGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGcc                                                   >  1:8284069/1‑150 (MQ=255)
                                                                                                   aGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGcc                                                   >  2:7390106/1‑150 (MQ=255)
                                                                                                       gCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCcagca                                               <  2:8782070/150‑1 (MQ=255)
                                                                                                        cAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCcagcag                                              >  2:1325227/1‑150 (MQ=255)
                                                                                                         aGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAgg                                             >  2:56897/1‑150 (MQ=255)
                                                                                                           gTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTc                                           >  1:4592898/1‑150 (MQ=255)
                                                                                                           gTCGCCGTGTACACGCATAAGAAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTc                                           >  1:205966/1‑150 (MQ=255)
                                                                                                             cGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGa                                         <  1:5509775/150‑1 (MQ=255)
                                                                                                              gCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAg                                        >  1:6513656/1‑150 (MQ=255)
                                                                                                                cGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGCCGAgcg                                      <  2:4934618/150‑3 (MQ=255)
                                                                                                                     acacGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTgg                                 >  1:2746910/1‑150 (MQ=255)
                                                                                                                       acGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGt                               >  2:3031469/1‑150 (MQ=255)
                                                                                                                       acGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGt                               >  1:396234/1‑150 (MQ=255)
                                                                                                                                atcgacACCCTGAGCAAGCAGGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcagttcagt                      >  2:3982306/6‑150 (MQ=255)
                                                                                                                                  ggCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGc                    >  2:1188841/1‑150 (MQ=255)
                                                                                                                                  cgacACCCTGAGCAAGCAGGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGc                    >  2:3584713/4‑150 (MQ=255)
                                                                                                                                   gCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCt                   <  2:5354836/150‑1 (MQ=255)
                                                                                                                                   gCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCt                   >  2:7357523/1‑150 (MQ=255)
                                                                                                                                    acACCCTGAGCAAGCAGGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTg                  >  1:1561087/2‑150 (MQ=255)
                                                                                                                                          tGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCtgt            <  1:5818474/150‑1 (MQ=255)
                                                                                                                                          tGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCtgt            >  1:2849757/1‑150 (MQ=255)
                                                                                                                                           gAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCtgtg           <  1:1074370/150‑1 (MQ=255)
                                                                                                                                                 gcttgctCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACAcc     >  1:8786515/1‑150 (MQ=255)
                                                                                                                                                  cttgctCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCa    <  1:3892765/150‑1 (MQ=255)
                                                                                                                                                  cttgctCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCa    >  1:6244480/1‑150 (MQ=255)
                                                                                                                                                    tgctCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAgc  >  2:2201410/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACACCAGC  >  NC_002947/4585449‑4585746

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: