Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 444,779 C→T 72.5% W440* (TGG→TAG)  cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273444,7790CT72.5% 100.7 / 35.5 74W440* (TGG→TAG) cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (8/13);  new base T (21/31);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCGAG  >  CP009273/444633‑444925
                                                                                                                                                  |                                                                                                                                                  
ggTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGctac                                                                                                                                                  <  2:860970/149‑1 (MQ=255)
 gTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCACAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGctacc                                                                                                                                                 >  1:823661/1‑149 (MQ=255)
  tCAATCTGCTGGCTCAAACGAGGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGccacca                                                                                                                                                <  2:693774/149‑1 (MQ=255)
  tCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGctaccc                                                                                                                                                >  2:81355/1‑148 (MQ=255)
  tCAATCTGCGGGCACAAACTACGGGTGATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGGTTGAAACCGAACGCTTTAGGctacca                                                                                                                                                >  2:1037895/1‑149 (MQ=255)
     aTCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGGACGAGCGGTTACCCCAGGTTTCGTGCAGTTTGAAACCGAACGCTGTAGGCTACAAGTa                                                                                                                                             >  2:996623/1‑149 (MQ=255)
         gctgGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCACCGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGCAAGCTTTAGGCCACCAGTAGGTc                                                                                                                                         >  1:817104/1‑149 (MQ=255)
         gctgGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTc                                                                                                                                         >  2:291776/1‑149 (MQ=255)
          ctgGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCAGACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCa                                                                                                                                        <  2:910323/149‑1 (MQ=255)
             gCTCAAACGACGGGTCATGCCCCTGAAGCCCAGCGAATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATcc                                                                                                                                     <  1:1037335/149‑1 (MQ=255)
              cTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATccc                                                                                                                                    >  1:586417/1‑149 (MQ=255)
               tCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCt                                                                                                                                   >  2:541703/1‑149 (MQ=255)
                cAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGCCATCCCTg                                                                                                                                  <  1:905900/149‑1 (MQ=255)
                    cgcggggTCATGTCCATGAAGCCCAGCGCCTACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGaa                                                                                                                              <  2:1043471/145‑1 (MQ=255)
                    cgacgGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGaa                                                                                                                              <  2:947143/149‑1 (MQ=255)
                    cgacgGCTCATGCGCATGAAGCCCAGCGCAGACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGaa                                                                                                                              <  2:966354/149‑1 (MQ=255)
                     gacgGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGCACCAGAACGCGCGTTCACCCCAGGTCTCGTTCAGCTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAag                                                                                                                             >  1:867113/1‑149 (MQ=255)
                       cgGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGACATCCAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAagca                                                                                                                           <  1:231482/149‑1 (MQ=255)
                        gggTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAagcag                                                                                                                          <  1:317753/149‑1 (MQ=255)
                         ggTCATGCCCATGAAGCCCAGCGCATACAGTCGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCTGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAagcagc                                                                                                                         <  1:831497/149‑1 (MQ=255)
                          gTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGcc                                                                                                                        >  1:113101/1‑149 (MQ=255)
                            cATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGa                                                                                                                      <  1:1051734/149‑1 (MQ=255)
                             aggCCCATGAAGCCCAGCGCAAACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACAAGAACGCGCGCTTACCCCTGGGTTCGTTCAGGTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGaa                                                                                                                     <  1:555763/147‑1 (MQ=255)
                             aTGCCCATAAAGCCCAGCGCATACAGGGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAATGCGCGTTTACCCCAGGTTTCTTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGaa                                                                                                                     <  1:120986/149‑1 (MQ=255)
                                cccATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGAc                                                                                                                  <  2:251867/149‑1 (MQ=255)
                                cccATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGAc                                                                                                                  <  2:884965/149‑1 (MQ=255)
                                      aaGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGCGGGCCAcgcc                                                                                                            >  1:529585/1‑149 (MQ=255)
                                      aaGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCAcgcc                                                                                                            <  2:1173929/149‑1 (MQ=255)
                                         cccAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACgccgcc                                                                                                         >  2:45715/1‑149 (MQ=255)
                                                  tACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCAc                                                                                                >  1:255751/1‑149 (MQ=255)
                                                    cggTGGCAAAAAGGCAACGAAGAATCCGATGAGCCAGAATCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCGTTATGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGt                                                                                              <  2:177314/147‑1 (MQ=255)
                                                      gTGGCATAAAGGCAACGAAGCAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTGCGTTTAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTa                                                                                            <  2:396270/149‑1 (MQ=255)
                                                        ggCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTCCCCCCGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCCTCCCTGCGAAGCAGCCGAAGACGACGCCGCCGATGATCACGTTATg                                                                                          <  2:421918/149‑1 (MQ=255)
                                                        ggCAGAAAGGCAACGAAGAAGCCGATGAGCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATg                                                                                          <  1:573753/149‑1 (MQ=255)
                                                           aTAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGaa                                                                                       <  2:954383/149‑1 (MQ=255)
                                                           aTAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGaa                                                                                       <  1:1237703/149‑1 (MQ=255)
                                                           aTAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGaa                                                                                       >  2:1039364/1‑149 (MQ=255)
                                                            tAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAg                                                                                      >  2:293645/1‑149 (MQ=255)
                                                             aaaGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGt                                                                                     <  2:831539/149‑1 (MQ=255)
                                                                ggCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAGGTgcg                                                                                  <  2:823661/149‑1 (MQ=255)
                                                                ggCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTgcg                                                                                  >  1:388650/1‑149 (MQ=255)
                                                                   aaCGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCaa                                                                               >  1:966118/1‑149 (MQ=255)
                                                                      gaagaaGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCCGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCCAAGCAGCCGAAGACCACGCCGCCGATGAGCGCGTTCTGGAAGTGCGCAATCa                                                                            >  2:292781/1‑149 (MQ=255)
                                                                      gaagaaGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCa                                                                            <  1:934486/149‑1 (MQ=255)
                                                                      gaagaaGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTCCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGAGCACGTTATggcgagccgccatgc                                                                            >  2:1176215/1‑136 (MQ=255)
                                                                       aagcagCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAg                                                                           <  1:164613/149‑1 (MQ=255)
                                                                          aagCCGATGATCGAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCCAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGaa                                                                        <  2:648712/149‑1 (MQ=255)
                                                                                atgacCCAGAACCAGCGCGCGCGTTTACCCGAGGTTTCGTTCAGTTTGAACCCGAACGCCTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCt                                                                  <  2:583416/149‑1 (MQ=255)
                                                                                   ttCCAGAACCAGACCGCGCGGTTGCCCCAGGTTTCTTTCAGTTCGAAACCGAACGCTTTAGGCTTCCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGtt                                                               <  2:855023/148‑1 (MQ=255)
                                                                                   atCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGCACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGACGATGATCAAGTTATGGAAGTGCGCAATTAGGAACAGGCTGtt                                                               <  2:298083/149‑1 (MQ=255)
                                                                                     ccagaaccagagCGCGCATTTACCCCAGGTTTCGTTCAGCTTGAAACCGAATGCTTTAGTCTCCCAGTAGGAAATCCCTCCGAAGCAGCCGGAGACCACCACGCAGAGGAACAAGTCATGGAAAGGCGCCTAAAGGAACAgccggctat                                                             >  2:56763/1‑141 (MQ=255)
                                                                                     ccagaaccagaaCGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTAt                                                             >  1:419592/1‑149 (MQ=255)
                                                                                         aaccagatcGGTCGGTTACCCCAGGTTTCGTGTAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAg                                                         <  2:257338/149‑1 (MQ=255)
                                                                                             agaacGCGCGTTTCCCGCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCTGGAACAGGCTGTTATGCAGAACg                                                     <  1:812904/149‑1 (MQ=255)
                                                                                                 cgcgtgTTTACCCCGTGTTTCGTTCCGTTTGAAACCGAACGCTTGAGGCTACCGGTAGGTCATCCCTGCGAAGCAGCCTAAGACCAGGGCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCGGTTATGCAGAACGAAGt                                                 <  2:15640/149‑1 (MQ=255)
                                                                                                   cgcgTTTACCCCAGGTTTCGGTCAGTTCGAAACCGAACGCTTTCAGGTACCAGTAGGTCATCGCAGAGCGTCTACCGAAGAAGAACTCTTGTATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTcc                                               <  1:293645/149‑1 (MQ=255)
                                                                                                         tACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGccc                                         <  1:1175340/149‑1 (MQ=255)
                                                                                                          aCCCCTGGTTTCGTTCAGTTTGAATCCGAACGCTTTAGGCTACCAGTAGGTCATCCGTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCg                                        <  2:117673/149‑1 (MQ=255)
                                                                                                            cccAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGACGTCCGCGCCCGGt                                      >  1:351917/1‑149 (MQ=255)
                                                                                                            cccAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGt                                      >  1:125719/1‑149 (MQ=255)
                                                                                                             ccAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCACTCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTa                                     >  1:56965/1‑149 (MQ=255)
                                                                                                                 gTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTAGGGAAGCGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGCCCGCGCCCGGCACGGc                                 >  2:57760/1‑149 (MQ=255)
                                                                                                                  tttCGTTCAGTTTGAAACCGAACGCATTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGCACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGcc                                >  1:352992/1‑149 (MQ=255)
                                                                                                                     cGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagc                             <  1:1037895/149‑1 (MQ=255)
                                                                                                                     cGTTCAGTTGGAAACCGAACGCTTTAGGCGACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagc                             <  1:81355/149‑1 (MQ=255)
                                                                                                                      gtttatttttactccGAATGCTCTCGGCCACCAGTATTTCATTCCGTCTAAGCATCCGAACACCTCGCCGCCGATGATCATGTTATGGAAGTGCGCAACCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagca                            <  2:653810/136‑1 (MQ=255)
                                                                                                                           gTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATGAGGAACAGGCTGTTATGCAGACCGAAGTCCGCGCCCGGAACGCCCGGCGGcgcg                       >  2:949472/1‑149 (MQ=255)
                                                                                                                              tGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCa                    <  1:226338/149‑1 (MQ=255)
                                                                                                                                 aaCCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGAACACGCCACTCCGGATCACTTAATGGATTACCCCAATCAGGAACAGGCTGTTATGCAGAACAAACACCGCCCCCGGTAAGGCCGCCGGCCCGCCAgcc                 >  2:140922/1‑147 (MQ=255)
                                                                                                                                     gAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATcc             <  1:1039364/149‑1 (MQ=255)
                                                                                                                                      aaCGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATccc            >  1:135504/1‑149 (MQ=255)
                                                                                                                                      aaCGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATccc            >  1:579375/1‑149 (MQ=255)
                                                                                                                                      aaCGCTTTAGGCCCCACGGGGCTCGTCCCTGCGAAGCAGCCGAAGACCCCGCCGCCGTTGTTGACGTTATGGAAGTGCGCAATCCGGAACAGGCTGTTATGCAGAACGAGGTCCGCGCCCGGTACGGCCACCAGCATGCCAGTCATccc            <  2:455035/149‑1 (MQ=255)
                                                                                                                                      aaCGCTTGAGGCCACCATTATGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATTCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATccc            <  1:1041984/149‑1 (MQ=255)
                                                                                                                                       aCGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCg           <  1:541703/149‑1 (MQ=255)
                                                                                                                                        cgcgTTAGGCCACCAGTAGGTCACCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGc          <  2:188600/149‑1 (MQ=255)
                                                                                                                                        cGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCGGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGc          <  2:209308/149‑1 (MQ=255)
                                                                                                                                             tAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCAcc     >  2:77159/1‑149 (MQ=255)
                                                                                                                                                gCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAGGTGCGCAATCAGGAACAGGCCGGTCTGCAGAACGAAGTCCCGGCCCGGGACAGCCAGCAGCACGCAAGCGATCCCGCCCCCCgcg  >  2:587699/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCGAG  >  CP009273/444633‑444925

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: