Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,679,402 T→A 6.0% intergenic (+31/‑155) ydgH → / → ydgI DUF1471 family periplasmic protein/putative arginine/ornithine antiporter transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,679,4020TA 6.0% 183.3 / 2.2 67intergenic (+31/‑155)ydgH/ydgIDUF1471 family periplasmic protein/putative arginine/ornithine antiporter transporter
Reads supporting (aligned to +/- strand):  ref base T (35/28);  new base A (2/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

TCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGG  >  NC_000913/1679258‑1679550
                                                                                                                                                |                                                                                                                                                    
tcCTATCAGGTTGCGAAACGTGCCGCGATGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgccg                                                                                                                                                  >  1:253269/1‑149 (MQ=255)
tcCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgccg                                                                                                                                                  >  1:536259/1‑149 (MQ=255)
 ccTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgccgc                                                                                                                                                 <  1:93331/149‑1 (MQ=255)
  cgaTCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGCCCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgccgcc                                                                                                                                                >  2:886712/3‑149 (MQ=255)
      cAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCATCAAGCCGCCttat                                                                                                                                            <  1:12786/149‑3 (MQ=255)
       aGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCATCAAGCCGCCttata                                                                                                                                           <  2:450235/149‑4 (MQ=255)
       aGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTCGCAGGAACGTGCTAATAACCTGACCGTCAGCCCAGATCTGTATAAATAACTATTCAGGTACTCCTATAAGGCGGCTTGATGCCGCCTTTCg                                                                                                                                           >  2:762843/1‑149 (MQ=255)
          ttGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCAt                                                                                                                                        >  1:815108/1‑149 (MQ=255)
            gCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCAttt                                                                                                                                      <  2:467703/149‑1 (MQ=255)
                aaCGTGCCGCGATGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAg                                                                                                                                  <  2:885847/149‑1 (MQ=255)
                      ccGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACaa                                                                                                                            >  2:33657/1‑149 (MQ=255)
                       cgcgAAGAAATGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACaaa                                                                                                                           <  2:884930/149‑1 (MQ=255)
                       cgcgAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACaaa                                                                                                                           >  1:741758/1‑149 (MQ=255)
                          gaagaaAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAAttt                                                                                                                        <  2:715155/149‑1 (MQ=255)
                           aagaaAAGTGCTAAGTATTACCACTTCACCCGCCAGTGGCAGGAACGTTGTAATAACCTGACCGTCAGCGCAGATCTGTATTAATAACTATTCAGTTACTCCTATAAGGCGGCTTTATGCCGCCTTTCGCATTTTTAGTCACAAATTTc                                                                                                                       <  2:629532/149‑1 (MQ=255)
                                  gTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTTACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGa                                                                                                                <  2:518741/149‑1 (MQ=255)
                                      tAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCAt                                                                                                            >  2:418715/1‑149 (MQ=255)
                                         gTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCATATCTGTATAAATAACTATTCAGATACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgc                                                                                                         <  1:155937/149‑1 (MQ=255)
                                         gTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgc                                                                                                         <  2:619412/149‑1 (MQ=255)
                                          tATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgca                                                                                                        >  1:19584/1‑149 (MQ=255)
                                             tACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTCCTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCCCAAATTTCCCCTCGCCCATtgcatgc                                                                                                     >  1:1063889/1‑149 (MQ=255)
                                              agcacaTCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATATGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgcatgca                                                                                                    <  1:573391/147‑1 (MQ=255)
                                              aCCACATCTCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTAGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgcatgca                                                                                                    <  1:578103/149‑1 (MQ=255)
                                                        ccGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGCCCATTGCATGCAGACACAtctc                                                                                          >  2:1128925/1‑149 (MQ=255)
                                                            cAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCg                                                                                      >  1:448799/1‑149 (MQ=255)
                                                                  cAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAAt                                                                                >  2:249046/1‑149 (MQ=255)
                                                                          tGGTAATAACCTGCCCGTCAGCGCAGAGCTGTATAAATAACTATTCAGTTTCTCCTATAAGGCGGGTGGATGCCGCCTTTCGGAGTTTTAGTCACAAATTTCCCCGCGACCCTTGCATGCAGACCCAGCTCTTCGTAAAATCCCgcgcg                                                                        >  2:644810/1‑148 (MQ=255)
                                                                          tGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCt                                                                        <  2:725160/149‑1 (MQ=255)
                                                                           ggTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCACGCGCtt                                                                       >  1:916941/1‑149 (MQ=255)
                                                                             taataaCCTGACCGTCAGCGCAGATGTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGCTGCAGACTTTCGCATTTTTAGTCACAAATTTCCCCGCGTCCAGTGCAGGCAGAAACATCTCGTCGTAAAATGCCGCGCTTCa                                                                     >  2:776102/1‑149 (MQ=255)
                                                                               ataaCCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCTCAAATTTCCCCGCGACCATTGCTTGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGt                                                                   <  1:886712/149‑1 (MQ=255)
                                                                               ataaCCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGt                                                                   <  2:10452/149‑1 (MQ=255)
                                                                                      gACCGTCTGAGCATATCTGTATTAAAAACTATTCAGCTTCTCCAATCTGGCGGCATTATGCGGCCTTTCGCATTTATAGTCACAAATTTCCCCGCGACCATAGCCTGCAGACACTTCTCTTCGTTATATCCCTCGCTTCAGTGGCattc                                                            <  2:515234/149‑4 (MQ=255)
                                                                                         cGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTcc                                                         <  2:981243/149‑1 (MQ=255)
                                                                                         cGTAAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTcc                                                         >  2:731602/1‑149 (MQ=255)
                                                                                          gTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTACa                                                        >  2:613106/1‑149 (MQ=255)
                                                                                               cgcAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCAttttt                                                   >  1:413270/1‑149 (MQ=255)
                                                                                               cgaaGAACTGTATAAATAACTATTCAGTTACTACTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGAAACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCAttttt                                                   >  2:567559/4‑149 (MQ=255)
                                                                                                 cAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTAt                                                 <  1:731602/149‑1 (MQ=255)
                                                                                                  aGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCTAACTCCATTTTTATg                                                >  1:262185/1‑149 (MQ=255)
                                                                                                         tataAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTg                                         <  1:418715/149‑1 (MQ=255)
                                                                                                          ataaataaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGc                                        >  1:1145055/1‑149 (MQ=255)
                                                                                                          ataaataaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGc                                        <  1:638876/149‑1 (MQ=255)
                                                                                                                aaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGAAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGGGTTTTGCCATatt                                  >  1:590851/1‑147 (MQ=255)
                                                                                                                   tATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGCCCCAGCGCTTCGTAAAAGCCCGCGCTGCCGTGGCAACCTCCATTTTGATGCGGTTTGACCTAAtatt                               >  1:769163/1‑149 (MQ=255)
                                                                                                                   tATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtatt                               >  2:561635/1‑149 (MQ=255)
                                                                                                                    aTTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtattt                              <  1:51759/149‑1 (MQ=255)
                                                                                                                     ttCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGCCCCTTGCCTGGAGACACATCACTTCGTAAAATCCCGCGCTTCCGTGGCAGACTCCTTATTTATGCGTTATGCCCTAAtatttt                             >  2:1094582/1‑148 (MQ=255)
                                                                                                                     ttCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTATAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtattta                             >  2:477375/1‑149 (MQ=255)
                                                                                                                        aGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAATCTCCATTTTTATGCGTTTTGCCCTAAtatttatta                          >  2:1111956/1‑148 (MQ=255)
                                                                                                                        aGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtatttattg                          >  1:83685/1‑148 (MQ=255)
                                                                                                                          ttACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTAttctt                        >  2:930284/1‑149 (MQ=255)
                                                                                                                             cTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTcttctt                     <  1:249046/149‑1 (MQ=255)
                                                                                                                                cTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATc                  <  2:336078/149‑1 (MQ=255)
                                                                                                                                  gaaaGGCGGCATCAAGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCAc                >  1:298299/3‑149 (MQ=255)
                                                                                                                                    aaGGCGGCTTGATGCCGCCTTTCGCAATTTTAGTCACAAATTTCCCCGCGACAATTGCCTGCAGAAACATATCTTCGTAAAATACAGCTATTCAGTGGCAAACTCGATTTTTATGCGTTTTGCCCTAATATTTATTGTTCTTATCACGt              >  2:202452/1‑149 (MQ=255)
                                                                                                                                    aaGGCGGCTTCATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACAACTCTTCGTAAAATCCCGCGCTTCCGTGGCAAACTCCATTTTTATGCGTTTTGCACTAATATTTAATCTGCTTATCACGt              >  2:772227/1‑149 (MQ=255)
                                                                                                                                     aGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCGCAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGtt             <  2:483530/149‑1 (MQ=255)
                                                                                                                                     aGGCGGCATCAAGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCAATGCATGCAGACACATCTCGTCGTAAAATCCCGCGATTCAGTGGCAAACGCCATATTTATGCGTTTTGCCCTAATATTTATTATTGTTATCACGtt             >  2:778142/1‑149 (MQ=255)
                                                                                                                                       gcggcTTGATGCCGCCTTTCGGATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCTGTGGCAAACTCCATTTATATGCGTTTTGCCCTAATTTTTATTCTTCTTATCACCtttt           >  2:930181/1‑149 (MQ=255)
                                                                                                                                       gcggcTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATGACGtttt           >  1:217521/1‑149 (MQ=255)
                                                                                                                                       gcggcTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGtttt           <  2:253269/149‑1 (MQ=255)
                                                                                                                                        gggcTTGATGCCGCCTTTCGGATTTTTAGTCACAAATTTCCACTCGACCATTGCATGCAGACAAATCTCTTCGTAAAATCCCGCGCTTCAGTCGCAAACTCCATTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTaa         <  1:778316/148‑1 (MQ=255)
                                                                                                                                          gcTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAAt        <  2:741758/149‑1 (MQ=255)
                                                                                                                                           cgtGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATc       <  2:815108/147‑1 (MQ=255)
                                                                                                                                           cTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATc       >  1:1170729/1‑149 (MQ=255)
                                                                                                                                              gATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACt    <  1:713213/149‑1 (MQ=255)
                                                                                                                                                tGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTgg  <  2:536511/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGG  >  NC_000913/1679258‑1679550

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: