Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 448,966 T→C 5.4% H559R (CAC→CGC) ‡ cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913448,9660TC 5.4% 358.7 / 5.5 107H559R (CAC→CGC) ‡cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base T (50/51);  new base C (2/4);  total (52/55)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

CCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAG  >  NC_000913/448818‑449114
                                                                                                                                                    |                                                                                                                                                    
ccAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGt                                                                                                                                                      <  2:810244/149‑1 (MQ=255)
 cAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTg                                                                                                                                                     <  1:302304/149‑1 (MQ=255)
       aTGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTg                                                                                                                                               >  3:414532/1‑149 (MQ=255)
        tGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGc                                                                                                                                              <  1:505151/149‑1 (MQ=255)
         ggCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTgat                                                                                                                                             <  3:717648/149‑3 (MQ=255)
         ggCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTgat                                                                                                                                             <  4:314552/149‑3 (MQ=255)
         ggCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTgat                                                                                                                                             <  4:273741/149‑3 (MQ=255)
                ccGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTAc                                                                                                                                      >  4:704865/1‑149 (MQ=255)
                 cGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACg                                                                                                                                     <  4:414532/149‑1 (MQ=255)
                  gAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACgg                                                                                                                                    >  1:718223/1‑149 (MQ=255)
                     gATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCaa                                                                                                                                 >  4:723666/1‑149 (MQ=255)
                      aTGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGGATTTCTTCGTAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCAACAAGTTCGTGAACGTGCGGCACTACGGCaaa                                                                                                                                <  4:567716/149‑1 (MQ=255)
                      aTGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCaaa                                                                                                                                <  4:702428/149‑1 (MQ=255)
                      aTGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCaaa                                                                                                                                <  3:680707/149‑1 (MQ=255)
                       tggtggAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAg                                                                                                                               <  1:899854/149‑1 (MQ=255)
                              gaAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGt                                                                                                                                                      <  2:246708/119‑1 (MQ=255)
                              gaAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGt                                                                                                                                                      >  1:246708/1‑119 (MQ=255)
                                   cTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGa                                                                                                                   <  1:465132/149‑1 (MQ=255)
                                      cAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGc                                                                                                                <  4:114545/149‑1 (MQ=255)
                                       aaTGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCCTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCg                                                                                                               >  2:768285/1‑149 (MQ=255)
                                        aTGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCgg                                                                                                              >  1:487753/1‑149 (MQ=255)
                                        aTGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCgg                                                                                                              <  1:921832/149‑1 (MQ=255)
                                         tGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCggg                                                                                                             >  4:104167/1‑149 (MQ=255)
                                           aCGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGa                                                                                                           <  3:165170/149‑1 (MQ=255)
                                            cGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGaa                                                                                                          <  3:780472/149‑1 (MQ=255)
                                            cGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGaa                                                                                                          <  3:455983/149‑1 (MQ=255)
                                            cGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGaa                                                                                                          <  2:487753/149‑1 (MQ=255)
                                               tACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAagag                                                                                                       <  1:514092/149‑1 (MQ=255)
                                               tACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAagag                                                                                                       <  4:626236/149‑1 (MQ=255)
                                                 ccTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCCTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGt                                                                                                     <  1:768285/149‑1 (MQ=255)
                                                  cTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGtt                                                                                                    <  1:115989/149‑1 (MQ=255)
                                                    gCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGc                                                                                                  <  2:290183/149‑1 (MQ=255)
                                                    gCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGc                                                                                                  <  3:828099/149‑1 (MQ=255)
                                                    gCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTgat                                                                                                                                             >  2:624593/1‑104 (MQ=255)
                                                     cACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGcc                                                                                                 <  3:519152/149‑1 (MQ=255)
                                                          cTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTc                                                                                            >  2:488374/1‑149 (MQ=255)
                                                                 tCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTa                                                                                     <  3:388370/149‑1 (MQ=255)
                                                                   ggCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACg                                                                                   <  1:380386/149‑1 (MQ=255)
                                                                   ggCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCATGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACg                                                                                   >  4:735159/1‑149 (MQ=255)
                                                                      atatGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCCCGTTCATGAACGTGCGGCACTACGGCAAAGTtata                                                                                                                          <  3:551024/107‑1 (MQ=255)
                                                                      atatGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTtata                                                                                                                          >  4:551024/1‑107 (MQ=255)
                                                                       tatGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCa                                                                               <  3:853504/149‑1 (MQ=255)
                                                                        atGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCAc                                                                              <  4:220444/149‑1 (MQ=255)
                                                                         tGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCAcc                                                                             >  2:381818/1‑149 (MQ=255)
                                                                             tttCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCAc                                                                         >  3:848980/1‑149 (MQ=255)
                                                                             tttCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCAc                                                                         >  4:723159/1‑149 (MQ=255)
                                                                                 ttcATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGt                                                                     >  2:651772/1‑149 (MQ=255)
                                                                                 ttcATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGt                                                                     >  1:924154/1‑149 (MQ=255)
                                                                                  tcATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTc                                                                    <  4:483754/149‑1 (MQ=255)
                                                                                  tcATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTc                                                                    <  4:696359/149‑1 (MQ=255)
                                                                                  tcATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTc                                                                    >  2:363536/1‑149 (MQ=255)
                                                                                   cATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCg                                                                   >  3:333672/1‑149 (MQ=255)
                                                                                    aTAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGc                                                                  >  2:570717/1‑149 (MQ=255)
                                                                                         ggTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTc                                                             >  2:269081/1‑149 (MQ=255)
                                                                                             aGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGt                                                         <  1:363536/149‑1 (MQ=255)
                                                                                               gCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCa                                                       >  1:748424/1‑149 (MQ=255)
                                                                                               gCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCa                                                       >  2:739118/1‑149 (MQ=255)
                                                                                                 tttttATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACg                                                     >  4:211493/1‑149 (MQ=255)
                                                                                                 tttttATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACg                                                     >  1:525436/1‑149 (MQ=255)
                                                                                                        aCGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTgg                                              >  4:514382/1‑149 (MQ=255)
                                                                                                        aCGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTgg                                              >  2:25007/1‑149 (MQ=255)
                                                                                                          gCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCggg                                                                                                             <  1:739480/84‑1 (MQ=255)
                                                                                                          gCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCggg                                                                                                             >  2:739480/1‑84 (MQ=255)
                                                                                                          gCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtc                                            <  1:101338/149‑1 (MQ=255)
                                                                                                          gCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtc                                            >  1:141385/1‑149 (MQ=255)
                                                                                                            ttCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtcgc                                          <  1:281481/149‑1 (MQ=255)
                                                                                                            ttCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtcgc                                          >  1:239362/1‑149 (MQ=255)
                                                                                                            ttCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtcgc                                          >  3:84186/1‑149 (MQ=255)
                                                                                                               gCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGGTCTggtcgcggt                                       >  1:286061/1‑149 (MQ=255)
                                                                                                                  tttCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtcgcggtcgc                                    <  3:735300/149‑1 (MQ=255)
                                                                                                                     ctctTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGaa                                 <  4:866075/149‑1 (MQ=255)
                                                                                                                     ctctTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGaa                                 <  4:796459/149‑1 (MQ=255)
                                                                                                                       ctTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGCCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATa                               <  2:791406/149‑1 (MQ=255)
                                                                                                                        tttCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTgg                                              <  1:327453/133‑1 (MQ=255)
                                                                                                                        tttCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTgg                                              >  2:327453/1‑133 (MQ=255)
                                                                                                                        tttCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAg                              >  4:100145/1‑149 (MQ=255)
                                                                                                                           cATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGaaa                           >  4:590678/1‑149 (MQ=255)
                                                                                                                            aTTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCCGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAAc                          >  4:365982/1‑149 (MQ=255)
                                                                                                                            aTTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAAc                          <  3:699226/149‑1 (MQ=255)
                                                                                                                            aTTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAAc                          >  4:623893/1‑149 (MQ=255)
                                                                                                                              ttCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGt                        <  3:100145/149‑1 (MQ=255)
                                                                                                                               tCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTa                       <  1:570717/149‑1 (MQ=255)
                                                                                                                                 ccAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTAc                                                                                    >  1:635910/1‑86 (MQ=255)
                                                                                                                                 ccAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTAc                                                                                    <  2:635910/86‑1 (MQ=255)
                                                                                                                                 ccAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACa                     >  2:363332/1‑149 (MQ=255)
                                                                                                                                     aaTGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTg                 >  1:171408/1‑149 (MQ=255)
                                                                                                                                         ccgcACGTTCACGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATa             >  4:683931/4‑149 (MQ=255)
                                                                                                                                          aTCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATaa            >  2:566868/1‑149 (MQ=255)
                                                                                                                                          aTCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATaa            <  3:253295/149‑1 (MQ=255)
                                                                                                                                          aTCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGCTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCt                                                >  1:802605/1‑113 (MQ=255)
                                                                                                                                          aTCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAAGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCt                                                <  2:802605/113‑1 (MQ=255)
                                                                                                                                           tCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACTGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAAc           >  4:227279/1‑149 (MQ=255)
                                                                                                                                            cACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACg          >  2:5535/1‑149 (MQ=255)
                                                                                                                                             aCGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGa         <  4:334633/149‑1 (MQ=255)
                                                                                                                                             aCGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGa         >  2:623995/1‑149 (MQ=255)
                                                                                                                                              cGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAg        <  2:396376/149‑1 (MQ=255)
                                                                                                                                               gTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAgg       >  3:891648/1‑149 (MQ=255)
                                                                                                                                                 tCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa     <  3:365982/149‑1 (MQ=255)
                                                                                                                                                   gTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCaga   >  1:383997/1‑149 (MQ=255)
                                                                                                                                                    tGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCagag  >  3:766881/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAG  >  NC_000913/448818‑449114

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: