Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,629,219 G→A 5.8% intergenic (‑281/‑71) yfgF ← / → yfgG cyclic‑di‑GMP phosphodiesterase, anaerobic/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,629,2190GA 5.8% 216.6 / 2.0 69intergenic (‑281/‑71)yfgF/yfgGcyclic‑di‑GMP phosphodiesterase, anaerobic/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (35/30);  new base A (2/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TATATCAAACAAAATAACCCTGATTAATGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTT  >  NC_000913/2629072‑2629366
                                                                                                                                                    |                                                                                                                                                   
tatatCAAACAAAATAACCCTGATTAATGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGa                                                                                                                                                    >  4:177594/1‑149 (MQ=255)
    tCAAACAAAATAACCCTGATTAATGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCAt                                                                                                                                                >  4:102685/1‑149 (MQ=255)
      aaacaaaATAACCCTGATTAATGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGATGAATTGACCAtt                                                                                                                                               <  8:227052/148‑1 (MQ=255)
       aacaaaATAACCCTGATTAATGAATTATTACGTTTATCATGTTAATTAATCATTATTACAACATCATTGTAAATAATTAAATTAACTTCCATAACATTAAATTATGTATCCACTGACGCTTTTTAACATAACTAAGAATTGACCAtttt                                                                                                                                             >  6:125085/1‑149 (MQ=255)
        acaaaATAACCCTGATTAATGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTg                                                                                                                                            >  5:229845/1‑149 (MQ=255)
              tAACCCTGATTAATGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGa                                                                                                                                                    <  4:220767/135‑1 (MQ=255)
              tAACCCTGATTAATGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGa                                                                                                                                                    >  3:220767/1‑135 (MQ=255)
                    tGATTAATGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCAACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGATATGc                                                                                                                                 >  4:41128/1‑148 (MQ=255)
                      aTTAATGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGcc                                                                                                                                <  7:143744/147‑1 (MQ=255)
                       ttAATGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCtt                                                                                                                              >  7:120172/1‑148 (MQ=255)
                       ttAATGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTa                                                                                                                             <  5:66452/149‑1 (MQ=255)
                          aTGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAAt                                                                                                                           <  8:47299/148‑1 (MQ=255)
                           tGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACAATTTTGTCCTGTTGTGCCTTAATg                                                                                                                          >  3:111755/1‑148 (MQ=255)
                               ttattaCGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAg                                                                                                                      >  6:89731/1‑148 (MQ=255)
                                  ataCGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTAcc                                                                                                                  >  1:70056/2‑149 (MQ=255)
                                   taCGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTAACCATTTTGTCCTGTGTTGCGTTAATTTAATTAcct                                                                                                                 >  7:38619/1‑148 (MQ=255)
                                   taCGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTAACCATTTTGTCCTGTGTTGCGTAAATGTAATTAcct                                                                                                                 >  7:42690/1‑148 (MQ=255)
                                      gTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTc                                                                                                               >  7:125685/1‑148 (MQ=255)
                                          aTCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCcaca                                                                                                           <  7:6714/148‑1 (MQ=255)
                                          aTCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAg                                                                                                          >  1:62518/1‑149 (MQ=255)
                                               gTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTgg                                                                                                     >  7:101172/1‑149 (MQ=255)
                                                ttAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTgg                                                                                                     <  8:125685/148‑1 (MQ=255)
                                                ttAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTgg                                                                                                     <  1:99870/148‑1 (MQ=255)
                                                ttAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTgg                                                                                                     <  4:111755/148‑1 (MQ=255)
                                                aataTTCATCATTTATTCATCATCATTGTATATAATTAATTTAACTTCGATAACATTAAAATATGTATCCACTGACGCTTTATAACATAACGATGAATTGACCATTTTGTCCTGTTTTGCCTTAATGTATGTACCGTCCTCTGCGTgg                                                                                                     <  6:43929/145‑1 (MQ=255)
                                                       atcatTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACt                                                                                             >  4:45454/1‑149 (MQ=255)
                                                       atcatTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCATGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACt                                                                                             >  7:243058/1‑149 (MQ=255)
                                                        tcatTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACt                                                                                             >  1:107906/1‑148 (MQ=255)
                                                          attattACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTc                                                                                           <  3:56142/148‑1 (MQ=255)
                                                          attattACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTc                                                                                           <  2:70056/148‑1 (MQ=255)
                                                           ttattaCATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTc                                                                                           >  1:240858/1‑147 (MQ=255)
                                                                 catcatcatTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCtgtt                                                                                                                                     >  4:203218/1‑99 (MQ=255)
                                                                 catcatcatTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCtgtt                                                                                                                                     <  3:203218/99‑1 (MQ=255)
                                                                 catcatcatTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAAcc                                                                                   <  3:102685/149‑1 (MQ=255)
                                                                 catcatcatTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAAc                                                                                    <  6:189016/148‑1 (MQ=255)
                                                                   tcatcatTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCtt                                                                                 <  3:177594/149‑1 (MQ=255)
                                                                        atTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTtgt                                                                             <  8:120172/148‑1 (MQ=255)
                                                                         ttGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTtgt                                                                             >  2:236422/1‑147 (MQ=255)
                                                                            tAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTc                                                                        <  2:213970/149‑1 (MQ=255)
                                                                             aaaTAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGt                                                                         >  4:196851/1‑147 (MQ=255)
                                                                              aataatTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGATGAATTGACCATTTTGTACTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAg                                                                      <  3:41128/149‑1 (MQ=255)
                                                                              aataatTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGt                                                                                                                         <  4:151538/98‑1 (MQ=255)
                                                                              aataatTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCa                                                                       <  6:229845/148‑1 (MQ=255)
                                                                              aataatTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAg                                                                      <  1:236422/149‑1 (MQ=255)
                                                                              aataatTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACAATTTTGTCCTGTTGTGCCTTAATGt                                                                                                                         >  3:151538/1‑98 (MQ=255)
                                                                                  attaaattaaCTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCaa                                                                  >  4:106564/1‑149 (MQ=255)
                                                                                   ttaaattaaCTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGAACATTTTGTCCTGTTGTGCCTGAATGTTCGTACCGTCCACAGATTGGGACATACTTCAAGGGACCTTTTGTGAGGCTGGCAAc                                                                 >  7:174577/1‑149 (MQ=255)
                                                                                      aattaaCTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAg                                                              >  7:158386/1‑149 (MQ=255)
                                                                                                caataaaTTAAAATATGGATTCACTTTAGGTTTTTTACATAAAGTATAAATAAACAATTTTGTCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGaaaa                                                    <  8:38619/143‑1 (MQ=255)
                                                                                                cATTACATTTAAATATTGTTTCACTGTCGGTTTTTTTCATAACGTACTATTAACCATTTTGGTCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGaaaa                                                    <  8:42690/149‑1 (MQ=255)
                                                                                                cATAACATTAAAAAATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATAATTCAAGGAACCATTTGTGAGACAGGCAACCAGTAAGAGaaaa                                                    >  3:191618/1‑149 (MQ=255)
                                                                                                    aCATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGa                                                 >  6:214344/1‑148 (MQ=255)
                                                                                                    aCATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGa                                                 <  2:62518/148‑1 (MQ=255)
                                                                                                       ttATAATATGTATACACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATg                                                          <  1:79750/136‑1 (MQ=255)
                                                                                                       ttAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATTTAATTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGccaca                                             >  7:146596/1‑148 (MQ=255)
                                                                                                       ttAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCAACAGACTGGGACATACTTCAAGGAACATTTTGTGAGTAAGGCAACCAGTATg                                                          >  2:79750/1‑136 (MQ=255)
                                                                                                         aaaaTATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGa                                           >  3:157503/1‑149 (MQ=255)
                                                                                                          aaaTATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGa                                           <  5:125085/148‑1 (MQ=255)
                                                                                                               tGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGaa                                                                                     <  2:215432/101‑1 (MQ=255)
                                                                                                               tGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGaa                                                                                     >  1:215432/1‑101 (MQ=255)
                                                                                                               tGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTaa                                      <  3:45454/148‑1 (MQ=255)
                                                                                                                          aCGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAATAGTCGCATGAc                          >  8:163400/1‑149 (MQ=255)
                                                                                                                              ttttttACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCg                       >  1:230550/1‑148 (MQ=255)
                                                                                                                              ttttttACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCg                       >  6:14320/1‑148 (MQ=255)
                                                                                                                                  ttAAATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGGTGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACAGATTTAACAGTCGCATGACCCGTAta                   <  8:146596/148‑2 (MQ=255)
                                                                                                                                          cGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCtc          <  3:196851/149‑1 (MQ=255)
                                                                                                                                             aagaaTTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCtcat        <  2:107906/149‑1 (MQ=255)
                                                                                                                                               gaaTTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCtcat        >  6:198133/1‑147 (MQ=255)
                                                                                                                                                  ttGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTATGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCt   >  4:52078/1‑149 (MQ=255)
                                                                                                                                                    gACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCtt  <  5:89731/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TATATCAAACAAAATAACCCTGATTAATGAATTATTACGTTTATCATGTTAATTCATCATTATTACATCATCATTGTAAATAATTAAATTAACTTCCATAACATTAAAATATGTATCCACTGACGCTTTTTTACATAACGAAGAATTGACCATTTTGTCCTGTTGTGCCTTAATGTAAGTACCGTCCACAGCGTGGGACATACTTCAAGGAACCTTTTGTGAGTCAGGCAACCAGTATGCGAAAACGACACCGATTTAACAGTCGCATGACCCGTATCGTACTGCTCATCAGCTT  >  NC_000913/2629072‑2629366

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: