Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 444,779 C→T 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,7790CT98.4% 196.6 / NA 65W440* (TGG→TAG) cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (27/35);  minor base G (0/1);  total (29/36)
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

AATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCG  >  CP009273/444637‑444923
                                                                                                                                              |                                                                                                                                                
aacctgctgGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGt                                                                                                                                            <  1:460362/146‑1 (MQ=255)
  tcggctgGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCAGACGGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAGCCGAACGCTTTAGGCTACCAGTAg                                                                                                                                          <  2:1218797/146‑1 (MQ=255)
        tgCTCAAACGACCGGTCATGCCCATGAAGCCCAGCGCATACAGTTGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCACGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATc                                                                                                                                    <  1:18686/148‑1 (MQ=255)
         gCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCACGTTTACCCCAGATTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATcc                                                                                                                                   >  1:110034/1‑149 (MQ=255)
          cTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCAGACAGTGGCATAAAGGCAACGAAGAAGTCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATccc                                                                                                                                  <  2:961193/149‑1 (MQ=255)
            cAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTg                                                                                                                                <  1:1262105/149‑1 (MQ=255)
               acgacgGGTCATGGCCATCAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCACGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGa                                                                                                                             <  1:559780/149‑1 (MQ=255)
               acgacgGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGa                                                                                                                             <  2:264440/149‑1 (MQ=255)
                cgacgGGTCATGCCCCTGCAGCCCAGCGCATACAGTGGCATACAGGCAACGAAGAAGCCGATGATCCAGAACCAGCACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGaa                                                                                                                            <  1:1101934/149‑1 (MQ=255)
                   cgGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAagca                                                                                                                         <  2:125343/149‑1 (MQ=255)
                   cgGGTCATGCCCATGAAGCCCAGCGCATACAGGGGCATAAAGGCAATGAAGTTGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAagca                                                                                                                         <  1:179776/149‑1 (MQ=255)
                   ccggTCATGCCCATGAAGCCCAGCCCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAagca                                                                                                                         <  1:337061/147‑1 (MQ=255)
                     ggTCATGCCCATGAAGCCCAGCGCATACGGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAagcagc                                                                                                                       <  1:356178/149‑1 (MQ=255)
                         aTGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGCCATCCCAGAGAAACAGCCGaa                                                                                                                   >  2:1122144/1‑149 (MQ=255)
                              cATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCa                                                                                                              >  1:1285158/1‑149 (MQ=255)
                               aTGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCAc                                                                                                             >  1:205894/1‑149 (MQ=255)
                                  aggCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCATCTGAAGACCAcgcc                                                                                                          <  2:1146969/147‑1 (MQ=255)
                                    gCCCAGCGGATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCAcgccgc                                                                                                        <  2:905233/149‑1 (MQ=255)
                                      caaGCGCATACAGTGGCATAAAGGCAACAAAGTCAGAAATGATCCAGAACCAGAACGCGCGCTTACACCAGGTTTCGTTCAGGTTGAAACCGAAGGCTTTCGGCTACCACGAGGTGATCCCTGCGAAGCAGCCGAAGACCACGaccccg                                                                                                      >  2:617226/3‑149 (MQ=255)
                                       ctgcccATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCga                                                                                                     <  1:514341/144‑1 (MQ=255)
                                       cAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGCTCCAGAACCAGAACGCGCGTTTACCCAATGTTTCGTTCAGTTTGAAACCGAACGCCTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCga                                                                                                     <  1:514056/149‑1 (MQ=255)
                                         gcgcATAGAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCgagg                                                                                                   >  2:905585/1‑147 (MQ=255)
                                           gcATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGACCCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCgatgat                                                                                                 <  1:555773/149‑1 (MQ=255)
                                                   tGGCATAAAGGCAACGAAGAAGCCGATAATGCAGCACCAGAACGCGCGTTTACACCCGGTACCGTTGAGTCGGCAACCGGACGCTTCAGGCCACGAGAAGGTCATCCCGGCGAAGCAGCCGAAGACCACGCCGCCGAGGATGACGtgag                                                                                         >  2:362426/1‑146 (MQ=255)
                                                    ggCATAAAGCGAACGAAGAGGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCCGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTATCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATg                                                                                        <  2:861986/149‑1 (MQ=255)
                                                     gCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGAAGCCGAAGACCACGCCGCCGATGCTCACGTTATgg                                                                                       >  1:7235/1‑149 (MQ=255)
                                                       aTAAAGGCAACGGAGAAGCCGATGAGCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGaa                                                                                     <  2:153267/149‑1 (MQ=255)
                                                       aTAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGaa                                                                                     >  2:898164/1‑149 (MQ=255)
                                                         aaaGCCAACGACGACGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGt                                                                                   <  2:757809/149‑1 (MQ=255)
                                                           aGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTTTTTCAGTTTGAAAGCGAACGCTTTAGGCGACCAGTAGGTCCTCCCGGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTgc                                                                                 <  1:387961/149‑1 (MQ=255)
                                                               aaCGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCGGCCGAGGACCACGCCGCCGATGATCACGTGATGGAAGTGGGCac                                                                             >  2:444953/1‑148 (MQ=255)
                                                               aaCGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGAGGATCACCGATTGGAAGTgcgcca                                                                             >  1:848849/1‑147 (MQ=255)
                                                                  gaagaaGCCGATTATCCAGAACTAGATCGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCa                                                                          <  1:1217535/149‑1 (MQ=255)
                                                                  gaagaaGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCa                                                                          >  2:1099881/1‑149 (MQ=255)
                                                                   aagaagCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAg                                                                         <  2:1213532/149‑1 (MQ=255)
                                                                      aagCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGaa                                                                      >  1:95872/1‑149 (MQ=255)
                                                                              gatCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGCACAGGCTgg                                                              >  1:120511/1‑148 (MQ=255)
                                                                              gatCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCAAAGCAGTGGAAGGCCACGCCGCCGATGACCCCGATTCGGAAATGCGACATCAGGAAGAGGcgga                                                              >  2:529828/1‑146 (MQ=255)
                                                                                   agaaccagaacGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACAATTTAAGCAAAAACCATATCATACCTGCGTAGCAGCCGAAGACCACGACGCCGATGATCACGTGATGGAAGTGCGCAATCAAGAACCGGCGGTTATGc                                                         >  1:315148/1‑149 (MQ=255)
                                                                                      accagaacGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGa                                                      >  1:442906/1‑149 (MQ=255)
                                                                                      accacaacGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGa                                                      <  2:1193935/149‑1 (MQ=255)
                                                                                        cagaacGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAAc                                                    <  1:1099881/149‑1 (MQ=255)
                                                                                           aacCCGCGTTTACTCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGaa                                                 <  1:405015/149‑1 (MQ=255)
                                                                                                gcgTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCcg                                            >  1:557144/1‑149 (MQ=255)
                                                                                                   tttACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCgcgc                                         >  1:1137694/1‑149 (MQ=255)
                                                                                                   tttACCCCAGGGTGCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCTCGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCgcgc                                         <  1:1285688/149‑1 (MQ=255)
                                                                                                     tACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGccc                                       <  2:686519/149‑1 (MQ=255)
                                                                                                     gaCCCCAGGGTTCGTTCCGTTTGAAACCGAACGCTTTAGGCTACCAGTCGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGACCACGTTATGGAAGTGCGCGATCAGGAACTGGCTGTTATGCAGAACGAAGTCCGCGccc                                       <  2:12388/148‑1 (MQ=255)
                                                                                                         ccAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTa                                   >  2:1170345/1‑149 (MQ=255)
                                                                                                          cctaaagcGTTCGGTTTCAAACTGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCggagc                                  >  1:443373/8‑146 (MQ=255)
                                                                                                                 cGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagc                           <  1:1057657/149‑1 (MQ=255)
                                                                                                                  gTTCAGTTTGAAGCCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAAAAGGCGGTTATGCAGAACGAGGTCGGCCCCCGGTACGGgcagca                          >  2:749352/1‑149 (MQ=255)
                                                                                                                  gTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcc                          >  1:533560/1‑148 (MQ=255)
                                                                                                                   ttCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag                         >  2:1279561/1‑149 (MQ=255)
                                                                                                                     gatgtggAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAgcagca                       <  1:956038/143‑1 (MQ=255)
                                                                                                                     cAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAgcagca                       <  2:940698/149‑1 (MQ=255)
                                                                                                                       gTTTGAAACCGAACGCTTTCGGCTACCGGTACGTCACCACGGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACg                     <  1:1024506/149‑1 (MQ=255)
                                                                                                                        tttGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGc                    >  1:376074/1‑149 (MQ=255)
                                                                                                                        tttGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGc                    <  1:540492/149‑1 (MQ=255)
                                                                                                                        tttGAAACCGAACGCGTGCGGCTTCCAGTGGGGCATCCCGGCCAAGCAGCCGAAGCCCACGCCGCCGATGATCACGTTCTGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGACGTCCGCGCCCGGTACGGCCAGCAGCACGc                    <  1:806037/149‑1 (MQ=255)
                                                                                                                           gAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAg                 >  1:957426/1‑149 (MQ=255)
                                                                                                                              aCCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGGCa              >  1:1170186/1‑149 (MQ=255)
                                                                                                                               ccGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTAAGGAAGTGCGCCATCAGGAACAGGCGGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAgacca             >  2:987890/1‑145 (MQ=255)
                                                                                                                                  aaCGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATccc          >  2:395622/1‑149 (MQ=255)
                                                                                                                                          aGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGTGATCAGGTATAGGCGGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCg  <  2:805518/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
AATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCG  >  CP009273/444637‑444923

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: ‑