Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,604,993 G→C 11.9% L61F (TTG→TTC) ‡ hyfD → hydrogenase 4 component D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,604,9930GC11.9% 243.2 / 22.9 84L61F (TTG→TTC) ‡hyfDhydrogenase 4 component D
Reads supporting (aligned to +/- strand):  ref base G (37/36);  new base C (6/4);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

CTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGG  >  NC_000913/2604845‑2605142
                                                                                                                                                    |                                                                                                                                                     
cTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTgg                                                                                                                                                      >  1:760128/1‑150 (MQ=255)
cTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTgg                                                                                                                                                      >  2:1635185/1‑150 (MQ=255)
cTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTgg                                                                                                                                                      >  2:1391874/1‑150 (MQ=255)
    tATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCaa                                                                                                                                                  >  1:1139736/1‑150 (MQ=255)
       cGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACcaacgtcaacgt                                                                                                                                               <  2:330395/150‑1 (MQ=255)
        ggCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACCAACGTCAACGTg                                                                                                                                              <  1:956483/150‑1 (MQ=255)
          cgcACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACCAACGTCAACGTgac                                                                                                                                            <  1:514436/150‑3 (MQ=255)
          cgcACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACCAACGTCAACGTgac                                                                                                                                            <  2:88394/150‑3 (MQ=255)
          cgcACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTATGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACCAACGTCAACGTgac                                                                                                                                            <  1:1077202/150‑3 (MQ=255)
             aCTGGTCGTTTCGTTTTCGCCACAACGGCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGa                                                                                                                                         <  2:180842/150‑1 (MQ=255)
                ggTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTTATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGAtgt                                                                                                                                      >  1:740503/1‑150 (MQ=255)
                 gTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGAtgtg                                                                                                                                     >  1:451618/1‑150 (MQ=255)
                   cgtttcgtttTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGc                                                                                                                                   >  1:321888/1‑150 (MQ=255)
                   cgtttcgtttTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGc                                                                                                                                   >  2:854605/1‑150 (MQ=255)
                     tttcgtttTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCgt                                                                                                                                 >  1:140492/1‑150 (MQ=255)
                       tcgtttTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCgttg                                                                                                                               >  1:784670/1‑150 (MQ=255)
                        cgtttTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCgttgt                                                                                                                              >  1:918486/1‑150 (MQ=255)
                            ttCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTgg                                                                                                                          <  1:977325/150‑1 (MQ=255)
                               gCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCt                                                                                                                       >  1:1639116/1‑150 (MQ=255)
                                ccACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTg                                                                                                                      >  2:1903939/1‑150 (MQ=255)
                                ccACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTg                                                                                                                      <  2:325095/150‑1 (MQ=255)
                                      cGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCAtt                                                                                                                <  2:1192489/150‑1 (MQ=255)
                                        tCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGa                                                                                                              <  2:1639116/150‑1 (MQ=255)
                                          gggCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACGACGCTGTGCATGTTGTCACTGATCTGCGCGTTTTATGAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGTTCATTGATc                                                                                                            <  2:676020/150‑1 (MQ=255)
                                          gggCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATc                                                                                                            >  1:1344471/1‑150 (MQ=255)
                                             ccgccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATcgcg                                                                                                         >  1:444511/1‑150 (MQ=255)
                                             ccgccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATcgcg                                                                                                         <  2:321888/150‑1 (MQ=255)
                                             ccgccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATcgcg                                                                                                         >  1:1558253/1‑150 (MQ=255)
                                               gccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTTTCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTCTGGCCTGGTCATTGATCGCGTg                                                                                                       <  1:1041257/150‑1 (MQ=255)
                                               gccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTTAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTg                                                                                                       <  1:1768498/150‑1 (MQ=255)
                                               gccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTg                                                                                                       >  2:226246/1‑150 (MQ=255)
                                               gccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTg                                                                                                       <  2:784670/150‑1 (MQ=255)
                                                     tGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACCGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCt                                                                                                                       >  1:1257317/1‑128 (MQ=255)
                                                     tGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACCGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCt                                                                                                                       <  2:1257317/128‑1 (MQ=255)
                                                      ggggggTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGc                                                                                                <  1:1200494/150‑1 (MQ=255)
                                                      ggggggTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGc                                                                                                >  2:1436444/1‑150 (MQ=255)
                                                      ggggggTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGc                                                                                                >  2:1429408/1‑150 (MQ=255)
                                                       gggggTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCt                                                                                               <  2:451618/150‑1 (MQ=255)
                                                        ggggTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTg                                                                                              <  1:854605/150‑1 (MQ=255)
                                                         gggTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGa                                                                                             <  2:444511/150‑1 (MQ=255)
                                                            ttttGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTc                                                                                          <  1:1675689/150‑1 (MQ=255)
                                                              ttGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTg                                                                                        <  1:7941/150‑1 (MQ=255)
                                                                gTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGtt                                                                                      >  1:211503/1‑150 (MQ=255)
                                                                gTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGtt                                                                                      >  1:437069/1‑150 (MQ=255)
                                                                   cgccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTtgt                                                                                   >  2:204875/1‑150 (MQ=255)
                                                                    gccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTtgtg                                                                                  >  2:1986724/1‑150 (MQ=255)
                                                                    gccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTtgtg                                                                                  >  2:1840492/1‑150 (MQ=255)
                                                                    gccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTtgtg                                                                                  <  2:309036/150‑1 (MQ=255)
                                                                            gACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGtt                                                                          >  1:895890/1‑150 (MQ=255)
                                                                                aCGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTtctc                                                                      <  1:794513/150‑1 (MQ=255)
                                                                                   cTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGt                                                                   <  1:215288/150‑1 (MQ=255)
                                                                                     gtgCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGTGGGATGTGGCGTTGTTTGGCATGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTTGTGTTTCTCTGttt                                                                 <  2:760128/150‑1 (MQ=255)
                                                                                       gCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGc                                                               <  2:140492/150‑1 (MQ=255)
                                                                                          tgttgtCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTgg                                                            >  1:344310/1‑150 (MQ=255)
                                                                                           gttgtCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGt                                                           >  1:1919708/1‑150 (MQ=255)
                                                                                            ttgtCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTc                                                          <  2:417304/150‑1 (MQ=255)
                                                                                             tgtCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCa                                                         <  2:170949/150‑1 (MQ=255)
                                                                                               tCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACg                                                       <  2:1344471/150‑1 (MQ=255)
                                                                                                   tGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCt                                                   <  1:1391874/150‑1 (MQ=255)
                                                                                                    gATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTa                                                  >  1:977468/1‑150 (MQ=255)
                                                                                                       ctcCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTc                                               >  2:353178/1‑150 (MQ=255)
                                                                                                        tcCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTcc                                              <  2:1016445/150‑1 (MQ=255)
                                                                                                           gcgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACg                                           >  1:1882232/1‑150 (MQ=255)
                                                                                                           gcgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACg                                           >  2:842891/1‑150 (MQ=255)
                                                                                                            cgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTTGTCAACGTGGGGGATGTGTCGTTGTTTGGCCTTGTCATTGATCGCGTGTGTACGCTGATTCTGTTTGTGGTGGTGTTTGTCGGTTTGCTGGTCACGATCTACTCCACgg                                          <  2:1139736/150‑1 (MQ=255)
                                                                                                                  aTCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATc                                    <  2:418050/150‑1 (MQ=255)
                                                                                                                   tCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCt                                   <  2:344310/150‑1 (MQ=255)
                                                                                                                   tCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCt                                   <  1:841112/150‑1 (MQ=255)
                                                                                                                    cAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTg                                  <  2:895890/150‑1 (MQ=255)
                                                                                                                     ggcctgATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGa                                 >  1:1640920/6‑150 (MQ=255)
                                                                                                                      ggCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                                >  1:986763/1‑150 (MQ=255)
                                                                                                                      ggCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                                <  2:1882232/150‑1 (MQ=255)
                                                                                                                      ggCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                                >  1:174875/1‑150 (MQ=255)
                                                                                                                      ggCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCAGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                                >  2:909309/1‑150 (MQ=255)
                                                                                                                        cctaTAAAGTTGCCGTGACGTTGTGGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGGTGATTCTGTTTGTGGTTGTGTTTCTCGGTTTGCTGGTGACGATCTACTCCACTGGTTATCTGActt                              <  2:9409/147‑3 (MQ=255)
                                                                                                                        ccGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACgg                              <  2:740503/150‑1 (MQ=255)
                                                                                                                        ccGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACgg                              >  1:1495060/1‑150 (MQ=255)
                                                                                                                            tAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATaa                          >  1:1226501/1‑150 (MQ=255)
                                                                                                                                    ccGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTGTCTGACGGATAAAAATcgcg                  <  2:918486/150‑1 (MQ=255)
                                                                                                                                    ccGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATcgcg                  <  1:1903939/150‑1 (MQ=255)
                                                                                                                                     ccccACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGa                 >  2:250630/4‑150 (MQ=255)
                                                                                                                                     ccccACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGa                 >  2:1569044/4‑150 (MQ=255)
                                                                                                                                     cGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGa                 >  1:1005006/1‑150 (MQ=255)
                                                                                                                                     cGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTCATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGa                 >  1:428962/1‑150 (MQ=255)
                                                                                                                                      cccACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGaa                >  1:1326641/3‑150 (MQ=255)
                                                                                                                                      cccACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGaa                >  1:1133474/3‑150 (MQ=255)
                                                                                                                                       ccACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAac               >  2:608399/2‑150 (MQ=255)
                                                                                                                                       ccACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAac               >  1:118739/2‑150 (MQ=255)
                                                                                                                                              gacgttgGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCa        <  1:226246/150‑1 (MQ=255)
                                                                                                                                                cgttgGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATa      >  1:913978/1‑150 (MQ=255)
                                                                                                                                                  ttgGTCAACGTGGGGGAGGTGGCGGTGTTTGGCCTGGTCATTGATCGCGTGAGGACGCGGATTCTGTTTGTGGGGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAAc    <  2:1254242/150‑1 (MQ=255)
                                                                                                                                                    ggTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACgg  <  2:437069/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGG  >  NC_000913/2604845‑2605142

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