Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 767,399 A→T 100% E162D (GAA→GAT cydA → cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273767,3990AT98.4% 617.7 / ‑0.2 190E162D (GAA→GATcydAcytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (100/87);  total (102/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-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 190 total aligned reads displayed.

CTTCTTCGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGT  >  CP009273/767255‑767543
                                                                                                                                                |                                                                                                                                                
   cttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATg                                                                                                                                            <  2:1709095/148‑1 (MQ=255)
     tcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCg                                                                                                                                          >  2:660189/1‑148 (MQ=255)
      cGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGt                                                                                                                                         <  2:1132088/148‑1 (MQ=255)
       ggTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTa                                                                                                                                        >  2:2243808/1‑148 (MQ=255)
            ggATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCAGGTTGCGAACGGCTGGATGCAAAACACAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGgag                                                                                                                                   >  2:3045124/1‑148 (MQ=255)
                  tCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTg                                                                                                                             >  2:1098638/1‑148 (MQ=255)
                  tCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCGGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTg                                                                                                                             >  2:2595587/1‑148 (MQ=255)
                    tGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAg                                                                                                                           <  1:1519663/148‑1 (MQ=255)
                     gggTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGc                                                                                                                          <  1:2281272/148‑1 (MQ=255)
                     gggTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGc                                                                                                                          <  2:1977468/148‑1 (MQ=255)
                      ggTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCt                                                                                                                         >  2:2419980/1‑148 (MQ=255)
                           aGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTcc                                                                                                                    >  1:1453148/1‑148 (MQ=255)
                              tCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAg                                                                                                                 <  1:952369/148‑1 (MQ=255)
                               catcATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTGCGCTCTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGc                                                                                                                <  1:3291486/145‑1 (MQ=255)
                               cAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGc                                                                                                                <  1:1380656/148‑1 (MQ=255)
                                 gCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTg                                                                                                              >  1:1811500/1‑148 (MQ=255)
                                  cATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTgg                                                                                                             >  2:3515904/1‑148 (MQ=255)
                                   atatGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGt                                                                                                            >  1:2616461/1‑148 (MQ=255)
                                    tatGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTg                                                                                                           <  1:2468958/148‑1 (MQ=255)
                                    tatGTGTGTCACCTGGCTGGTGGCGCTAGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACTGCTGGATGCAAAACACAATCGCGTCCGATTTAAACATTGATACTATGCGTATGGAGATGGTGAGCTTCTACGAGATtgtg                                                                                                           >  1:348366/1‑148 (MQ=255)
                                    tatGTGTGGCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTg                                                                                                           <  1:687248/148‑1 (MQ=255)
                                     atGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGc                                                                                                          <  1:2011123/148‑1 (MQ=255)
                                     atGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGc                                                                                                          >  1:2231324/1‑148 (MQ=255)
                                     atGTGTGTCACCTGGCTGGTGGCGCTCCGTTCAAAACAGTCCGCACTGTGGATTCTGGTTGCGAAAGGCTGGATGCAAAACCAAATCGCGTCCGATTTCAACTTTGATACTATGCGAATGGAGATGGTTAGCTTCTACGAGCTGGTGc                                                                                                          >  1:158784/1‑148 (MQ=255)
                                       gtgtgtCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTAGGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCtt                                                                                                        >  2:128730/1‑148 (MQ=255)
                                        tgtgtCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTa                                                                                                       <  2:2353457/148‑1 (MQ=255)
                                        tgtgtCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTa                                                                                                       >  2:1238444/1‑148 (MQ=255)
                                         gtgtCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTaa                                                                                                      <  1:2307662/148‑1 (MQ=255)
                                                  ggctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGc                                                                                             >  1:1988800/1‑148 (MQ=255)
                                                   gctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCt                                                                                            >  2:3383812/1‑148 (MQ=255)
                                                   gctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCt                                                                                            >  1:3024383/1‑148 (MQ=255)
                                                   gctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCt                                                                                            >  1:2133670/1‑148 (MQ=255)
                                                   gctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCt                                                                                            >  2:2991858/1‑148 (MQ=255)
                                                    ctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTc                                                                                           <  1:2043277/148‑1 (MQ=255)
                                                          gcgcTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGa                                                                                     >  1:290465/1‑148 (MQ=255)
                                                           cgcTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGaa                                                                                    <  1:1714162/148‑1 (MQ=255)
                                                             cTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAAt                                                                                  <  1:698981/148‑1 (MQ=255)
                                                              tCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAAtt                                                                                 >  1:1746668/1‑148 (MQ=255)
                                                                ggTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCg                                                                               >  2:1784577/1‑148 (MQ=255)
                                                                ggTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCg                                                                               >  2:3199528/1‑148 (MQ=255)
                                                                ggTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCg                                                                               <  1:1655166/148‑1 (MQ=255)
                                                                 gTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGt                                                                              >  1:2089573/1‑148 (MQ=255)
                                                                 gTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGt                                                                              <  1:3499107/148‑1 (MQ=255)
                                                                 gTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGt                                                                              <  2:1459541/148‑1 (MQ=255)
                                                                 gTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGt                                                                              <  1:660189/148‑1 (MQ=255)
                                                                  ttCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGtt                                                                             >  1:1197101/1‑148 (MQ=255)
                                                                   tCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTc                                                                            <  2:3152019/148‑1 (MQ=255)
                                                                    cAAACCTGTCCGCACTGTGGAGTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTca                                                                           <  2:1492511/148‑1 (MQ=255)
                                                                     aaaCCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTcac                                                                          <  2:2568442/148‑1 (MQ=255)
                                                                      aaCCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTcaca                                                                         <  1:817625/148‑1 (MQ=255)
                                                                        ccTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACt                                                                       >  2:238870/1‑148 (MQ=255)
                                                                        ccTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACt                                                                       >  2:2545302/1‑148 (MQ=255)
                                                                        ccTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACt                                                                       >  2:3158587/1‑148 (MQ=255)
                                                                          tGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGt                                                                     <  2:320272/148‑1 (MQ=255)
                                                                            tCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAg                                                                   <  2:288854/148‑1 (MQ=255)
                                                                               gCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGt                                                                <  1:1904599/148‑1 (MQ=255)
                                                                               gCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGt                                                                <  1:2165392/148‑1 (MQ=255)
                                                                                  cTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTg                                                             >  1:601678/1‑148 (MQ=255)
                                                                                  cTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTg                                                             <  2:2854051/148‑1 (MQ=255)
                                                                                   tgtgGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCAGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTgg                                                            >  2:1317649/1‑148 (MQ=255)
                                                                                   tgtgGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGCGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTgg                                                            >  2:1317668/1‑148 (MQ=255)
                                                                                     tgGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGtt                                                          >  1:2278471/1‑148 (MQ=255)
                                                                                     tgGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGtt                                                          >  1:1255857/1‑148 (MQ=255)
                                                                                     tgGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGtt                                                          >  1:491218/1‑148 (MQ=255)
                                                                                     tgGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGtt                                                          >  1:2278959/1‑148 (MQ=255)
                                                                                      ggATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTa                                                         >  1:1000945/1‑148 (MQ=255)
                                                                                      ggATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTa                                                         >  1:1000904/1‑148 (MQ=255)
                                                                                      ggATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTa                                                         >  1:1614408/1‑148 (MQ=255)
                                                                                          tCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtgtg                                                     <  2:1582329/148‑1 (MQ=255)
                                                                                          tCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtgtg                                                     <  1:2072007/148‑1 (MQ=255)
                                                                                          tCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtgtg                                                     >  1:1109929/1‑148 (MQ=255)
                                                                                           cTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGa                                                    <  2:2078545/148‑1 (MQ=255)
                                                                                           cTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGa                                                    >  1:211004/1‑148 (MQ=255)
                                                                                               ttGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTgg                                                >  1:1946476/1‑148 (MQ=255)
                                                                                               ttGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTgg                                                <  1:3157365/148‑1 (MQ=255)
                                                                                               ttGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTgg                                                >  1:2743628/1‑148 (MQ=255)
                                                                                                      cGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATg                                         <  2:827918/148‑1 (MQ=255)
                                                                                                        gCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGtt                                       >  2:1452119/1‑148 (MQ=255)
                                                                                                         cTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTc                                      >  1:3031324/1‑148 (MQ=255)
                                                                                                         cTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTc                                      >  1:3031180/1‑148 (MQ=255)
                                                                                                            gATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATc                                   >  1:423964/1‑148 (MQ=255)
                                                                                                             aTGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATcc                                  <  1:1854975/148‑1 (MQ=255)
                                                                                                               gCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTc                                >  1:3413798/1‑148 (MQ=255)
                                                                                                               gCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTc                                >  1:1044982/1‑148 (MQ=255)
                                                                                                                    aCCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTAt                           >  2:2475262/1‑148 (MQ=255)
                                                                                                                    aCCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTAt                           <  2:2639905/148‑1 (MQ=255)
                                                                                                                       cAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAg                        <  2:1993182/148‑1 (MQ=255)
                                                                                                                         aTCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAgcg                      >  2:2055080/1‑148 (MQ=255)
                                                                                                                          tCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAgcgc                     >  2:2696426/1‑148 (MQ=255)
                                                                                                                                ccGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGta               <  1:1231480/148‑1 (MQ=255)
                                                                                                                                 cGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGtat              <  1:1446540/148‑1 (MQ=255)
                                                                                                                                  gATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGtata             <  2:902272/148‑1 (MQ=255)
                                                                                                                                         aCTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaa      >  1:2216828/1‑148 (MQ=255)
                                                                                                                                         aCTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaa      <  2:1811500/148‑1 (MQ=255)
                                                                                                                                          cTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaaa     >  2:709812/1‑148 (MQ=255)
                                                                                                                                          cTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaaa     >  2:901425/1‑148 (MQ=255)
                                                                                                                                           tttGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAg    >  2:3333326/1‑148 (MQ=255)
                                                                                                                                           tttGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAg    <  1:2419980/148‑1 (MQ=255)
                                                                                                                                             tGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGt  >  2:1649742/1‑148 (MQ=255)
                                                                                                                                             tGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGt  >  2:1170897/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CTTCTTCGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGGT  >  CP009273/767255‑767543

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: