Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,604,990 G→C 11.5% T60T (ACG→ACC hyfD → hydrogenase 4 component D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,604,9900GC11.5% 260.9 / 22.7 87T60T (ACG→ACChyfDhydrogenase 4 component D
Reads supporting (aligned to +/- strand):  ref base G (40/37);  new base C (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 = 2.36e-01

CTGCCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATA  >  NC_000913/2604841‑2605138
                                                                                                                                                     |                                                                                                                                                    
ctgcCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCacgttgacg                                                                                                                                                      >  1:1552268/1‑150 (MQ=255)
    cTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTgg                                                                                                                                                  >  1:760128/1‑150 (MQ=255)
    cTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTgg                                                                                                                                                  >  2:1391874/1‑150 (MQ=255)
    cTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTgg                                                                                                                                                  >  2:1635185/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                                                                                                                                        <  2:88394/150‑3 (MQ=255)
              cgcACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACCAACGTCAACGTgac                                                                                                                                        <  1:514436/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:325095/150‑1 (MQ=255)
                                    ccACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTg                                                                                                                  >  2:1903939/1‑150 (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                                                                                                     >  1:1558253/1‑150 (MQ=255)
                                                 ccgccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATcgcg                                                                                                     <  2:321888/150‑1 (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: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                               <  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:174875/1‑150 (MQ=255)
                                                                                                                          ggCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                            <  2:1882232/150‑1 (MQ=255)
                                                                                                                          ggCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                            >  1:986763/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)
                                                                                                                                                     |                                                                                                                                                    
CTGCCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATA  >  NC_000913/2604841‑2605138

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: