Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 755,646 C→A 38.2% Y29* (TAC→TAA)  sdhD → succinate dehydrogenase, membrane subunit, binds cytochrome b556

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913755,6460CA38.2% 33.7 / 50.7 60Y29* (TAC→TAA) sdhDsuccinate dehydrogenase, membrane subunit, binds cytochrome b556
Reads supporting (aligned to +/- strand):  ref base C (20/17);  new base A (10/13);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

GCTCCGCCAAAATCTCCTTTGTTATTACTGTCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTCTATCTTGATCCATGCCTGGATCGGCATGTGG  >  NC_000913/755499‑755790
                                                                                                                                                   |                                                                                                                                                
gCTCCGCCAAAATCTCCTTTGTTATTACTGTCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAca                                                                                                                                                 >  4:291620/1‑149 (MQ=255)
    cGCCAAAATCTCCTTTGTTATTACTGTCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACTCTCTAcatcat                                                                                                                                             >  6:183178/1‑149 (MQ=255)
    cGCCAAAATCTCCTTTGTTATTACTGTCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAcatcat                                                                                                                                             <  3:313845/149‑1 (MQ=255)
      ccAAAATCTCCTTTGTTATTACTGTCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACTCTCTACATCAttt                                                                                                                                           <  5:183178/149‑1 (MQ=255)
      ccAAAATCTCCTTTGTTATTACTGTCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGAAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCAtt                                                                                                                                            >  6:110361/1‑148 (MQ=255)
         aaaTCTCCTTTGTTATTACTGTCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTtata                                                                                                                                        <  3:113535/149‑1 (MQ=255)
            tctcCTTTGTTATTACTGTCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATgg                                                                                                                                     <  1:363628/149‑1 (MQ=255)
                cTTTGTTATTACTGTCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGt                                                                                                                                    <  2:49702/146‑1 (MQ=255)
                cTTTGTTATTACTGTCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGt                                                                                                                                    <  7:204816/146‑1 (MQ=255)
                       attaCTGTCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTc                                                                                                                          >  6:213274/1‑149 (MQ=255)
                           cTGTCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTa                                                                                                                      >  3:207324/1‑149 (MQ=255)
                              tCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCa                                                                                                                   <  8:397889/149‑1 (MQ=255)
                                 tGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTg                                                                                                                <  3:247779/149‑1 (MQ=255)
                                    tttCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGc                                                                                                              >  6:288879/1‑148 (MQ=255)
                                       cACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAg                                                                                                           <  4:142251/148‑1 (MQ=255)
                                          ttCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGa                                                                                                       >  1:273006/1‑149 (MQ=255)
                                            ctcGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGAc                                                                                                      <  3:321349/148‑1 (MQ=255)
                                            ctcGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTAGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGa                                                                                                       <  5:110361/147‑1 (MQ=255)
                                                      tCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCt                                                                                           >  4:298684/1‑149 (MQ=255)
                                                      tCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCt                                                                                           >  1:122441/1‑149 (MQ=255)
                                                           gTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGAt                                                                                       <  4:37272/148‑1 (MQ=255)
                                                             aTGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCgg                                                                                    >  3:47334/1‑149 (MQ=255)
                                                                gTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGttt                                                                                 >  5:91130/1‑149 (MQ=255)
                                                                       aCGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGcc                                                                          >  5:246522/1‑149 (MQ=255)
                                                                        cGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGaa                                                                                               >  8:23797/1‑127 (MQ=255)
                                                                           cTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTg                                                                      >  8:67692/1‑149 (MQ=255)
                                                                            tCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCTGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGc                                                                           <  5:151412/143‑1 (MQ=255)
                                                                            tCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGc                                                                           >  6:151412/1‑143 (MQ=255)
                                                                                  ttAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGa                                                                                                >  4:97408/1‑116 (MQ=255)
                                                                                  ttAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGa                                                                                                <  3:97408/116‑1 (MQ=255)
                                                                                      gatgCATTGCCGTACATGATATCATCCTCGTGCGCGCTACCGCTATCGTCCTGACGCGCTCCATCAGTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCaaa                                                           <  6:33973/146‑1 (MQ=255)
                                                                                      gACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCaaa                                                           <  8:338959/149‑1 (MQ=255)
                                                                                           aaTGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGtt                                                      >  5:370406/1‑149 (MQ=255)
                                                                                             tGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCa                                                    >  2:346356/1‑149 (MQ=255)
                                                                                                  tACATCATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCatg                                               >  6:105961/1‑149 (MQ=255)
                                                                                                        aTTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGgc                                          <  5:213274/148‑1 (MQ=255)
                                                                                                             aTCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCaaa                                                           >  6:27261/1‑126 (MQ=255)
                                                                                                             aTCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCaaa                                                           <  5:27261/126‑1 (MQ=255)
                                                                                                             aTCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTg                                                                                                                <  1:34852/73‑1 (MQ=255)
                                                                                                             aTCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTg                                                                                                                >  2:34852/1‑73 (MQ=255)
                                                                                                              tCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGttt                                   <  2:122441/149‑1 (MQ=255)
                                                                                                                cTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTc                                 >  2:305540/1‑149 (MQ=255)
                                                                                                                cTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACAGATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTc                                 >  1:374870/1‑149 (MQ=255)
                                                                                                                cTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGtttt                                  <  6:91130/148‑1 (MQ=255)
                                                                                                                cTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTc                                 <  5:105961/149‑1 (MQ=255)
                                                                                                                 tCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGtttt                                  <  5:36786/147‑1 (MQ=255)
                                                                                                                 tCGTTCGCGCTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGtttt                                  >  6:36786/1‑147 (MQ=255)
                                                                                                                   gTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTCTAt                              >  1:249931/1‑149 (MQ=255)
                                                                                                                    ttCGCGCTACCGCTATCGTCCTGACGCTCTCCATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGCTCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGGTTTCTATc                             >  7:140036/1‑149 (MQ=255)
                                                                                                                         gcTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTCTATCTTGAt                        <  5:144175/149‑1 (MQ=255)
                                                                                                                          cTACCGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTCTATCTTGATc                       <  7:67692/149‑1 (MQ=255)
                                                                                                                           tACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTCTATCTTGAt                        <  2:249931/147‑1 (MQ=255)
                                                                                                                              cGCTATCGTCCTGACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTCTATCTTGATCCATg                   <  4:47334/149‑1 (MQ=255)
                                                                                                                                cTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCctg                                               <  2:338893/119‑1 (MQ=255)
                                                                                                                                cTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCctg                                               >  1:338893/1‑119 (MQ=255)
                                                                                                                                      tCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTCTATCTTGATCCATGCCTGGATc           <  5:288879/149‑1 (MQ=255)
                                                                                                                                          gACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTCTATCTTGATCCATGCCTGGATCGGCa       >  3:231891/1‑149 (MQ=255)
                                                                                                                                          gACGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTCTATCTTGATCCATGCCTGGATCGGCa       >  4:83337/1‑149 (MQ=255)
                                                                                                                                           aCGCTCTAAATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTACTGGCGCTGTTTTCTATCTTGATCCATGCCTGGATCGGCAt      >  6:228305/1‑149 (MQ=255)
                                                                                                                                            cGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTCTATCTTGATCCATGCCTGGATCGGCAtg     <  1:305540/149‑1 (MQ=255)
                                                                                                                                             gCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTCTATCTTGATCCATGCCTGGATCGGCAtgt    <  1:346356/149‑1 (MQ=255)
                                                                                                                                             gCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACACTGCTGGCGCTGTTTTCTATCTTGATCCATGCCTGGATCGGCAtgt    >  3:250202/1‑149 (MQ=255)
                                                                                                                                               tctACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTCTATCTTGATCCATGCCTGGATCGGCATGTgg  >  8:178754/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCTCCGCCAAAATCTCCTTTGTTATTACTGTCGTGCTTTCACTTCTCGCAGGAGTCCTCGTATGGTAAGCAACGCCTCCGCATTAGGACGCAATGGCGTACATGATTTCATCCTCGTTCGCGCTACCGCTATCGTCCTGACGCTCTACATCATTTATATGGTCGGTTTTTTCGCTACCAGTGGCGAGCTGACATATGAAGTCTGGATCGGTTTCTTCGCCTCTGCGTTCACCAAAGTGTTCACCCTGCTGGCGCTGTTTTCTATCTTGATCCATGCCTGGATCGGCATGTGG  >  NC_000913/755499‑755790

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