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.1% 694.9 / 1.4 214E162D (GAA→GATcydAcytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base T (101/109);  total (102/112)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 214 total aligned reads displayed.

TCTTCTTCGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGG  >  CP009273/767254‑767542
                                                                                                                                                 |                                                                                                                                               
tcttcttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATAc                                                                                                                                               >  1:1393286/1‑148 (MQ=255)
tcttcttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATAc                                                                                                                                               >  2:2382305/1‑148 (MQ=255)
 cttcttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACt                                                                                                                                              >  2:1210206/1‑148 (MQ=255)
  ttcttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTa                                                                                                                                             <  1:649301/148‑1 (MQ=255)
    cttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATg                                                                                                                                           >  1:1501492/1‑148 (MQ=255)
      tcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCg                                                                                                                                         >  2:560230/1‑148 (MQ=255)
       cGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGt                                                                                                                                        <  1:3094889/148‑1 (MQ=255)
       cGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGt                                                                                                                                        <  1:2670378/148‑1 (MQ=255)
           tGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATgg                                                                                                                                    <  2:2553945/148‑1 (MQ=255)
              gATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGgaga                                                                                                                                 <  2:1744524/148‑1 (MQ=255)
              gATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGgaga                                                                                                                                 <  2:2231921/148‑1 (MQ=255)
              gATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGgaga                                                                                                                                 <  1:1356901/148‑1 (MQ=255)
               aTCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGAt                                                                                                                                >  1:1606697/1‑148 (MQ=255)
               aTCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGAt                                                                                                                                >  2:1365978/1‑148 (MQ=255)
                tCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATg                                                                                                                               >  1:1576978/1‑148 (MQ=255)
                tCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGAGTATGTAGATg                                                                                                                               >  1:1576969/1‑148 (MQ=255)
                 cGTCTGGGTACAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAACCCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATgg                                                                                                                              <  2:1831697/148‑1 (MQ=255)
                 cGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATgg                                                                                                                              >  2:116190/1‑148 (MQ=255)
                        gTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCtt                                                                                                                       <  1:694070/148‑1 (MQ=255)
                         tAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTtc                                                                                                                      <  2:3108024/148‑1 (MQ=255)
                          aaaGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTtct                                                                                                                     <  2:147501/148‑1 (MQ=255)
                           aaGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTtctc                                                                                                                    >  1:95978/1‑148 (MQ=255)
                           aaGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTtctc                                                                                                                    >  2:151829/1‑148 (MQ=255)
                                cAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGc                                                                                                               <  1:829084/148‑1 (MQ=255)
                                 aGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCt                                                                                                              <  2:1606697/148‑1 (MQ=255)
                                 aGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCt                                                                                                              <  1:1339683/148‑1 (MQ=255)
                                  gCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTg                                                                                                             >  1:354387/1‑148 (MQ=255)
                                   cATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTgg                                                                                                            >  1:625119/1‑148 (MQ=255)
                                      atGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGc                                                                                                         <  1:3165856/148‑1 (MQ=255)
                                      atGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGc                                                                                                         >  2:2025877/1‑148 (MQ=255)
                                       tgtgtgTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCt                                                                                                        <  2:1250794/148‑1 (MQ=255)
                                       tgtgtgTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCt                                                                                                        <  1:464310/148‑1 (MQ=255)
                                          gtgtCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTaa                                                                                                     >  1:3020468/1‑148 (MQ=255)
                                           tgtCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAAc                                                                                                    >  2:2877737/1‑148 (MQ=255)
                                             tCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAAccc                                                                                                  <  2:1883929/148‑1 (MQ=255)
                                                ccTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGt                                                                                               <  2:2913478/148‑1 (MQ=255)
                                                   ggctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGc                                                                                            <  1:2839431/148‑1 (MQ=255)
                                                    gctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCt                                                                                           >  2:1422339/1‑148 (MQ=255)
                                                    gctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCt                                                                                           <  2:3454909/148‑1 (MQ=255)
                                                     ctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTc                                                                                          >  2:3309268/1‑148 (MQ=255)
                                                        gtggCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAgg                                                                                       <  1:1066318/148‑1 (MQ=255)
                                                        gtggCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAgg                                                                                       <  1:1951900/148‑1 (MQ=255)
                                                         tggCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTTCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACATTGATAATATTCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGt                                                                                      >  1:605961/1‑148 (MQ=255)
                                                         tggCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGt                                                                                      >  2:3105993/1‑148 (MQ=255)
                                                              cTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAAt                                                                                 >  1:2026206/1‑148 (MQ=255)
                                                                   ttCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGtt                                                                            <  1:1428103/148‑1 (MQ=255)
                                                                      aaaCCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTcac                                                                         >  1:397921/1‑148 (MQ=255)
                                                                      aaaCCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTcac                                                                         >  1:3377308/1‑148 (MQ=255)
                                                                      aaaCCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTcac                                                                         >  1:1191737/1‑148 (MQ=255)
                                                                       aaCCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTcaca                                                                        <  2:812055/148‑1 (MQ=255)
                                                                        aCCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCacac                                                                       >  1:574267/1‑148 (MQ=255)
                                                                         ccTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACt                                                                      >  1:2777890/1‑148 (MQ=255)
                                                                          cTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTg                                                                     <  2:850450/148‑1 (MQ=255)
                                                                          cTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTg                                                                     <  2:444890/148‑1 (MQ=255)
                                                                                gCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGt                                                               <  2:3341415/148‑1 (MQ=255)
                                                                                  aCTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCt                                                             <  2:3267235/148‑1 (MQ=255)
                                                                                  aCTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCt                                                             >  1:2692077/1‑148 (MQ=255)
                                                                                  aCTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCt                                                             >  1:22325/1‑148 (MQ=255)
                                                                                      tgGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGtt                                                         <  2:1311026/148‑1 (MQ=255)
                                                                                        gATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAt                                                       >  1:2452347/1‑148 (MQ=255)
                                                                                        gATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAt                                                       >  1:1005073/1‑148 (MQ=255)
                                                                                         aTTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtg                                                      <  2:409914/148‑1 (MQ=255)
                                                                                         aTTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtg                                                      <  1:2221085/148‑1 (MQ=255)
                                                                                           tCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtgtg                                                    >  2:3441172/1‑148 (MQ=255)
                                                                                           tCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCGGGTTAtgtg                                                    >  2:1657154/1‑148 (MQ=255)
                                                                                            cTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGa                                                   <  2:3183736/148‑1 (MQ=255)
                                                                                              ggTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACt                                                 >  2:2614404/1‑148 (MQ=255)
                                                                                               gTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTg                                                <  1:714543/148‑1 (MQ=255)
                                                                                                ttGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTgg                                               >  1:728228/1‑148 (MQ=255)
                                                                                                ttGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTgg                                               >  1:1790896/1‑148 (MQ=255)
                                                                                                     aaCGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGa                                          <  1:2446670/148‑1 (MQ=255)
                                                                                                      aCGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGAt                                         >  1:3015995/1‑148 (MQ=255)
                                                                                                       cGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATg                                        >  1:2348741/1‑148 (MQ=255)
                                                                                                       cGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATg                                        >  2:2663035/1‑148 (MQ=255)
                                                                                                       cGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATg                                        <  1:1716498/148‑1 (MQ=255)
                                                                                                          cTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTc                                     <  1:772688/148‑1 (MQ=255)
                                                                                                            ggATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCAt                                   >  2:784470/1‑148 (MQ=255)
                                                                                                             gATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATc                                  >  2:759866/1‑148 (MQ=255)
                                                                                                                   aaaCCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGt                            >  1:3463160/1‑148 (MQ=255)
                                                                                                                    aaCCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTa                           <  2:2744118/148‑1 (MQ=255)
                                                                                                                     aCCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTAt                          <  2:2193951/148‑1 (MQ=255)
                                                                                                                      cccAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATc                         <  1:2736093/148‑1 (MQ=255)
                                                                                                                       ccAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCa                        >  1:541986/1‑148 (MQ=255)
                                                                                                                       ccAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGATCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATACTCGGTATCa                        <  2:467302/148‑1 (MQ=255)
                                                                                                                           tCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAgcgc                    >  2:3477566/1‑148 (MQ=255)
                                                                                                                            cgcgTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCa                   <  1:3011829/148‑1 (MQ=255)
                                                                                                                            cgcgTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCa                   >  1:2379670/1‑148 (MQ=255)
                                                                                                                            cgcgTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCa                   >  1:968226/1‑148 (MQ=255)
                                                                                                                               gTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATgg                <  2:2313602/148‑1 (MQ=255)
                                                                                                                                tCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGt               <  1:2967796/148‑1 (MQ=255)
                                                                                                                                 ccGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGta              >  1:1641612/1‑148 (MQ=255)
                                                                                                                                 ccGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGta              >  1:2569748/1‑148 (MQ=255)
                                                                                                                                   gATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGtata            >  2:647965/1‑148 (MQ=255)
                                                                                                                                    aTTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTatat           >  1:2312293/1‑148 (MQ=255)
                                                                                                                                     tttCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATg          <  2:3527085/148‑1 (MQ=255)
                                                                                                                                     tttCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATg          >  2:2719245/1‑148 (MQ=255)
                                                                                                                                          aCTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaa     <  2:391248/148‑1 (MQ=255)
                                                                                                                                          aCTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaa     <  2:3020468/148‑1 (MQ=255)
                                                                                                                                           cTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaaa    >  1:3050904/1‑148 (MQ=255)
                                                                                                                                             ttGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAgg  <  2:1039665/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TCTTCTTCGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGG  >  CP009273/767254‑767542

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: