Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,847,034 C→T 5.1% A127T (GCG→ACG)  hycD ← hydrogenase 3, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,847,0340CT 5.1% 253.7 / 2.9 79A127T (GCG→ACG) hycDhydrogenase 3, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (40/34);  new base T (2/2);  total (43/36)
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.34e-01

CCAGCGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGAACCGACGGTCACCACCGGCAGCGCAGTAGCGA  >  NC_000913/2846887‑2847180
                                                                                                                                                   |                                                                                                                                                  
ccAGCGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTcgc                                                                                                                                                    <  4:128426/148‑1 (MQ=255)
ccAGCGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTcgc                                                                                                                                                    <  1:12234/148‑1 (MQ=255)
 cAGCGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCgcgc                                                                                                                                                  <  2:65232/149‑1 (MQ=255)
    cGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGa                                                                                                                                               >  6:350864/1‑149 (MQ=255)
    cGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGa                                                                                                                                               <  3:366916/149‑1 (MQ=255)
    cGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGa                                                                                                                                               >  8:150149/1‑149 (MQ=255)
    cGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGa                                                                                                                                               >  2:149927/1‑149 (MQ=255)
       gATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAg                                                                                                                                            >  5:224898/1‑149 (MQ=255)
          gCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCgg                                                                                                                                         >  1:163454/1‑149 (MQ=255)
           ctctGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGGCCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGt                                                                                                                                        <  4:43205/149‑1 (MQ=255)
           ctctGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCgg                                                                                                                                         >  8:365129/1‑148 (MQ=255)
            tctGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGt                                                                                                                                        >  3:346532/1‑148 (MQ=255)
                gCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACg                                                                                                                                   <  2:357666/149‑1 (MQ=255)
                    aGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACggg                                                                                                                                 <  8:295422/147‑1 (MQ=255)
                    aGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACggg                                                                                                                                 <  6:224898/147‑1 (MQ=255)
                      cGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTAc                                                                                                                             >  5:4283/1‑149 (MQ=255)
                        gCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTTGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCg                                                                                                                           <  2:163454/149‑1 (MQ=255)
                            gTGATAAACGGTGTCGTTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTAt                                                                                                                       <  3:29008/149‑1 (MQ=255)
                            gTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTAt                                                                                                                       >  4:109171/1‑149 (MQ=255)
                             tGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATc                                                                                                                      <  3:244933/149‑1 (MQ=255)
                             tGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATc                                                                                                                      <  8:326433/149‑1 (MQ=255)
                               aTAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCCCAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCCAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGGATc                                                                                                                      >  2:245819/1‑147 (MQ=255)
                               aTAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATcca                                                                                                                    <  8:264614/149‑1 (MQ=255)
                               aTAAACGGTGTCGGGGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATcca                                                                                                                    <  2:242097/149‑1 (MQ=255)
                                 aaaCGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATccaga                                                                                                                  >  2:87496/1‑149 (MQ=255)
                                  aaCGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCACCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCgg                                                                                                                          >  2:42260/1‑140 (MQ=255)
                                  aaCGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCACCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCgg                                                                                                                          <  1:42260/140‑1 (MQ=255)
                                       tgtCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGaa                                                                                                            <  1:76941/149‑1 (MQ=255)
                                           ggTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATgg                                                                                                        >  1:75403/1‑149 (MQ=255)
                                           ggTGATGTTGCTCATGTTGGTGGAACCGGCACCCTGTGCGCCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATg                                                                                                         >  6:316926/1‑148 (MQ=255)
                                              gATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCa                                                                                                      >  3:318854/1‑148 (MQ=255)
                                               aTGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCaaa                                                                                                    <  2:75403/149‑1 (MQ=255)
                                               aTGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCa                                                                                                      >  8:143122/1‑147 (MQ=255)
                                                tggtGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAaag                                                                                                   <  1:387516/146‑1 (MQ=255)
                                                         tGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAAcgc                                                                                          >  2:143248/1‑149 (MQ=255)
                                                          gTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAAcgcg                                                                                         >  2:197563/1‑149 (MQ=255)
                                                              gtggAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGa                                                                                      <  3:323160/148‑1 (MQ=255)
                                                               tggAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATg                                                                                    <  7:150149/149‑1 (MQ=255)
                                                                 gAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGc                                                                                  <  1:143248/149‑1 (MQ=255)
                                                                  aaCCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGCCCAGAAATGGCAAAGAAGAGACGCGCGATGGc                                                                                  >  4:248508/1‑148 (MQ=255)
                                                                  aaCCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGc                                                                                  >  3:19545/1‑148 (MQ=255)
                                                                   aCCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGGCAAGCATCGCGTCACGGCTCGCGCCGATCGCGGGAAACGGGCGACCGGTAGCCAGCCCAGAAAGGGCAAAGAAGACACGCGCGATGGCa                                                                                 >  4:54186/1‑148 (MQ=255)
                                                                   aCCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCaa                                                                                <  5:350864/149‑1 (MQ=255)
                                                                    ccGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCaaa                                                                               <  5:136932/149‑1 (MQ=255)
                                                                    ccGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCa                                                                                 >  6:136932/1‑147 (MQ=255)
                                                                         aaCCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAgaga                                                                           >  3:12282/1‑148 (MQ=255)
                                                                         aaCCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGAt                                                                          <  7:143122/149‑1 (MQ=255)
                                                                                 cGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAg                                                                  <  4:19545/149‑1 (MQ=255)
                                                                                 cGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAg                                                                  <  1:24093/149‑1 (MQ=255)
                                                                                 cGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAg                                                                  >  5:356099/1‑149 (MQ=255)
                                                                                     aaCCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGa                                                              <  2:384320/149‑1 (MQ=255)
                                                                                     aaCCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGAGCCGCTAGCGGTAAACGTGCTACCGGTATCCAGACCACAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAg                                                                  >  5:163510/1‑145 (MQ=255)
                                                                                          acaGACCAAGCAGCAGCATCGGTTCGACCAGTACGCCAAGCATCGCTTCACGGCTCGTGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCaaa                                                         >  3:240067/1‑149 (MQ=255)
                                                                                          acaGACCAAGCAGCAGCATCGGTTCGACCAGTACGCCAAGCATCGCTTCACGGCTCGTGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCaa                                                          >  7:363143/1‑148 (MQ=255)
                                                                                           caGACCAAGCAGCAGCATCGGTTCGACCAGTACGCCAAGCATCGCTTCACGGCTCGTGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAAt                                                        <  4:240067/149‑1 (MQ=255)
                                                                                           caGACCAAGCAGCAGCATCGGTTCGACCAGTACGCCAAGCATCGCTTCACGGCTCGTGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAAt                                                        <  8:363143/149‑1 (MQ=255)
                                                                                            aGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATc                                                       <  7:365129/149‑1 (MQ=255)
                                                                                                cAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCAcc                                                    >  5:186340/1‑148 (MQ=255)
                                                                                                 aagaagCAGCATCGGTTCGACCATCAACCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACAACAAATGGCAAAGAAGCAACGCGCGATCGCAAAGACATACAGTAAGGTGATCaaa                                                         >  4:157612/1‑142 (MQ=255)
                                                                                                 aaGCAGCAGCATCGTTTCGACCAGCACGCCAATCATCGCTTCAATTATCGCGCCGATAGCGGTAAACTGGCTAACGTTATCCAGACCAGAAATGGCTAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCaaa                                                         <  3:157612/142‑1 (MQ=255)
                                                                                                   gcagcagcaTCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGt                                                >  1:24009/1‑149 (MQ=255)
                                                                                                      gcagcaTCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATc                                                            >  5:176884/1‑134 (MQ=255)
                                                                                                      gcagcaTCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATc                                                            <  6:176884/134‑1 (MQ=255)
                                                                                                      gcagcaTCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGc                                             >  8:102621/1‑149 (MQ=255)
                                                                                                         gcaTCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCgg                                           >  6:48115/1‑148 (MQ=255)
                                                                                                           aTCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAg                                        >  2:268330/1‑149 (MQ=255)
                                                                                                               gTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCg                                      >  1:87136/1‑147 (MQ=255)
                                                                                                               gTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGga                                    <  3:109171/149‑1 (MQ=255)
                                                                                                                ttCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGgag                                   <  4:346532/149‑1 (MQ=255)
                                                                                                                  cGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGaa                                 >  1:191235/1‑149 (MQ=255)
                                                                                                                   gACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGaa                                 <  7:70191/148‑1 (MQ=255)
                                                                                                                   gACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGaa                                 <  1:197563/148‑1 (MQ=255)
                                                                                                                   gACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGAAc                                <  6:186340/149‑1 (MQ=255)
                                                                                                                      cAGCACGCCAAGCATCGCTTCACGGCTCGCTCCGATAGCGGTAAACGGGCTACCGGTATCCTGACCAGAAATGTCAAAGAACAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGAACCGa                             <  5:250818/149‑1 (MQ=255)
                                                                                                                       aGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGCGCTACCGGTATCCAGACCAGAAATCGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGAACCTAc                            >  5:141781/1‑149 (MQ=255)
                                                                                                                            gCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGAACCGACGGTca                       <  1:87496/149‑1 (MQ=255)
                                                                                                                              cAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGCAATGGCAAAGAAGACACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGACCCGACGGTcccc                     >  6:7105/1‑149 (MQ=255)
                                                                                                                              cAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCACATCACCCACTTGCGGCAGCGGAGAACCGACGGTcacc                     >  6:375459/1‑149 (MQ=255)
                                                                                                                              cAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGAACCGACGGTcacc                     >  6:170056/1‑149 (MQ=255)
                                                                                                                              cAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGAACCGACGGTcac                      >  8:96842/1‑148 (MQ=255)
                                                                                                                                   aTCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGAACCGACGGTCACCACCgg                >  2:352509/1‑149 (MQ=255)
                                                                                                                                               cTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGAACCGACGGTCACCACCGGCAGCGCAGTAGc    >  2:81190/1‑149 (MQ=255)
                                                                                                                                                  gcgcCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGAACCGACGGTCACCACCGGCAGCGCAGTAGCGa  >  3:226955/1‑148 (MQ=255)
                                                                                                                                                  gcgcCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGAACCGACGGTCACCACCGGCAGCGCAGTAGCGa  <  5:141694/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCAGCGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGTTGGTGGAACCGGCAACCTGTGCGGCAACCCACAGACCAAGCAGCAGCATCGGTTCGACCAGCACGCCAAGCATCGCTTCACGGCTCGCGCCGATAGCGGTAAACGGGCTACCGGTATCCAGACCAGAAATGGCAAAGAAGAAACGCGCGATGGCAAAGAGATACAGTAAGGTGATCAAATCACCCAGTTGCGGCAGCGGAGAACCGACGGTCACCACCGGCAGCGCAGTAGCGA  >  NC_000913/2846887‑2847180

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