Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,684,105 C→A 100% intergenic (+97/‑4) manA → / → ydgA mannose‑6‑phosphate isomerase/DUF945 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,684,1050CA98.5% 646.3 / 0.6 201intergenic (+97/‑4)manA/ydgAmannose‑6‑phosphate isomerase/DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (107/90);  total (108/93)
Fisher's exact test for biased strand distribution p-value = 9.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-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 203 total aligned reads displayed.

ACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGA  >  CP009273/1683961‑1684247
                                                                                                                                                |                                                                                                                                              
aCTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATTACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAAtt                                                                                                                                             <  2:1237521/148‑1 (MQ=255)
 cTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTa                                                                                                                                            <  1:323196/148‑1 (MQ=255)
  tGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTAt                                                                                                                                           <  2:1088376/148‑1 (MQ=255)
  tGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTAt                                                                                                                                           <  2:2864305/148‑1 (MQ=255)
   gTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATg                                                                                                                                          <  1:1526044/148‑1 (MQ=255)
       aaGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATa                                                                                                                                      <  1:344748/148‑1 (MQ=255)
        aGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATaa                                                                                                                                     <  1:2565659/148‑1 (MQ=255)
         ggCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATaaa                                                                                                                                    <  2:36273/148‑1 (MQ=255)
             aCGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGc                                                                                                                                <  2:3528271/148‑1 (MQ=255)
                gCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTgg                                                                                                                             <  1:1009726/148‑1 (MQ=255)
                 ccGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGt                                                                                                                            >  2:1583473/1‑148 (MQ=255)
                 ccGTGTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAAATAACATCTGTTGCTAAGCGTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCGGTGTTGTTTTTTTGGCTAATGATTTTGTATAAATCGCTGGt                                                                                                                            >  1:330900/1‑148 (MQ=255)
                    tttAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGc                                                                                                                         >  2:2048981/1‑148 (MQ=255)
                    tttAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGc                                                                                                                         >  2:538867/1‑148 (MQ=255)
                       aGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGt                                                                                                                      <  2:2482345/148‑1 (MQ=255)
                        gcgcGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTa                                                                                                                     <  1:2743437/148‑1 (MQ=255)
                             tgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGt                                                                                                                >  1:2766904/1‑148 (MQ=255)
                              gtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTc                                                                                                               <  1:2805609/148‑1 (MQ=255)
                                ttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAt                                                                                                             >  1:295035/1‑148 (MQ=255)
                                ttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAt                                                                                                             >  1:294310/1‑148 (MQ=255)
                                ttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAt                                                                                                             >  2:232386/1‑148 (MQ=255)
                                     acaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgttg                                                                                                        >  2:1867306/1‑148 (MQ=255)
                                        aGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                                                                                                     >  2:2412215/1‑148 (MQ=255)
                                        aGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                                                                                                     >  2:1441017/1‑148 (MQ=255)
                                         gCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCt                                                                                                    >  2:1913786/1‑148 (MQ=255)
                                          cTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTa                                                                                                   >  2:1106380/1‑148 (MQ=255)
                                          cTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTa                                                                                                   >  2:1033960/1‑148 (MQ=255)
                                          cTGTAAGAGCTTACTGAAAAAATTAACATCTCTTCCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTa                                                                                                   >  2:2240424/1‑148 (MQ=255)
                                           tGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAg                                                                                                  <  1:1390147/148‑1 (MQ=255)
                                           tGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAg                                                                                                  <  1:1390344/148‑1 (MQ=255)
                                            gTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                                                                                                 <  1:190542/148‑1 (MQ=255)
                                            gTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGATTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTAATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                                                                                                 >  1:3369070/1‑148 (MQ=255)
                                             tAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGc                                                                                                >  2:1185622/1‑148 (MQ=255)
                                                gagCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTa                                                                                             >  2:2130672/1‑148 (MQ=255)
                                                         gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACa                                                                                    >  1:2282428/1‑148 (MQ=255)
                                                               aTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcggc                                                                              >  2:2365522/1‑148 (MQ=255)
                                                                ttAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGgcggcg                                                                             >  1:3037960/1‑148 (MQ=255)
                                                                ttAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGgcggcg                                                                             >  2:3095142/1‑148 (MQ=255)
                                                                ttAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGgcggcg                                                                             >  2:1929642/1‑148 (MQ=255)
                                                                  aaCATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCa                                                                           >  2:3449500/1‑148 (MQ=255)
                                                                  aaCATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCa                                                                           >  2:2500475/1‑148 (MQ=255)
                                                                   aCATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAACTCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCAt                                                                          >  1:2166419/1‑148 (MQ=255)
                                                                      tctcTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGt                                                                       <  1:632091/148‑1 (MQ=255)
                                                                       ctcgTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGta                                                                      <  1:3030210/148‑1 (MQ=255)
                                                                        tctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGtat                                                                     <  1:1756486/148‑1 (MQ=255)
                                                                           tGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACa                                                                  <  2:1975479/148‑1 (MQ=255)
                                                                           tGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACa                                                                  >  1:2098935/1‑148 (MQ=255)
                                                                             cTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAgg                                                                <  2:1352854/148‑1 (MQ=255)
                                                                             cTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAgg                                                                <  2:1351057/148‑1 (MQ=255)
                                                                               aaGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCa                                                              <  2:187627/148‑1 (MQ=255)
                                                                               aaGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCa                                                              <  2:1046165/148‑1 (MQ=255)
                                                                                aGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCaa                                                             >  2:3496784/1‑148 (MQ=255)
                                                                                 gCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCaag                                                            >  1:2541547/1‑148 (MQ=255)
                                                                                     attaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAagaaga                                                        <  2:66250/148‑1 (MQ=255)
                                                                                      ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAt                                                       >  2:2697791/1‑148 (MQ=255)
                                                                                      ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGAGAATAATTATGAATAAATCGCAGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAatat                                                       >  1:3465729/1‑148 (MQ=255)
                                                                                       taAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAtt                                                      >  1:2116090/1‑148 (MQ=255)
                                                                                        aaaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTg                                                     <  2:623513/148‑1 (MQ=255)
                                                                                          aGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGaa                                                   >  1:3113921/1‑148 (MQ=255)
                                                                                          aGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGaa                                                   <  1:538867/148‑1 (MQ=255)
                                                                                           ggCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGaaa                                                  >  2:1233479/1‑148 (MQ=255)
                                                                                           ggCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGaaa                                                  >  1:2441913/1‑148 (MQ=255)
                                                                                             cTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAAcc                                                >  1:130644/1‑148 (MQ=255)
                                                                                               tataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCa                                              >  2:1503400/1‑148 (MQ=255)
                                                                                               tataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCa                                              <  1:407723/148‑1 (MQ=255)
                                                                                               tataACAACTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCa                                              <  1:2586081/148‑1 (MQ=255)
                                                                                                ataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAt                                             <  2:2417067/148‑1 (MQ=255)
                                                                                                  aaCACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtct                                           <  1:1225908/148‑1 (MQ=255)
                                                                                                     acCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGa                                        >  1:1710424/1‑148 (MQ=255)
                                                                                                      ccTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGaa                                       >  1:2268737/1‑148 (MQ=255)
                                                                                                      ccTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGaa                                       >  1:1793594/1‑148 (MQ=255)
                                                                                                       cTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAg                                      >  1:2073138/1‑148 (MQ=255)
                                                                                                        ttCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGa                                     >  1:3509400/1‑148 (MQ=255)
                                                                                                         tCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGAc                                    >  2:391784/1‑148 (MQ=255)
                                                                                                          cAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACa                                   <  1:101774/148‑1 (MQ=255)
                                                                                                          cAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACa                                   <  1:2048981/148‑1 (MQ=255)
                                                                                                           aGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACAt                                  <  2:2161089/148‑1 (MQ=255)
                                                                                                            ggcggcCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGACACCCATCTCGAAGACATg                                 >  2:3249063/1‑148 (MQ=255)
                                                                                                                 ccAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTcgc                            >  2:2301666/1‑148 (MQ=255)
                                                                                                                     tCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAg                        >  1:2648996/1‑148 (MQ=255)
                                                                                                                     tCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAg                        <  1:3139500/148‑1 (MQ=255)
                                                                                                                     tCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAg                        <  2:2547130/148‑1 (MQ=255)
                                                                                                                       cGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGc                      <  1:1930747/148‑1 (MQ=255)
                                                                                                                       cGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGc                      >  2:251987/1‑148 (MQ=255)
                                                                                                                        gCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCg                     <  1:2241342/148‑1 (MQ=255)
                                                                                                                         ccTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGa                    >  1:1067282/1‑148 (MQ=255)
                                                                                                                         ccTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGa                    <  1:3334156/148‑1 (MQ=255)
                                                                                                                          cTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGCTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGaa                   <  2:1730919/148‑1 (MQ=255)
                                                                                                                           tGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAAc                  >  2:1583789/1‑148 (MQ=255)
                                                                                                                               ttCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACgcgc              >  2:1644491/1‑148 (MQ=255)
                                                                                                                               ttCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACgcgc              >  2:3467572/1‑148 (MQ=255)
                                                                                                                                   ttttATTGATAATAATTATGAATAAATCGCTGGTAGGTGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGTCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACt          <  2:3120207/148‑1 (MQ=255)
                                                                                                                                    tttATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTc         >  2:1301434/1‑148 (MQ=255)
                                                                                                                                       aTGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCaaa      <  2:2282450/148‑1 (MQ=255)
                                                                                                                                       aTGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCaaa      >  1:58417/1‑148 (MQ=255)
                                                                                                                                        tGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAAc     >  1:478946/1‑148 (MQ=255)
                                                                                                                                          gATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTg   <  1:1652611/148‑1 (MQ=255)
                                                                                                                                          gATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTg   >  1:778893/1‑148 (MQ=255)
                                                                                                                                           ataataATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGa  <  1:391784/148‑1 (MQ=255)
                                                                                                                                           ataataATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGa  >  1:1569743/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                              
ACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGA  >  CP009273/1683961‑1684247

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: