Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 450,001 T→C 6.3% K214R (AAG→AGG)  cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913450,0010TC 6.3% 192.6 / 3.3 64K214R (AAG→AGG) cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base T (26/33);  new base C (2/2);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GGTGCCCAGATAGCGATCCAGGGTCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGC  >  NC_000913/449859‑450146
                                                                                                                                              |                                                                                                                                                 
ggTGCCCAGATAGCGATCCAGGGTGAACAACGCGACGGTTACCGTCAGAATTGGGAAGGAACCAAAAATCAGTACGATCGCGCAGAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGa                                                                                                                                              <  7:223646/148‑1 (MQ=255)
ggTGCCCAGATAGCGATCCAGGGTCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGa                                                                                                                                              <  1:141132/148‑1 (MQ=255)
ggTGCCCAGATAGCGATCCAGGGTCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGa                                                                                                                                              <  8:85708/148‑1 (MQ=255)
      cAGATAGCGATCCAGGGTCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGt                                                                                                                                         >  4:319416/1‑147 (MQ=255)
          tAGCGATCCAGGGTCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACga                                                                                                                                    <  3:119263/148‑1 (MQ=255)
                tCCAGGGTCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGt                                                                                                                              <  3:146570/148‑1 (MQ=255)
                tCCAGGGTCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGt                                                                                                                              <  6:153348/148‑1 (MQ=255)
                tCCAGGGTCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGt                                                                                                                              <  7:88274/148‑1 (MQ=255)
                tCCAGGGTCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGt                                                                                                                              <  2:117482/148‑1 (MQ=255)
                  cAGGGTCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTg                                                                                                                            <  3:149888/148‑1 (MQ=255)
                    gggTCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGAt                                                                                                                          <  1:326991/148‑1 (MQ=255)
                     ggTCAACAACGCGTCGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGGCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAACAAGTTGATa                                                                                                                         <  3:200575/148‑1 (MQ=255)
                     ggTCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGATTGGTAACGAAGAAGTTGATa                                                                                                                         <  2:167053/148‑1 (MQ=255)
                       tCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATAcc                                                                                                                       >  1:284221/1‑148 (MQ=255)
                        caacaaCGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCg                                                                                                                      >  8:313708/1‑148 (MQ=255)
                           caacGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCATTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGt                                                                                                                    >  1:362377/1‑147 (MQ=255)
                           caacGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGt                                                                                                                    >  1:41792/1‑147 (MQ=255)
                           caacGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGt                                                                                                                    >  8:115327/1‑147 (MQ=255)
                            aacGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGt                                                                                                                    >  4:75758/1‑146 (MQ=255)
                                cgACGGGAACCGTCAGAATTGGGAAGGAAGCAATACTCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGcgt                                                                                                              <  6:281264/148‑1 (MQ=255)
                                      tAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCAGGTGCGCGCATCTTCAGCATGGTAACGAAGCAGGTGATACCGGTAAGCGTCGTagc                                                                                                        >  2:12873/1‑146 (MQ=255)
                                        aCCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCt                                                                                                       <  7:115327/147‑1 (MQ=255)
                                              aGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCgg                                                                                                <  1:347808/148‑1 (MQ=255)
                                                      gaaggaagCAATAATCAGTACGTTCGCGCACAGTGCTGCCCAGGTAAATACTGTCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTgg                                                                                        <  1:12873/148‑1 (MQ=255)
                                                       aaggaagCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACCGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGga                                                                                       <  2:21858/148‑1 (MQ=255)
                                                        aggaagCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTg                                                                                         >  3:355242/1‑145 (MQ=255)
                                                        aggaagCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGgag                                                                                      <  7:313708/148‑1 (MQ=255)
                                                          gaagCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGAc                                                                                    >  3:132371/1‑148 (MQ=255)
                                                           aagCAATAATCAGTACGTTCGCGCACTGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATTCCCGGGGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACt                                                                                   <  2:270224/148‑1 (MQ=255)
                                                             gCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGGCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTcc                                                                                 <  8:94903/148‑1 (MQ=255)
                                                               aataatCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCa                                                                                >  4:145946/1‑147 (MQ=255)
                                                                ataatCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCata                                                                              >  3:24280/1‑148 (MQ=255)
                                                                 taatCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCatat                                                                             >  1:276101/1‑148 (MQ=255)
                                                                   atCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATcc                                                                           <  3:319416/148‑1 (MQ=255)
                                                                   atCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATcc                                                                           <  8:289450/148‑1 (MQ=255)
                                                                   atCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATc                                                                            >  3:109670/1‑147 (MQ=255)
                                                                   atCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATc                                                                            <  3:145946/147‑1 (MQ=255)
                                                                   atCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATc                                                                            <  4:24280/147‑1 (MQ=255)
                                                                    tCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCa                                                                          <  7:159675/148‑1 (MQ=255)
                                                                      aGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGACGTACCTATACCGGATAGCAGGAGACTCCATATCCAg                                                                         >  8:350236/1‑147 (MQ=255)
                                                                         aCGTTCGCGCACAGTTATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAAt                                                                     >  2:9269/1‑148 (MQ=255)
                                                                         aCGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAAt                                                                     >  4:167423/1‑148 (MQ=255)
                                                                                  cacaGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGa                                                            >  7:253996/1‑148 (MQ=255)
                                                                                    caGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCGATATCCAGTAATCGACACCGACt                                                          >  5:288948/1‑148 (MQ=255)
                                                                                    caGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCCTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACt                                                          >  2:58964/1‑148 (MQ=255)
                                                                                            cccAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACt                                                  >  5:348211/1‑148 (MQ=255)
                                                                                                        cTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTc                                       <  5:332847/147‑1 (MQ=255)
                                                                                                           gCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCg                                    <  2:276101/147‑1 (MQ=255)
                                                                                                               cTTGCACATGGTCATGCCCGCTGCGCGCCTCTTAAGTATGGTAACGAAGAATTTGATACCGGTAAGCGTCGTACATATACAGGAAAGCTTGATACTCCATATCCAGTAATCTACACCGACTCCCGTACTGTACTCTATTCCAGATatc                               >  6:259583/1‑146 (MQ=255)
                                                                                                                  gAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCgg                             <  2:362377/147‑1 (MQ=255)
                                                                                                                  gAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCgg                             <  2:284221/147‑1 (MQ=255)
                                                                                                                   aaCATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACCAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGt                            <  8:84973/147‑1 (MQ=255)
                                                                                                                    aCATGGTCATGCCCGGTGCGCGCATCCTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTg                           >  2:352853/1‑147 (MQ=255)
                                                                                                                     cATGGTCATGCCCGGTGCGCGCATCCTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGa                         <  1:352853/148‑1 (MQ=255)
                                                                                                                      aTGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGa                         >  1:65734/1‑147 (MQ=255)
                                                                                                                      aTGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGAt                        <  2:41792/148‑1 (MQ=255)
                                                                                                                      aTGGTCATGCCCGGTGCGCGCATCCTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGAt                        <  1:58964/148‑1 (MQ=255)
                                                                                                                       tGGTCATGCCCGGTGCGCGCATCGTCAGAATGGGAACGAAGAAGTTGATACCGGTAAGCGTCGTCCCTATACCGGATAGCTGGAGACTCCAGATCCAGTAAGCGCCGCCGACTCCCGGACGGGACTCGATTCCCGATAGCGGGGGa                         >  8:333995/1‑146 (MQ=255)
                                                                                                                          tCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGgc                    <  4:355242/148‑1 (MQ=255)
                                                                                                                           cATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGgcc                   >  6:373737/1‑148 (MQ=255)
                                                                                                                                       gcgcATCTTCAGTATGGTATCGAAGAAGTTGATACCTGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATCTCCAGTAATCGTCACCGACTCCCGGACTGTACTCTATTCCCGATATGGGTGGATAGGCCAGCCAGCCGGTc       <  6:288948/148‑1 (MQ=255)
                                                                                                                                         gtggcTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTg     >  4:339860/5‑148 (MQ=255)
                                                                                                                                             cTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTgcgc  >  4:325958/1‑147 (MQ=255)
                                                                                                                                             cTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTgcgc  >  5:69052/1‑147 (MQ=255)
                                                                                                                                              ttCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTccc                                     <  6:174360/111‑1 (MQ=255)
                                                                                                                                              ttCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTccc                                     >  5:174360/1‑111 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GGTGCCCAGATAGCGATCCAGGGTCAACAACGCGACGGTAACCGTCAGAATTGGGAAGGAAGCAATAATCAGTACGTTCGCGCACAGTGATGCCCAGGTAAATACTGGCATCTTGAACATGGTCATGCCCGGTGCGCGCATCTTCAGAATGGTAACGAAGAAGTTGATACCGGTAAGCGTCGTACCTATACCGGATAGCTGGAGACTCCATATCCAGTAATCGACACCGACTCCCGGACTGTACTCTATTCCCGATAGCGGTGGATAGGCCAGCCAGCCGGTCTGCGC  >  NC_000913/449859‑450146

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