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,3990AT97.8% 579.9 / 1.9 180E162D (GAA→GATcydAcytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base T (84/92);  total (84/96)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-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 180 total aligned reads displayed.

TCTTCTTCGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGG  >  CP009273/767254‑767542
                                                                                                                                                 |                                                                                                                                               
tcttcttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATAc                                                                                                                                               >  2:432804/1‑148 (MQ=255)
tcttcttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATAc                                                                                                                                               >  2:2870218/1‑148 (MQ=255)
 cttcttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACt                                                                                                                                              <  2:2045809/148‑1 (MQ=255)
  ttcttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCTCGTCCGATTTCAACTTTGATACTa                                                                                                                                             <  2:2584161/148‑1 (MQ=255)
  ttcttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTa                                                                                                                                             <  2:2584240/148‑1 (MQ=255)
  ttcttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTa                                                                                                                                             >  2:3545222/1‑148 (MQ=255)
  ttcgtcGGTTGGGTTCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCACACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTa                                                                                                                                             <  1:1430972/148‑1 (MQ=255)
    cttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATg                                                                                                                                           <  2:2249037/148‑1 (MQ=255)
    cttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATg                                                                                                                                           >  2:496513/1‑148 (MQ=255)
    cttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATg                                                                                                                                           <  2:3360387/148‑1 (MQ=255)
     ttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGc                                                                                                                                          >  1:1939935/1‑148 (MQ=255)
       cGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGt                                                                                                                                        <  2:2080117/148‑1 (MQ=255)
       cGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGt                                                                                                                                        <  2:2761010/148‑1 (MQ=255)
       cGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGt                                                                                                                                        <  1:3062374/148‑1 (MQ=255)
        ggTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTa                                                                                                                                       <  1:775432/148‑1 (MQ=255)
         gTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTAt                                                                                                                                      <  1:827241/148‑1 (MQ=255)
           tGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATgg                                                                                                                                    >  2:3636878/1‑148 (MQ=255)
              gATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGgaga                                                                                                                                 >  2:1544523/1‑148 (MQ=255)
              gATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGgaga                                                                                                                                 <  2:1209086/148‑1 (MQ=255)
               aTCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGAt                                                                                                                                >  1:153953/1‑148 (MQ=255)
               aTCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGAt                                                                                                                                >  1:2805807/1‑148 (MQ=255)
                tCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATg                                                                                                                               >  2:80233/1‑148 (MQ=255)
                   tCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTg                                                                                                                            <  1:3036312/148‑1 (MQ=255)
                      gggTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGc                                                                                                                         <  1:3376239/148‑1 (MQ=255)
                      gggTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGc                                                                                                                         >  2:1473378/1‑148 (MQ=255)
                           aaGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTtctc                                                                                                                    >  2:845977/1‑148 (MQ=255)
                            aggtCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTcc                                                                                                                   >  1:3156260/4‑148 (MQ=255)
                             gTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCg                                                                                                                  <  1:252729/148‑1 (MQ=255)
                             gTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCg                                                                                                                  <  2:3402064/148‑1 (MQ=255)
                             gTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCg                                                                                                                  <  1:382092/148‑1 (MQ=255)
                             gTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGAGTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCg                                                                                                                  <  2:3402040/148‑1 (MQ=255)
                               tCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAg                                                                                                                >  2:661261/1‑148 (MQ=255)
                                  gCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTg                                                                                                             >  1:3801586/1‑148 (MQ=255)
                                   cATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTgg                                                                                                            <  2:1136449/148‑1 (MQ=255)
                                   cATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTgg                                                                                                            >  1:1079071/1‑148 (MQ=255)
                                      atGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGc                                                                                                         >  1:2514306/1‑148 (MQ=255)
                                      atGTGTGACACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGc                                                                                                         >  2:87678/1‑148 (MQ=255)
                                        gtgtgtCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCtt                                                                                                       >  2:236168/1‑148 (MQ=255)
                                         tgtgtCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTa                                                                                                      <  2:757134/148‑1 (MQ=255)
                                          gtgtCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTaa                                                                                                     <  1:2698704/148‑1 (MQ=255)
                                           tgtCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAAc                                                                                                    <  2:2443458/148‑1 (MQ=255)
                                            gtCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAAcc                                                                                                   <  2:1338598/148‑1 (MQ=255)
                                             tCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAAccc                                                                                                  >  1:1371889/1‑148 (MQ=255)
                                             tCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAAAGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAAccc                                                                                                  >  1:1370649/1‑148 (MQ=255)
                                               aCCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCgg                                                                                                <  1:3671422/148‑1 (MQ=255)
                                               aCCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCgg                                                                                                <  1:1779901/148‑1 (MQ=255)
                                                 ctggctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGtt                                                                                              >  2:238939/1‑148 (MQ=255)
                                                  tggctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTg                                                                                             >  2:2414771/1‑148 (MQ=255)
                                                  tggctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTg                                                                                             >  1:774818/1‑148 (MQ=255)
                                                      tggtggCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCa                                                                                         <  1:2238855/148‑1 (MQ=255)
                                                             gcTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGaaa                                                                                  <  1:91413/148‑1 (MQ=255)
                                                                cGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTc                                                                               >  2:914628/1‑148 (MQ=255)
                                                                cGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTc                                                                               >  2:415654/1‑148 (MQ=255)
                                                                   ttCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGtt                                                                            <  2:294332/148‑1 (MQ=255)
                                                                      aaaCCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTcac                                                                         <  2:2745398/148‑1 (MQ=255)
                                                                        aCCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCacac                                                                       >  1:1045235/1‑148 (MQ=255)
                                                                        aCCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCacac                                                                       >  1:3445594/1‑148 (MQ=255)
                                                                                 cACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTc                                                              >  1:677093/1‑148 (MQ=255)
                                                                                    tgtgGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTgg                                                           >  1:2997537/1‑148 (MQ=255)
                                                                                       ggATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTa                                                        <  1:3590309/148‑1 (MQ=255)
                                                                                         aTTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtg                                                      <  1:661007/148‑1 (MQ=255)
                                                                                         aTTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtg                                                      <  1:713372/148‑1 (MQ=255)
                                                                                         aTTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtg                                                      <  2:38708/148‑1 (MQ=255)
                                                                                         aTTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtg                                                      <  1:2098229/148‑1 (MQ=255)
                                                                                           tCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtgtg                                                    <  2:1612216/148‑1 (MQ=255)
                                                                                             tGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGAc                                                  <  2:3394365/148‑1 (MQ=255)
                                                                                               gTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTg                                                >  2:1653283/1‑148 (MQ=255)
                                                                                                ttGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTgg                                               >  2:3266706/1‑148 (MQ=255)
                                                                                                ttGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTgg                                               >  2:3145650/1‑148 (MQ=255)
                                                                                                     aaCGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGa                                          <  1:496513/148‑1 (MQ=255)
                                                                                                      aCGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGAt                                         >  1:298392/1‑148 (MQ=255)
                                                                                                       cGGCCGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATg                                        <  2:417562/148‑1 (MQ=255)
                                                                                                          cTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACAATGTAGCGTCTGGTTATGTGACTGGCGCGATGTTc                                     >  2:640393/1‑148 (MQ=255)
                                                                                                           tGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACGATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCa                                    <  1:2939820/148‑1 (MQ=255)
                                                                                                                  aaacCCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCgg                             <  2:2928412/148‑1 (MQ=255)
                                                                                                                    aaCCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTa                           >  2:2316365/1‑148 (MQ=255)
                                                                                                                      cccAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATc                         >  1:1366224/1‑148 (MQ=255)
                                                                                                                            cgcgTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCa                   >  1:2699502/1‑148 (MQ=255)
                                                                                                                             gcgTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCAt                  <  1:87678/148‑1 (MQ=255)
                                                                                                                             gcgTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCAt                  <  2:1233953/148‑1 (MQ=255)
                                                                                                                              cgTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATg                 <  1:1544523/148‑1 (MQ=255)
                                                                                                                              cgTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATg                 <  1:2311934/148‑1 (MQ=255)
                                                                                                                                tCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGt               >  1:3675004/1‑148 (MQ=255)
                                                                                                                                  cGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGtat             <  1:1278994/148‑1 (MQ=255)
                                                                                                                                  cGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGtat             <  2:3472674/148‑1 (MQ=255)
                                                                                                                                  cGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGtat             >  1:1383206/1‑148 (MQ=255)
                                                                                                                                    aTTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTatat           >  2:1993359/1‑148 (MQ=255)
                                                                                                                                    aTTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTatat           >  2:2286491/1‑148 (MQ=255)
                                                                                                                                    aTTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTatat           >  2:363223/1‑148 (MQ=255)
                                                                                                                                      ttCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGc         >  1:2491366/1‑148 (MQ=255)
                                                                                                                                        cAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTg       <  2:1282358/148‑1 (MQ=255)
                                                                                                                                          aCTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaa     <  1:2101027/148‑1 (MQ=255)
                                                                                                                                          aCTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaa     >  1:1968521/1‑148 (MQ=255)
                                                                                                                                          aCTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaa     <  2:1323602/148‑1 (MQ=255)
                                                                                                                                           cTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaaa    >  1:3494468/1‑148 (MQ=255)
                                                                                                                                           cTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaaa    <  2:763344/148‑1 (MQ=255)
                                                                                                                                           cTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaaa    <  2:2831731/148‑1 (MQ=255)
                                                                                                                                             ttGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAgg  >  1:741825/1‑148 (MQ=255)
                                                                                                                                             ttGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAgg  >  1:2754670/1‑148 (MQ=255)
                                                                                                                                             ttGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAgg  >  2:108719/1‑148 (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: