Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,835,127 A→C 8.2% intergenic (+82/+46) norW → / ← hypF NADH:flavorubredoxin oxidoreductase/carbamoyl phosphate phosphatase and [NiFe] hydrogenase maturation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,835,1270AC 8.2% 256.5 / 11.3 94intergenic (+82/+46)norW/hypFNADH:flavorubredoxin oxidoreductase/carbamoyl phosphate phosphatase and [NiFe] hydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base A (41/45);  new base C (2/6);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GCGCCTTTGTGGTCAGTGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCA  >  NC_000913/2834978‑2835276
                                                                                                                                                     |                                                                                                                                                     
gcgcCTTTGTGGTCAGTGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATaa                                                                                                                                                       >  2:1594302/1‑150 (MQ=255)
 cgcCTTTGTGGTCAGTGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAg                                                                                                                                                      >  1:933600/1‑150 (MQ=255)
  gcCTTTGTGGTCAGTGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAgg                                                                                                                                                     <  1:1101573/150‑1 (MQ=255)
       tgtgGTCAGAGAGGATCGGATGAAAGAGGAATTTGGATTGTTTAAAACATTGCCGATGTAGGTGGGCTACTGTTCCTAAAATGTCGGACGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGttt                                                                                                                                                <  2:1230526/150‑1 (MQ=38)
          ggTCAGTGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                             >  1:593777/1‑150 (MQ=255)
          ggTCAGTGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                             >  1:1360218/1‑150 (MQ=255)
           ggcattGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGGTTAcgc                                                                                                                                            <  2:1496778/145‑1 (MQ=255)
            tCAGTGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGCCCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                           >  2:1589043/1‑150 (MQ=255)
             cAGGGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                          <  2:294175/150‑1 (MQ=255)
                 gAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                      >  2:682393/1‑150 (MQ=255)
                  aGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                     <  2:1506967/150‑1 (MQ=38)
                    gATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                   <  2:1398014/150‑1 (MQ=37)
                    gATCGGATGAAAGAGGCATTTGGATGGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                   <  2:1547804/150‑1 (MQ=37)
                     ttCGGATGAAAGATGCATTTGGATTGTTGAAAACATCGCCGATGTTGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATGCGGCGTAAACGCCTTATCCggg                                                                                                                                  <  1:917468/149‑2 (MQ=25)
                      tCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCTTGTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                 <  1:530821/150‑2 (MQ=37)
                      tCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                 <  2:198341/150‑2 (MQ=37)
                        ggATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAt                                                                                                                               >  2:475063/1‑150 (MQ=255)
                        ggATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAt                                                                                                                               >  1:1658993/1‑150 (MQ=255)
                         gATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATg                                                                                                                              >  1:39281/1‑150 (MQ=255)
                         gATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGGGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATg                                                                                                                              <  1:1411160/150‑1 (MQ=255)
                            gAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGTGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgtg                                                                                                                           <  1:713831/150‑1 (MQ=255)
                             aaaGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgtgt                                                                                                                          >  2:1016936/1‑150 (MQ=255)
                             aaaGAGGCATTTGGATTGTTGAAAACATTGCCGATGGAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTAGTATCCGACCTACGGGGCCGCATGTGTAGGCCGGATACGGCGTTTACGCCGCATCCGGCAAgggggg                                                                                                                          >  2:277083/1‑147 (MQ=25)
                                gagGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCa                                                                                                                       >  2:1718214/1‑150 (MQ=255)
                                    cATTTGGATTGTTGAAAACATTGCCGATGTAGGTTGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCGACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATg                                                                                                                   <  2:755323/150‑1 (MQ=255)
                                    cATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATg                                                                                                                   <  1:291933/150‑1 (MQ=255)
                                    cATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATg                                                                                                                   <  2:830636/150‑1 (MQ=255)
                                        tGGATTGTTGAAAACATAGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAac                                                                                                               <  1:310210/150‑1 (MQ=255)
                                          gATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAacac                                                                                                             >  2:684841/1‑150 (MQ=255)
                                              gttgAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGttt                                                                                                         <  2:472786/150‑1 (MQ=255)
                                              gttgAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGttt                                                                                                         <  2:1002641/150‑1 (MQ=255)
                                              gttgAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGttt                                                                                                         <  1:748029/150‑1 (MQ=255)
                                               ttgAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGtttt                                                                                                        <  1:14532/150‑1 (MQ=255)
                                                tgAAAACATTGCCGATGTAGGTGGGCTACTGTGCCCAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATTCAACACGTTTTa                                                                                                       >  2:1152655/1‑150 (MQ=255)
                                                   aaaCATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATcc                                                                                                    <  2:420091/150‑1 (MQ=255)
                                                      cATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCTACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCCATGGTGTCCAAATTCAACACGTGTTATCCGtt                                                                                                 >  2:1571111/1‑150 (MQ=255)
                                                      cATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGATGGCGCGTCTTTTCCGACCTACGGTGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACGCGTTTTAGCCGtt                                                                                                 >  2:216369/1‑150 (MQ=37)
                                                       aTTGGCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTc                                                                                                <  2:1694689/150‑1 (MQ=255)
                                                       aTTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTc                                                                                                <  2:37803/150‑1 (MQ=255)
                                                          gCCGATGTAGGTGGGCTACTGTGCCTTATATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTgg                                                                                             <  2:395418/150‑1 (MQ=255)
                                                          ccccatgtatgtGCGCTACTGTGCCTAAAATGTCGTATTCGACGCTGGCGCGTCTTATCCGACCTAAGGGGACTCATGTGTAGGCCGGATAAGGCGTTTACGACGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTgg                                                                                             <  2:76871/146‑1 (MQ=25)
                                                           ccGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGa                                                                                            >  1:602771/1‑150 (MQ=255)
                                                            cGATGTAGGTGGGCTACTGTTCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGTCCTACGTTGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGAc                                                                                           <  1:834036/150‑1 (MQ=255)
                                                            cGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGAc                                                                                           <  2:264388/150‑1 (MQ=255)
                                                             gATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACt                                                                                          >  2:520448/1‑150 (MQ=255)
                                                             gATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACt                                                                                          >  2:190534/1‑150 (MQ=255)
                                                             gATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGGCCAAATGCAACACGTTTTATCCGTTCTGGACt                                                                                          <  1:126365/150‑1 (MQ=255)
                                                                 tAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCCACACGTTTTATCCGTTCTGGACTTCAc                                                                                      >  2:569363/1‑150 (MQ=255)
                                                                  aGGGGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCAcc                                                                                     <  2:393707/150‑1 (MQ=255)
                                                                     tGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTATATCCGTTCTGGACATCACCCGc                                                                                  >  1:17472/1‑150 (MQ=255)
                                                                       ggCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTa                                                                                <  2:25407/150‑1 (MQ=255)
                                                                        gCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACCCGGTTTATCCGTTCTGGACTTCGCCCGCGaa                                                                               >  1:1708270/1‑150 (MQ=255)
                                                                          tACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAAcc                                                                             <  2:588964/150‑1 (MQ=255)
                                                                                 ccTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACgcgc                                                                      >  1:398922/1‑150 (MQ=255)
                                                                                  cTAAAATGTCGGATGCGACGCTGGCGCGTCTTATACGACCTACGGGGACGAATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAACTGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGcgcc                                                                     >  1:1498763/1‑150 (MQ=255)
                                                                                   tAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGcgccg                                                                    <  2:1443571/150‑1 (MQ=255)
                                                                                     aaaTGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCgca                                                                  >  1:1473870/1‑150 (MQ=255)
                                                                                          tCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAAt                                                             <  2:1150475/150‑1 (MQ=255)
                                                                                           cGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAAGACGTTTTATCCGTTCTGGACTTCGCCCGCGAACCAACGCGCCGCAGCAATa                                                            >  1:1133262/1‑150 (MQ=255)
                                                                                             gATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATTCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCTGCAATAAc                                                          <  1:1152655/150‑1 (MQ=255)
                                                                                              aTGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAAcc                                                         <  2:990649/150‑1 (MQ=255)
                                                                                               tGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAAccc                                                        >  2:479745/1‑150 (MQ=255)
                                                                                               tGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAAccc                                                        <  1:144902/150‑1 (MQ=255)
                                                                                                 cGACGCTGGCGCGTCTTATACGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCATCACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAAccccc                                                      >  2:636380/1‑150 (MQ=255)
                                                                                                          cgcgTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCaga                                             >  2:1665732/1‑150 (MQ=255)
                                                                                                          cgcgTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGGCCAAATGCCCCACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCagc                                             >  2:1506757/1‑149 (MQ=255)
                                                                                                           gcgTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCagag                                            <  2:475738/150‑1 (MQ=255)
                                                                                                             gTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGAc                                          <  2:1360218/150‑1 (MQ=255)
                                                                                                             gTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGAc                                          >  2:1379140/1‑150 (MQ=255)
                                                                                                             gTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGAc                                          <  2:593777/150‑1 (MQ=255)
                                                                                                              tCTTATCCGACCTACGGGGACGCATGTGTAGGTCGGTTTAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACa                                         <  2:132363/150‑1 (MQ=255)
                                                                                                                  aTCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAccc                                     >  1:115868/1‑150 (MQ=255)
                                                                                                                  aTCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAAcc                                     <  1:682393/150‑1 (MQ=255)
                                                                                                                     cGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGGTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGcc                                  <  2:730714/150‑1 (MQ=255)
                                                                                                                          tACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCAc                             >  1:395265/1‑150 (MQ=255)
                                                                                                                           aCGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCAcc                            <  1:1383610/150‑1 (MQ=255)
                                                                                                                           aCGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCAcc                            >  2:1461873/1‑150 (MQ=255)
                                                                                                                            cGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCGTCCCCCAGAGCCAACCCGCCAACcccc                           >  2:819361/1‑150 (MQ=38)
                                                                                                                               ggACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACccgcc                        <  1:981300/150‑1 (MQ=255)
                                                                                                                                gACGCATGTGTAGGCCGGATAATGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACccgccg                       <  2:1215126/150‑1 (MQ=255)
                                                                                                                                gACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACccgccg                       <  2:1658993/150‑1 (MQ=255)
                                                                                                                                    cATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCTTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTaa                   <  2:296684/150‑1 (MQ=255)
                                                                                                                                    cATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTaa                   <  2:39281/150‑1 (MQ=255)
                                                                                                                                           aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCtgt            >  1:1629116/1‑150 (MQ=255)
                                                                                                                                           aGGCCGGATAAGGCGTTTACGCCGCAGCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCtgt            <  1:347750/150‑1 (MQ=255)
                                                                                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGCCAATCCGCCATCCCCCGCCGGGAAACTCtggg           >  1:1709451/1‑148 (MQ=255)
                                                                                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCtgtg           >  1:1716120/1‑150 (MQ=255)
                                                                                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCtgtg           >  1:1334208/1‑150 (MQ=255)
                                                                                                                                             gCCGGATGCGGCGTAAACGCCTTATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATACCCCCCTGCCCCAGAGCCAAACCGCCATCACCCGCCGGTAAACTCTGTgg          >  2:879253/1‑150 (MQ=255)
                                                                                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCGGAGACAAACCGCCCTCACCCGCCGGTAAACTCTGTgg          >  2:60645/1‑150 (MQ=255)
                                                                                                                                              ccGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGCCTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGa         >  1:1659000/1‑150 (MQ=255)
                                                                                                                                               cGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTGCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCAGCCGGTAAACTCTGTGGaa        >  2:321355/1‑150 (MQ=255)
                                                                                                                                                ggATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGaaa       >  2:1371963/1‑150 (MQ=255)
                                                                                                                                                 gATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAag      <  2:933600/150‑1 (MQ=255)
                                                                                                                                                  aTAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAaga     <  1:1377917/150‑1 (MQ=255)
                                                                                                                                                     aGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCCCCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCACAGAGACAAACCGACATCACCCGCCGGTAAACTCTGTGGATAGAGCa  >  1:1685405/1‑150 (MQ=255)
                                                                                                                                                     aGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCCTCACCCGCCGGGAAACTCTGGGGAAAGAGCa  >  2:1076725/1‑150 (MQ=255)
                                                                                                                                                     aGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCa  >  2:1285487/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGCCTTTGTGGTCAGTGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCA  >  NC_000913/2834978‑2835276

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