Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 757,972:1 +A 7.2% coding (284/717 nt) sdhB → succinate dehydrogenase, FeS subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913757,9721.A 7.2% 313.7 / 22.3 83coding (284/717 nt)sdhBsuccinate dehydrogenase, FeS subunit
Reads supporting (aligned to +/- strand):  ref base . (39/38);  new base A (3/3);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

TCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGA  >  NC_000913/757826‑758119
                                                                                                                                                   |                                                                                                                                                    
tCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑a                                                                                                                                                    <  8:104214/148‑1 (MQ=255)
      ccTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTa                                                                                                                                              <  1:222437/148‑1 (MQ=255)
          tCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTAcc                                                                                                                                            <  5:206672/146‑1 (MQ=255)
          tCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCgg                                                                                                                                          <  1:228412/148‑1 (MQ=255)
            gTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGt                                                                                                                                         <  3:321004/147‑1 (MQ=255)
               gcgccgcTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAATATTCGTCTGGCCTGTATTACCCCGATTTCGGCACGCAACCTGCCGGGCAAGAAGATTGGGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGAt                                                                                                                                      <  5:82562/146‑1 (MQ=255)
               ccgccgCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGAt                                                                                                                                      <  3:414726/147‑1 (MQ=255)
                cgccgcTCCTGCCGTGAAGGTGGGGGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATcc                                                                                                                                    <  6:21838/148‑1 (MQ=255)
                   cgcTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑a                                                                                                                                                    >  7:84730/1‑129 (MQ=255)
                   cgcTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑a                                                                                                                                                    <  8:84730/129‑1 (MQ=255)
                       ccTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGAttt                                                                                                                             <  1:75557/148‑1 (MQ=255)
                          gCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGAAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGAttt                                                                                                                             >  4:137401/1‑145 (MQ=255)
                           ccGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTtg                                                                                                                            <  3:137401/145‑1 (MQ=255)
                                aaGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAg                                                                                                                     <  2:54819/147‑1 (MQ=255)
                                 aGGTGTGGGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGAc                                                                                                                   <  6:53767/148‑1 (MQ=255)
                                   gtgtgtGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACAt                                                                                                                 <  1:281355/148‑1 (MQ=255)
                                    tgtgtgCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATg                                                                                                                <  5:425196/148‑1 (MQ=255)
                                    tgtgtgCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATg                                                                                                                <  8:254204/148‑1 (MQ=255)
                                         gCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGAc                                                                                                            >  1:351753/1‑147 (MQ=255)
                                          cGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACAtaacagc                                                                                                           >  1:166059/1‑141 (MQ=255)
                                          cGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACa                                                                                                           <  7:446656/147‑1 (MQ=255)
                                           ggTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAAt                                                                                                         <  1:160702/148‑1 (MQ=255)
                                               ccGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCt                                                                                                      <  7:22159/147‑1 (MQ=255)
                                               ccGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATACGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCt                                                                                                      <  6:140321/147‑1 (MQ=255)
                                               ccGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTGCGGCACGCTACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGCTCCGCGATTTGGTGGTAGACATGGGACAATTCTa                                                                                                     <  5:96307/148‑1 (MQ=255)
                                                cGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCt                                                                                                      <  8:87765/148‑1 (MQ=255)
                                                       cTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAat                                                                                             <  4:234334/148‑1 (MQ=255)
                                                        tGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCt                                                                                                      <  3:21701/140‑1 (MQ=255)
                                                        tGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCt                                                                                                      >  4:21701/1‑140 (MQ=255)
                                                         gAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACATTTCTATGCGAAatat                                                                                           >  6:52190/1‑148 (MQ=255)
                                                          aaCATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATg                                                                                          <  2:113713/148‑1 (MQ=255)
                                                           aCATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAatat                                                                                           >  6:142723/1‑146 (MQ=255)
                                                               gAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGaa                                                                                      >  6:304388/1‑147 (MQ=255)
                                                                  cGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAAt                                                                                   >  3:220879/1‑147 (MQ=255)
                                                                   ggCAAGAATGGGCTGGCCTGTATTACCCCGATTTCGGCACTCATCCAGGCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTTGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTa                                                                                 <  6:224037/148‑1 (MQ=255)
                                                                     cAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAg                                                                               <  2:351753/148‑1 (MQ=255)
                                                                      aataaTGTTCTGGCCTGTATTACCCCGATTTCTGCACTCAACCAGCCGTTCAAGAAGATTGTGATTCGCCCGCTGCA‑‑AGGTTTACAAGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGAGCAATATGAGAAAATTAAg                                                                               <  5:185674/144‑1 (MQ=255)
                                                                      aaGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTa                                                                                 >  8:332242/1‑145 (MQ=255)
                                                                      aaGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTa                                                                                 <  7:332242/145‑1 (MQ=255)
                                                                        gAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCt                                                                            >  5:227593/1‑148 (MQ=255)
                                                                          aTGTTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCTGTGATCCTCGATTTGGTGGTAGACATGGGACAATTCTATTATCAATATGAGAAAATTAAGCAtt                                                                           >  4:102932/1‑147 (MQ=255)
                                                                               cTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGt                                                                     >  8:267541/1‑148 (MQ=255)
                                                                                  gCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGtt                                                                    <  5:142723/146‑1 (MQ=255)
                                                                                  gCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTg                                                                   >  8:276531/1‑147 (MQ=255)
                                                                                      gTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGaataa                                                              >  7:53859/1‑148 (MQ=255)
                                                                                       tATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGaata                                                               <  5:304388/146‑1 (MQ=255)
                                                                                          tACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTAcc                                                                        <  8:249343/134‑1 (MQ=255)
                                                                                          tACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTAcc                                                                        >  7:249343/1‑134 (MQ=255)
                                                                                          tACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGa                                                          >  6:187951/1‑148 (MQ=255)
                                                                                          tACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGAGTTGGTGGTAGACATGGGACACTTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGa                                                          >  3:224554/1‑148 (MQ=255)
                                                                                            ccccGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGAc                                                         >  6:457622/1‑147 (MQ=255)
                                                                                             cccGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATgg                                                           >  1:133162/1‑144 (MQ=255)
                                                                                             cccGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATgg                                                           <  2:133162/144‑1 (MQ=255)
                                                                                             cccGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACa                                                        >  5:28774/1‑147 (MQ=255)
                                                                                              ccGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACaa                                                       >  7:447885/1‑147 (MQ=255)
                                                                                              ccGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACaa                                                       >  4:40296/1‑147 (MQ=255)
                                                                                               cGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACa                                                        >  3:437890/1‑147 (MQ=255)
                                                                                               cGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACaaa                                                      >  8:465165/1‑147 (MQ=255)
                                                                                                   ttCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATc                                                   >  5:321095/1‑148 (MQ=255)
                                                                                                      ggCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCg                                                 <  6:321095/147‑1 (MQ=255)
                                                                                                       gCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAg                                             >  5:333206/1‑148 (MQ=255)
                                                                                                        cACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTTATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGa                                                                  <  3:131656/126‑1 (MQ=255)
                                                                                                         aCTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATTGGACAATTCTATGCGCAATATGAGAAAATTAATCCTTACCTGTAGaa                                                                 >  4:131656/1‑126 (MQ=255)
                                                                                                           tCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTcg                                         >  5:358023/1‑148 (MQ=255)
                                                                                                           tCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTcg                                         >  8:383451/1‑148 (MQ=255)
                                                                                                            cAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTcgc                                        >  8:305556/1‑148 (MQ=255)
                                                                                                            cAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTcgc                                        >  1:216605/1‑148 (MQ=255)
                                                                                                            cAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAAGATGAGAAAATTAAGCCGTACCTGGTGAAGAATGGACAAAATCCGCCAGCTcgc                                        >  1:223060/1‑148 (MQ=255)
                                                                                                              aCCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTcgc                                        >  6:94803/1‑146 (MQ=255)
                                                                                                               ccAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTtggtg                                                                                                                         >  1:133490/1‑64 (MQ=255)
                                                                                                               ccAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTtggtg                                                                                                                         <  2:133490/64‑1 (MQ=255)
                                                                                                                      ggCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTAc                              <  5:187951/148‑1 (MQ=255)
                                                                                                                       gCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTa                               >  3:280330/1‑146 (MQ=255)
                                                                                                                        cAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTagag                            <  6:227593/148‑1 (MQ=255)
                                                                                                                        cAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAg                            <  3:40296/148‑1 (MQ=255)
                                                                                                                        cAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAg                            <  7:465165/148‑1 (MQ=255)
                                                                                                                          agaagaTTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGAATTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTGACAGAt                          >  4:273678/1‑148 (MQ=255)
                                                                                                                            aagaTTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTTATCCGAGATTTGTTTGAAGACATGGGACAATTATATGCGCAATATGACATAATTAACCCTTACCTGTTGAATAATTGACAAAATCCGACAGCTCGCGAGAATTTAAAGATg                         >  3:353887/1‑147 (MQ=255)
                                                                                                                              gaTTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGc                        >  6:295879/1‑146 (MQ=255)
                                                                                                                                ttGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCa                      >  7:271122/1‑146 (MQ=255)
                                                                                                                                        cGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTAc                                                                         >  7:367732/1‑87 (MQ=255)
                                                                                                                                        cGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAgcgc             <  4:280330/147‑1 (MQ=255)
                                                                                                                                          cccGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCACAGCAGCGCGa           >  1:122263/1‑147 (MQ=255)
                                                                                                                                           ccGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGaa          >  4:342194/1‑147 (MQ=255)
                                                                                                                                            cGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGaaaa        <  7:383451/148‑1 (MQ=255)
                                                                                                                                             gCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGGAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGaaaaa       <  6:28774/148‑1 (MQ=255)
                                                                                                                                                 cc‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGaaaactcg    <  6:358023/147‑5 (MQ=255)
                                                                                                                                                 cc‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACt     <  6:333206/146‑1 (MQ=255)
                                                                                                                                                 cc‑‑AGGTTAACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACt     <  5:457622/146‑1 (MQ=255)
                                                                                                                                                  c‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGa  >  1:48261/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGGTCTGAACATGAACGGCAAGAATGGTCTGGCCTGTATTACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCC‑‑AGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGA  >  NC_000913/757826‑758119

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