Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,604,992 T→A 11.4% L61* (TTG→TAG) ‡ hyfD → hydrogenase 4 component D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,604,9920TA11.4% 242.6 / 17.3 87L61* (TTG→TAG) ‡hyfDhydrogenase 4 component D
Reads supporting (aligned to +/- strand):  ref base T (40/37);  new base A (6/4);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAAC  >  NC_000913/2604845‑2605140
                                                                                                                                                   |                                                                                                                                                    
cTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTgg                                                                                                                                                    >  2:1635185/1‑150 (MQ=255)
cTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTgg                                                                                                                                                    >  2:1391874/1‑150 (MQ=255)
cTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTgg                                                                                                                                                    >  1:760128/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                                                                                                       <  2:321888/150‑1 (MQ=255)
                                             ccgccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATcgcg                                                                                                       >  1:444511/1‑150 (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:784670/150‑1 (MQ=255)
                                               gccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTg                                                                                                     >  2:226246/1‑150 (MQ=255)
                                                     tGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACCGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCt                                                                                                                     <  2:1257317/128‑1 (MQ=255)
                                                     tGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACCGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCt                                                                                                                     >  1:1257317/1‑128 (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:437069/1‑150 (MQ=255)
                                                                gTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGtt                                                                                    >  1:211503/1‑150 (MQ=255)
                                                                   cgccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTtgt                                                                                 >  2:204875/1‑150 (MQ=255)
                                                                    gccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTtgtg                                                                                <  2:309036/150‑1 (MQ=255)
                                                                    gccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTtgtg                                                                                >  2:1986724/1‑150 (MQ=255)
                                                                    gccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTtgtg                                                                                >  2:1840492/1‑150 (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                                         >  2:842891/1‑150 (MQ=255)
                                                                                                           gcgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACg                                         >  1:1882232/1‑150 (MQ=255)
                                                                                                            cgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTTGTCAACGTGGGGGATGTGTCGTTGTTTGGCCTTGTCATTGATCGCGTGTGTACGCTGATTCTGTTTGTGGTGGTGTTTGTCGGTTTGCTGGTCACGATCTACTCCACgg                                        <  2:1139736/150‑1 (MQ=255)
                                                                                                                  aTCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATc                                  <  2:418050/150‑1 (MQ=255)
                                                                                                                   tCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCt                                 <  1:841112/150‑1 (MQ=255)
                                                                                                                   tCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCt                                 <  2:344310/150‑1 (MQ=255)
                                                                                                                    cAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTg                                <  2:895890/150‑1 (MQ=255)
                                                                                                                     ggcctgATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGa                               >  1:1640920/6‑150 (MQ=255)
                                                                                                                      ggCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                              <  2:1882232/150‑1 (MQ=255)
                                                                                                                      ggCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                              >  1:986763/1‑150 (MQ=255)
                                                                                                                      ggCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                              >  1:174875/1‑150 (MQ=255)
                                                                                                                      ggCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCAGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                              >  2:909309/1‑150 (MQ=255)
                                                                                                                        cctaTAAAGTTGCCGTGACGTTGTGGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGGTGATTCTGTTTGTGGTTGTGTTTCTCGGTTTGCTGGTGACGATCTACTCCACTGGTTATCTGActt                            <  2:9409/147‑3 (MQ=255)
                                                                                                                        ccGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACgg                            >  1:1495060/1‑150 (MQ=255)
                                                                                                                        ccGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACgg                            <  2:740503/150‑1 (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)
                                                                                                                                                   |                                                                                                                                                    
CTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAAC  >  NC_000913/2604845‑2605140

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: