Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 450,780 T→C T19A (ACT→GCT)  cyoA ← cytochrome o ubiquinol oxidase subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913450,7800TC98.0% 623.2 / 0.1 198T19A (ACT→GCT) cyoAcytochrome o ubiquinol oxidase subunit II
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (96/98);  total (97/101)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 202 total aligned reads displayed.

CAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTT  >  NC_000913/450690‑450869
                                                                                          |                                                                                         
cAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAgcgc                                                                                         >  1:246520/1‑93 (MQ=255)
 aaaTGCCGTCAGTATCAGTGAGCGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGcc                                                                                        <  1:970981/93‑1 (MQ=255)
 aaaTGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGcc                                                                                        >  2:2930730/1‑93 (MQ=255)
 aaaTGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGcc                                                                                        >  2:264940/1‑93 (MQ=255)
 aaaTGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGcc                                                                                        <  1:146569/93‑1 (MQ=255)
 aaaTGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGcc                                                                                        <  1:933399/93‑1 (MQ=255)
 aaaTGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGcc                                                                                        >  1:906571/1‑93 (MQ=255)
 aaaTGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGcc                                                                                        <  1:4290941/93‑1 (MQ=255)
  aaTGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCt                                                                                       >  1:107674/1‑93 (MQ=255)
  aaTGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCt                                                                                       <  2:2585665/93‑1 (MQ=255)
    tGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGc                                                                                     <  1:4509758/93‑1 (MQ=255)
    tGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGc                                                                                     >  1:651285/1‑93 (MQ=255)
      ccGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCaa                                                                                   <  1:190480/93‑1 (MQ=255)
      ccGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCaa                                                                                   >  2:3029085/1‑93 (MQ=255)
       cGTCGGTATCAGTGAGCGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCgaa                                                                                  <  2:4218589/93‑1 (MQ=255)
        gTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAaat                                                                                 >  2:2340482/1‑93 (MQ=255)
            gTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATg                                                                             >  1:4215602/1‑93 (MQ=255)
              aTCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGAc                                                                           <  1:1758184/93‑1 (MQ=255)
              aTCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGAc                                                                           <  1:889188/93‑1 (MQ=255)
               tCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACa                                                                          >  2:3900812/1‑93 (MQ=255)
               tCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACa                                                                          <  1:3902184/93‑1 (MQ=255)
               tCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACa                                                                          <  1:1419142/93‑1 (MQ=255)
               tCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACa                                                                          <  2:2622517/93‑1 (MQ=255)
                 aGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAAc                                                                        >  2:4378938/1‑93 (MQ=255)
                  gTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAAcc                                                                       >  1:1306162/1‑93 (MQ=255)
                  gTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAAcc                                                                       >  1:1675508/1‑93 (MQ=255)
                  gTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAAcc                                                                       <  1:2111776/93‑1 (MQ=255)
                    gAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCAt                                                                     <  2:395731/93‑1 (MQ=255)
                       cGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATccc                                                                  <  1:3785336/93‑1 (MQ=255)
                         ttGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCaa                                                                >  1:1296598/1‑93 (MQ=255)
                         ttGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCaa                                                                >  1:4392212/1‑93 (MQ=255)
                            cTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACt                                                             >  2:1927826/1‑93 (MQ=255)
                            cTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACt                                                             <  2:1385974/93‑1 (MQ=255)
                            cTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACt                                                             >  2:548007/1‑93 (MQ=255)
                            cTCAAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACt                                                             <  2:3735214/93‑1 (MQ=255)
                             tCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACtt                                                            >  1:208417/1‑93 (MQ=255)
                             cccAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACtt                                                            <  1:3236163/92‑1 (MQ=255)
                              ccAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACttt                                                           <  1:328119/93‑1 (MQ=255)
                               cAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACtttt                                                          <  1:500669/93‑1 (MQ=255)
                                aGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTa                                                         <  2:209479/93‑1 (MQ=255)
                                aGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTa                                                         >  1:3001786/1‑93 (MQ=255)
                                aGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTa                                                         >  1:2685202/1‑93 (MQ=255)
                                 gACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTAt                                                        >  2:2108137/1‑93 (MQ=255)
                                 gACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTAt                                                        >  1:3341575/1‑93 (MQ=255)
                                  aCCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTAtt                                                       <  2:1332191/93‑1 (MQ=255)
                                  aCCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTAtt                                                       <  1:2754866/93‑1 (MQ=255)
                                   ccAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTg                                                      <  1:2043383/93‑1 (MQ=255)
                                   ccAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTg                                                      >  1:222462/1‑93 (MQ=255)
                                   ccAAGCTGTTCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTg                                                      <  2:1191031/93‑1 (MQ=255)
                                    cAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGt                                                     <  1:759454/93‑1 (MQ=255)
                                    cAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGt                                                     >  1:872917/1‑93 (MQ=255)
                                     aaTCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCACACTTTTATTGtt                                                    >  1:3604485/1‑92 (MQ=255)
                                       tCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTAtt                                                  <  1:2581865/93‑1 (MQ=255)
                                        cTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTAttt                                                 <  1:2827093/93‑1 (MQ=255)
                                          gTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTcc                                               <  1:1539228/93‑1 (MQ=255)
                                           tCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCt                                              >  2:2081027/1‑93 (MQ=255)
                                            ccTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTg                                             >  1:436864/1‑93 (MQ=255)
                                              tttGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAg                                           >  2:3245481/1‑93 (MQ=255)
                                                tGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGtc                                         <  2:2176286/93‑1 (MQ=255)
                                                  ggATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGtctc                                       <  1:2333972/93‑1 (MQ=255)
                                                  ggATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGtctc                                       <  2:3473529/93‑1 (MQ=255)
                                                  ggATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGtctc                                       <  1:3668685/93‑1 (MQ=255)
                                                     tCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCAtt                                    <  1:3854723/93‑1 (MQ=255)
                                                       tAACAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTa                                  <  2:2245070/93‑1 (MQ=255)
                                                         aCAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAAc                                <  1:4458994/93‑1 (MQ=255)
                                                          cAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACg                               <  2:3403259/93‑1 (MQ=255)
                                                          cAGCGCAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACg                               >  1:2152506/1‑93 (MQ=255)
                                                             cgcAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGAcc                            <  2:3378598/93‑1 (MQ=255)
                                                             cgcAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGAcc                            >  1:3164859/1‑93 (MQ=255)
                                                             cgcAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGAcc                            >  2:4690992/1‑93 (MQ=255)
                                                              gcAGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCt                           <  2:1785911/93‑1 (MQ=255)
                                                                aGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCa                         >  1:3386192/1‑93 (MQ=255)
                                                                aGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCa                         >  1:3971023/1‑93 (MQ=255)
                                                                aGAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCa                         >  2:3312850/1‑93 (MQ=255)
                                                                 gAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCaa                        >  1:3874948/1‑93 (MQ=255)
                                                                 gAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCaa                        <  2:2480612/93‑1 (MQ=255)
                                                                 gAATTACAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCaa                        >  2:4591734/1‑93 (MQ=255)
                                                                      aCAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCa                   >  1:677946/1‑93 (MQ=255)
                                                                      aCAGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCa                   <  1:1943800/93‑1 (MQ=255)
                                                                        aGCCACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACg                 <  1:3490840/93‑1 (MQ=255)
                                                                          ccACTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACggg               <  2:4299736/93‑1 (MQ=255)
                                                                            aCTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGAt             >  2:3164662/1‑93 (MQ=255)
                                                                            aCTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGAt             >  2:4124092/1‑93 (MQ=255)
                                                                            aCTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGAt             <  2:2137301/93‑1 (MQ=255)
                                                                            aCTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGAt             <  2:1087015/93‑1 (MQ=255)
                                                                             cTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATc            >  2:1713877/1‑93 (MQ=255)
                                                                             cTGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGAATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATc            >  1:2626363/1‑93 (MQ=255)
                                                                              tGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCt           >  2:3942396/1‑93 (MQ=255)
                                                                              tGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCt           >  1:1125871/1‑93 (MQ=255)
                                                                              tGAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCt           <  1:460634/93‑1 (MQ=255)
                                                                               gAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCtg          <  1:2015548/93‑1 (MQ=255)
                                                                               gAGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCtg          <  1:184925/93‑1 (MQ=255)
                                                                                aGCAATACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCtgg         <  2:3526864/93‑1 (MQ=255)
                                                                                   aaTACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCtggtgg      <  2:1908747/93‑1 (MQ=255)
                                                                                    aTACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGc     >  2:3715677/1‑93 (MQ=255)
                                                                                    aTACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGc     >  2:2194311/1‑93 (MQ=255)
                                                                                    aTACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGc     >  2:747567/1‑93 (MQ=255)
                                                                                    aTACAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGc     <  2:859502/93‑1 (MQ=255)
                                                                                       cAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGtt  <  2:1946745/93‑1 (MQ=255)
                                                                                       cAGCGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGtt  <  2:1225057/93‑1 (MQ=255)
                                                                                          |                                                                                         
CAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTT  >  NC_000913/450690‑450869

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