Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,679,400 G→C 7.2% intergenic (+29/‑157) 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,4000GC 7.2% 177.8 / 4.8 66intergenic (+29/‑157)ydgH/ydgIDUF1471 family periplasmic protein/putative arginine/ornithine antiporter transporter
Reads supporting (aligned to +/- strand):  ref base G (34/25);  new base C (3/2);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACT  >  NC_000913/1679258‑1679548
                                                                                                                                              |                                                                                                                                                    
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)
                                                                                                                                              |                                                                                                                                                    
TCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACT  >  NC_000913/1679258‑1679548

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: