Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,679,398 T→A 7.3% intergenic (+27/‑159) 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,3980TA 7.3% 171.2 / 3.2 65intergenic (+27/‑159)ydgH/ydgIDUF1471 family periplasmic protein/putative arginine/ornithine antiporter transporter
Reads supporting (aligned to +/- strand):  ref base T (35/25);  new base A (2/3);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATC  >  NC_000913/1679258‑1679545
                                                                                                                                            |                                                                                                                                                   
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)
                                                                                                                                           cTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATc  >  1:1170729/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
TCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATC  >  NC_000913/1679258‑1679545

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: