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.2% 729.1 / 1.7 223E162D (GAA→GATcydAcytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (4/0);  new base T (113/106);  total (117/106)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-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 224 total aligned reads displayed.

TTCTTCTTCGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGG  >  CP009273/767253‑767542
                                                                                                                                                  |                                                                                                                                               
    tcttcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTAt                                                                                                                                            >  1:698674/1‑148 (MQ=255)
       tcGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCg                                                                                                                                         <  1:2509630/148‑1 (MQ=255)
        cGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGt                                                                                                                                        >  2:2740962/1‑148 (MQ=255)
            tGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGAACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTATGATACTATTCGTATgt                                                                                                                                    >  2:3151581/1‑147 (MQ=255)
              ggATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATg                                                                                                                                           >  2:378645/1‑139 (MQ=255)
                 tCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATg                                                                                                                               <  2:1107008/148‑1 (MQ=255)
                 tCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATg                                                                                                                               >  1:626598/1‑148 (MQ=255)
                 tCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATg                                                                                                                               <  1:1467485/148‑1 (MQ=255)
                       gggTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGc                                                                                                                         <  2:1362400/148‑1 (MQ=255)
                        ggTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCt                                                                                                                        <  1:99895/148‑1 (MQ=255)
                          tAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTtc                                                                                                                      <  2:2457563/148‑1 (MQ=255)
                           aaaGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTtct                                                                                                                     <  1:540756/148‑1 (MQ=255)
                             aGTTCAGCTTATGTGTGTCACCTGGCTGGTGGCGCTCGGGTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCATTCGCGTCCGATTTCAACTTTGATAGTATGCGTATGGAGTTGGTGAGCTTCTcc                                                                                                                   <  1:1427152/148‑1 (MQ=255)
                             aGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTcc                                                                                                                   >  1:1820121/1‑148 (MQ=255)
                             aGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTcc                                                                                                                   <  1:3621/148‑1 (MQ=255)
                              gTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCg                                                                                                                  <  2:2300703/148‑1 (MQ=255)
                               ttCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGa                                                                                                                 >  1:1700238/1‑148 (MQ=255)
                               ttCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGa                                                                                                                 >  1:1404985/1‑148 (MQ=255)
                               ttCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGa                                                                                                                 >  1:3461830/1‑148 (MQ=255)
                                tCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAg                                                                                                                <  2:700523/148‑1 (MQ=255)
                                    cATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTgg                                                                                                            >  2:2846332/1‑148 (MQ=255)
                                    cATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTgg                                                                                                            >  2:598283/1‑148 (MQ=255)
                                    cATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTgg                                                                                                            >  2:1621478/1‑148 (MQ=255)
                                       atGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGc                                                                                                         <  1:2666416/148‑1 (MQ=255)
                                         gtgtgtCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGAGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCtt                                                                                                       >  2:1878505/1‑148 (MQ=255)
                                                aCCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCgg                                                                                                <  1:446377/148‑1 (MQ=255)
                                                   tggctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTg                                                                                             >  2:2863110/1‑148 (MQ=255)
                                                     gctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGAGGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCt                                                                                           >  2:3436836/1‑148 (MQ=255)
                                                      ctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTc                                                                                          >  2:2530728/1‑148 (MQ=255)
                                                      ctggTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTc                                                                                          <  2:2167732/148‑1 (MQ=255)
                                                        ggtggCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAg                                                                                        <  1:2650451/148‑1 (MQ=255)
                                                         gtggCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAgg                                                                                       >  2:1186778/1‑148 (MQ=255)
                                                         gtggCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAgg                                                                                       <  2:3261731/148‑1 (MQ=255)
                                                               cTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAAt                                                                                 <  2:67779/148‑1 (MQ=255)
                                                               cTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAAt                                                                                 <  2:297687/148‑1 (MQ=255)
                                                               cTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAAt                                                                                 <  1:1608575/148‑1 (MQ=255)
                                                               cTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAAt                                                                                 <  1:3348295/148‑1 (MQ=255)
                                                                 cGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTc                                                                               <  2:1575773/148‑1 (MQ=255)
                                                                  ggTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCg                                                                              <  2:3242081/148‑1 (MQ=255)
                                                                   gTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGt                                                                             >  1:2893610/1‑148 (MQ=255)
                                                                    ttCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGtt                                                                            <  1:2562079/148‑1 (MQ=255)
                                                                    ttCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGtt                                                                            >  1:1913222/1‑148 (MQ=255)
                                                                     tCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTCTGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTc                                                                           >  2:3321856/1‑148 (MQ=255)
                                                                     tCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAGCCCGGTTGCTCAGGTGAAATTCGTTc                                                                           <  2:2439931/148‑1 (MQ=255)
                                                                        aaCCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTcaca                                                                        >  2:2768490/1‑148 (MQ=255)
                                                                           cTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTg                                                                     <  2:2697137/148‑1 (MQ=255)
                                                                            tGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGt                                                                    >  1:1210261/1‑148 (MQ=255)
                                                                             gTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTa                                                                   <  2:282210/148‑1 (MQ=255)
                                                                             gTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTa                                                                   >  1:1636429/1‑148 (MQ=255)
                                                                              tCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAg                                                                  >  1:961702/1‑148 (MQ=255)
                                                                                 ggaCTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGt                                                               <  2:3038630/146‑1 (MQ=255)
                                                                                 gCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGt                                                               <  2:1965725/148‑1 (MQ=255)
                                                                                 gCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGt                                                               <  2:793495/148‑1 (MQ=255)
                                                                                 gCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGt                                                               <  1:3349914/148‑1 (MQ=255)
                                                                                   aCTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCt                                                             <  1:217955/148‑1 (MQ=255)
                                                                                     tgtgGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTgg                                                           >  2:1689413/1‑148 (MQ=255)
                                                                                     tgtgGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTgg                                                           <  1:2291663/148‑1 (MQ=255)
                                                                                     tgtgGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTgg                                                           <  1:2109861/148‑1 (MQ=255)
                                                                                       tgGATTCTGGTTGCGAACGGGTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGtt                                                         >  2:1347085/1‑148 (MQ=255)
                                                                                       tgGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGtt                                                         >  2:2606327/1‑148 (MQ=255)
                                                                                        ggATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTa                                                        <  2:1136791/148‑1 (MQ=255)
                                                                                        ggATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTa                                                        <  2:2265977/148‑1 (MQ=255)
                                                                                          aTTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtg                                                      >  1:40078/1‑148 (MQ=255)
                                                                                          aTTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtg                                                      <  1:3192472/148‑1 (MQ=255)
                                                                                            tCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTTTGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtgtg                                                    <  1:3151581/148‑1 (MQ=255)
                                                                                            tCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTAtgtg                                                    >  1:1012202/1‑148 (MQ=255)
                                                                                               ggTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTCTGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACt                                                 <  2:2773084/148‑1 (MQ=255)
                                                                                                  tGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGgc                                              >  2:3285351/1‑148 (MQ=255)
                                                                                                     gAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGcgcg                                           >  1:2724028/1‑148 (MQ=255)
                                                                                                      aaCGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGa                                          >  1:1787083/1‑148 (MQ=255)
                                                                                                         ggCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGt                                       >  1:1187246/1‑148 (MQ=255)
                                                                                                           cTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTc                                     <  2:1296609/148‑1 (MQ=255)
                                                                                                             ggATGCAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCAt                                   <  2:1737959/148‑1 (MQ=255)
                                                                                                                  cAAAACCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCg                              >  2:356124/1‑148 (MQ=255)
                                                                                                                     aaCCCAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTa                           >  2:715809/1‑148 (MQ=255)
                                                                                                                       cccAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATc                         >  2:917703/1‑148 (MQ=255)
                                                                                                                       cccAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATc                         >  1:2012726/1‑148 (MQ=255)
                                                                                                                       cccAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATc                         >  1:2565768/1‑148 (MQ=255)
                                                                                                                         cAATCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAg                       >  2:2948344/1‑148 (MQ=255)
                                                                                                                          aaTCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAgc                      >  1:512111/1‑148 (MQ=255)
                                                                                                                          aaTCGCGTCCGATTTCAACTTTGATACTATGCGTATGGAGATGGGGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAgc                      >  1:1102115/1‑148 (MQ=255)
                                                                                                                                  ccGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGta              >  1:263191/1‑148 (MQ=255)
                                                                                                                                   cGATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGtat             >  1:1138606/1‑148 (MQ=255)
                                                                                                                                    gATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGtata            <  1:3400602/148‑1 (MQ=255)
                                                                                                                                    gATTTCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGtata            >  2:2189403/1‑148 (MQ=255)
                                                                                                                                      tttCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTCTATg          >  2:1467263/1‑148 (MQ=255)
                                                                                                                                      tttCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATg          <  1:1366167/148‑1 (MQ=255)
                                                                                                                                      tttCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATg          >  1:3437261/1‑148 (MQ=255)
                                                                                                                                      tttCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATg          <  1:258349/148‑1 (MQ=255)
                                                                                                                                      tttCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATg          <  2:1808470/148‑1 (MQ=255)
                                                                                                                                       ttCAACTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGc         >  1:1596056/1‑148 (MQ=255)
                                                                                                                                           aCTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaa     <  2:476461/148‑1 (MQ=255)
                                                                                                                                           aCTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaa     <  2:519277/148‑1 (MQ=255)
                                                                                                                                            cTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaaa    >  1:539155/1‑148 (MQ=255)
                                                                                                                                            cTTTGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGaaa    <  2:1858100/148‑1 (MQ=255)
                                                                                                                                             tttGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAg   <  1:77881/148‑1 (MQ=255)
                                                                                                                                              ttGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAgg  >  1:2251302/1‑148 (MQ=255)
                                                                                                                                              ttGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAgg  <  2:347390/148‑1 (MQ=255)
                                                                                                                                              ttGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAgg  >  1:484982/1‑148 (MQ=255)
                                                                                                                                              ttGATACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAgg  <  2:1924506/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TTCTTCTTCGGTTGGGATCGTCTGGGTAAAGTTCAGCATATGTGTGTCACCTGGCTGGTGGCGCTCGGTTCAAACCTGTCCGCACTGTGGATTCTGGTTGCGAACGGCTGGATGCAAAACCCAATCGCGTCCGATTTCAACTTTGAAACTATGCGTATGGAGATGGTGAGCTTCTCCGAGCTGGTGCTTAACCCGGTTGCTCAGGTGAAATTCGTTCACACTGTAGCGTCTGGTTATGTGACTGGCGCGATGTTCATCCTCGGTATCAGCGCATGGTATATGCTGAAAGG  >  CP009273/767253‑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: