Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 449,925 A→G I239I (ATT→ATC)  cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913449,9250AG99.2% 523.0 / NA 134I239I (ATT→ATC) cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (68/65);  total (68/66)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

GTAAACTTCCGGGTGGCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATC  >  NC_000913/449778‑450072
                                                                                                                                                     |                                                                                                                                                   
gTAAACTTCCGGGTGGCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATc                                                                                                                                                  <  2:402564/151‑1 (MQ=255)
 tAAACTTCCGGGTGGCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATc                                                                                                                                                  >  2:515689/1‑150 (MQ=255)
  aaaCTTCCGGGTGGCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATc                                                                                                                                                  >  1:11006/1‑149 (MQ=255)
  aaaCTTCCGGGTGGCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATc                                                                                                                                                  >  4:242311/1‑149 (MQ=255)
    aCTTCCGGGTGGCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTa                                                                                                                                              <  2:293470/151‑1 (MQ=255)
    aCTTCCGGGTGGCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTa                                                                                                                                              <  3:89440/151‑1 (MQ=255)
       tCCGGGTGGCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACg                                                                                                                                            <  3:92409/150‑1 (MQ=255)
         cGGGTGGCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTc                                                                                                                                         >  1:1121251/1‑151 (MQ=255)
           ggaggCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGAGCCCAGATAGCGAGCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTcgc                                                                                                                                       <  1:45756/148‑1 (MQ=255)
              ggCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCgcgc                                                                                                                                     >  1:461572/1‑150 (MQ=255)
              ggCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCgcgc                                                                                                                                     <  1:1364562/150‑1 (MQ=255)
               gCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGcac                                                                                                                                   <  1:1708042/151‑1 (MQ=255)
                ccccAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGcaca                                                                                                                                  <  2:1091194/151‑1 (MQ=255)
                  ccAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGt                                                                                                                                <  2:1547736/151‑1 (MQ=255)
                  ccAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGt                                                                                                                                <  2:353380/151‑1 (MQ=255)
                     ggCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGAt                                                                                                                             <  1:1119000/151‑1 (MQ=255)
                      gCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATg                                                                                                                            <  2:1641006/151‑1 (MQ=255)
                      gCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATg                                                                                                                            >  1:563895/1‑151 (MQ=255)
                       cccAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGc                                                                                                                           <  3:128767/151‑1 (MQ=255)
                        ccAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGcc                                                                                                                          >  4:346532/1‑151 (MQ=255)
                        ccAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGcc                                                                                                                          >  2:376392/1‑151 (MQ=255)
                         cAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGccc                                                                                                                         >  2:1303764/1‑151 (MQ=255)
                          aaaTCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGccc                                                                                                                         >  1:1150625/1‑150 (MQ=255)
                          aaaTCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGccc                                                                                                                         >  1:255372/1‑150 (MQ=255)
                          aaaTCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGccc                                                                                                                         >  2:1640035/1‑150 (MQ=255)
                          aaaTCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGccc                                                                                                                         >  2:1186893/1‑150 (MQ=255)
                           aaTCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAg                                                                                                                       <  2:1234402/151‑1 (MQ=255)
                              cAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTa                                                                                                                    <  3:324930/151‑1 (MQ=255)
                               aGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTa                                                                                                                    >  2:1477808/1‑150 (MQ=255)
                               aGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTa                                                                                                                    >  1:424564/1‑150 (MQ=255)
                                ggTTGTTGTACAACATCATGTTTCCACCCATATCGTTTGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTa                                                                                                                    >  1:1077445/1‑149 (MQ=255)
                                ggTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTaa                                                                                                                   >  2:772602/1‑150 (MQ=255)
                                   tGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATa                                                                                                                <  1:1048163/150‑1 (MQ=255)
                                     aTGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGTTCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTg                                                                                                             <  1:1312046/151‑1 (MQ=255)
                                        tACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGc                                                                                                           >  2:1624710/1‑150 (MQ=255)
                                         aCATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCAt                                                                                                         <  2:183002/151‑1 (MQ=255)
                                         aCATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCAt                                                                                                         <  2:1565727/151‑1 (MQ=255)
                                         aCATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCAt                                                                                                         <  2:335191/151‑1 (MQ=255)
                                         aCATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCAt                                                                                                         <  2:1436053/151‑1 (MQ=255)
                                           atcatcatGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCt                                                                                                       >  1:48005/1‑151 (MQ=255)
                                            tcatcatGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCtt                                                                                                      >  2:655675/1‑151 (MQ=255)
                                             catcatGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTg                                                                                                     <  1:1723519/151‑1 (MQ=255)
                                              atcatGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGa                                                                                                    <  4:563942/151‑1 (MQ=255)
                                               tcatGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGa                                                                                                    <  2:554973/150‑1 (MQ=255)
                                                   gTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATg                                                                                               <  2:255414/151‑1 (MQ=255)
                                                    ttGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATgg                                                                                              >  3:547584/1‑151 (MQ=255)
                                                     tGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATgg                                                                                              <  1:165864/150‑1 (MQ=255)
                                                       ccACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTc                                                                                            <  1:796712/150‑1 (MQ=255)
                                                       ccACCCAT‑TCCGTGGTAAAGGAATTGGTGCCCAGGAAGCCATCCAGGGTTAACAACGGGGACGTTACCGGTC‑GAATTGGGAAAGAAAGAAATAA‑CCGTACCGTCCCGGACCGGTATTGCCAGGTAAATACTGGCATCtt                                                                                                      <  1:516644/139‑1 (MQ=255)
                                                        cACCCATATCGTTGGTAAAGAAATGGGCGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCAt                                                                                          >  1:1417274/1‑151 (MQ=255)
                                                            cATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGccc                                                                                      <  3:489209/151‑1 (MQ=255)
                                                               atCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGt                                                                                   <  1:1673749/151‑1 (MQ=255)
                                                                tCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTg                                                                                  >  4:158053/1‑151 (MQ=255)
                                                                   ttGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTgcgc                                                                               >  2:1355651/1‑151 (MQ=255)
                                                                   ttGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATG‑CCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTgcgcg                                                                              >  3:257685/1‑151 (MQ=255)
                                                                    tGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTgcgcg                                                                              >  4:138279/1‑151 (MQ=255)
                                                                    tGGTAAAGAAATGGGTGCCAAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCACGGTgcg                                                                                >  1:51333/1‑149 (MQ=255)
                                                                       tAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCa                                                                            >  1:7211/1‑150 (MQ=255)
                                                                        aaaGAACTGGGGGCCCAGATAGCGATCCAGGGTCAACAAC‑GCCACGGTAACCGACAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATCCTGGCATCTTGAACCTGGTCATGCCCGGAGCGCGCATc                                                                          >  4:343152/1‑151 (MQ=255)
                                                                           gAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTc                                                                       >  1:215922/1‑151 (MQ=255)
                                                                            aaaTGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTc                                                                       >  2:1363389/1‑150 (MQ=255)
                                                                            aaaTGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCa                                                                      <  2:306140/151‑1 (MQ=255)
                                                                            aaaTGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCa                                                                      <  2:1545413/151‑1 (MQ=255)
                                                                              aTGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCGGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGa                                                                    <  1:312706/151‑1 (MQ=255)
                                                                               tGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGa                                                                    >  1:1639767/1‑150 (MQ=255)
                                                                                gggTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAAt                                                                  >  3:172011/1‑151 (MQ=255)
                                                                                 ggTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATg                                                                 <  2:1668156/151‑1 (MQ=255)
                                                                                   tGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGt                                                               >  1:776309/1‑151 (MQ=255)
                                                                                     cccAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTaa                                                             <  4:488391/151‑1 (MQ=255)
                                                                                     cccAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTa                                                              >  2:907683/1‑150 (MQ=255)
                                                                                       cAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACg                                                           <  1:1101292/151‑1 (MQ=255)
                                                                                           tAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTAAGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTaatg                                                           >  1:367782/1‑145 (MQ=255)
                                                                                            aGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAGTGGTAACgaagaa                                                      <  1:256208/151‑1 (MQ=255)
                                                                                            aGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACgaaga                                                       >  1:344657/1‑150 (MQ=255)
                                                                                              cGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGt                                                    <  1:1477649/151‑1 (MQ=255)
                                                                                               gATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGtt                                                   >  2:772628/1‑151 (MQ=255)
                                                                                                 tCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTg                                                  >  3:367966/1‑150 (MQ=255)
                                                                                                  ccAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCTCGCACAGTGATGCCCATGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCCCCTTCAGAATAGTAACGTAGAAGTTGa                                                 <  1:177384/150‑1 (MQ=255)
                                                                                                  ccAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGAt                                                <  4:391666/151‑1 (MQ=255)
                                                                                                   cAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGAt                                                >  2:677499/1‑150 (MQ=255)
                                                                                                     gggTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATAcc                                             <  1:640195/151‑1 (MQ=255)
                                                                                                      ggTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCg                                            <  3:158074/151‑1 (MQ=255)
                                                                                                         caacaaC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGt                                          >  2:496857/1‑150 (MQ=255)
                                                                                                         caacaaC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGt                                          <  2:464524/150‑1 (MQ=255)
                                                                                                           acaac‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAg                                       <  2:1031975/151‑1 (MQ=255)
                                                                                                              ac‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGcgt                                    <  2:653151/151‑1 (MQ=255)
                                                                                                                 gcgACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGcgtcg                                  <  1:376340/151‑1 (MQ=255)
                                                                                                                         aaCCGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCtata                          <  3:18855/151‑1 (MQ=255)
                                                                                                                           ccGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATAcc                        >  3:534355/1‑151 (MQ=255)
                                                                                                                            cGTCAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCg                       >  1:160355/1‑151 (MQ=255)
                                                                                                                               cAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGAt                    <  2:1325649/151‑1 (MQ=255)
                                                                                                                               cAGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCAAGGTAAAAACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGAt                    >  4:111091/1‑151 (MQ=255)
                                                                                                                                aGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGAt                    >  1:1143191/1‑150 (MQ=255)
                                                                                                                                aGAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATa                   <  2:272094/151‑1 (MQ=255)
                                                                                                                                 gAATTGGG‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGAt                    >  2:1459275/1‑149 (MQ=255)
                                                                                                                                      ggg‑AAGGAAGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTgg              <  2:7212/150‑1 (MQ=255)
                                                                                                                                       gg‑aaggaaGCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGgag            <  2:1060565/151‑1 (MQ=255)
                                                                                                                                        g‑aaggaagCACTGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGCACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTgg              >  1:708532/1‑148 (MQ=255)
                                                                                                                                          aaggaagCAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGAc          <  2:1077575/151‑1 (MQ=255)
                                                                                                                                                 cAATGATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCata    <  1:907570/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTAAACTTCCGGGTGGCCCCAGGCCCAAATCAGGTTGATGTACATCATCATGTTGCCACCCATATCGTTGGTAAAGAAATGGGTGCCCAGATAGCGATCCAGGGTCAACAAC‑GCGACGGTAACCGTCAGAATTGGG‑AAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATC  >  NC_000913/449778‑450072

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑