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

AGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATC  >  NC_000913/1679254‑1679545
                                                                                                                                                  |                                                                                                                                                 
aGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGa                                                                                                                                                  <  1:674638/148‑1 (MQ=255)
  tctcCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgc                                                                                                                                               >  3:234338/1‑149 (MQ=255)
   ctcCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgcc                                                                                                                                              <  4:365727/149‑1 (MQ=255)
      cTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgccgcc                                                                                                                                           >  3:446709/1‑149 (MQ=255)
       tATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCt                                                                                                                                          >  3:140561/1‑149 (MQ=255)
       tATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCATCAAGCCGCCt                                                                                                                                          <  4:724725/149‑1 (MQ=255)
       tATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCATCAAGCCGCCt                                                                                                                                          >  4:631475/1‑149 (MQ=255)
            ggTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCATCAAGCCGCCttatag                                                                                                                                     <  3:631475/149‑5 (MQ=255)
                   aaaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCAttttt                                                                                                                               <  2:514088/148‑1 (MQ=255)
                   aaaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTa                                                                                                                              <  4:140561/149‑1 (MQ=255)
                   aaaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTa                                                                                                                              <  4:234338/149‑1 (MQ=255)
                       gTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGAca                                                                                                                          <  3:599801/149‑1 (MQ=255)
                           cgcgAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCGGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACaaa                                                                                                                      <  3:771660/149‑1 (MQ=255)
                             cgAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAAt                                                                                                                     >  2:297059/1‑148 (MQ=255)
                                agaaAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTcc                                                                                                                 <  3:402233/149‑1 (MQ=255)
                                agaaAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTc                                                                                                                  <  1:862301/148‑1 (MQ=255)
                                           aaGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCAt                                                                                                       >  1:209723/1‑148 (MQ=255)
                                            aGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCAtt                                                                                                      >  1:102484/1‑148 (MQ=255)
                                             gTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtg                                                                                                     <  1:462132/148‑1 (MQ=255)
                                                ttACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgcatg                                                                                                 >  3:258802/1‑149 (MQ=255)
                                                  aCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgcatgca                                                                                               <  4:446709/149‑1 (MQ=255)
                                                  aCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgcatgca                                                                                               >  3:618630/1‑149 (MQ=255)
                                                   ccACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgcatgca                                                                                               <  2:209723/148‑1 (MQ=255)
                                                         cACCCGTCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGAcaca                                                                                         <  2:634743/148‑1 (MQ=255)
                                                          aCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACAtc                                                                                       >  3:601445/1‑149 (MQ=255)
                                                            ccGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACAtct                                                                                      >  2:424145/1‑148 (MQ=255)
                                                               ccAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCtt                                                                                   <  1:488141/148‑1 (MQ=255)
                                                                   tGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTaa                                                                              >  4:518372/1‑149 (MQ=255)
                                                                     gCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTaaa                                                                             <  2:282723/148‑1 (MQ=255)
                                                                                gTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTc                                                                 <  4:618630/149‑1 (MQ=255)
                                                                                 taataaCCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGAc                                                                                                          <  1:314908/107‑1 (MQ=255)
                                                                                 taataaCCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGAc                                                                                                          >  2:314908/1‑107 (MQ=255)
                                                                                 taataaCCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTc                                                                 >  2:724031/1‑148 (MQ=255)
                                                                                  aataaCCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTc                                                                                                                  <  2:607519/98‑1 (MQ=255)
                                                                                  aataaCCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTc                                                                                                                  >  1:607519/1‑98 (MQ=255)
                                                                                  aataaCCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAg                                                               >  4:714649/1‑149 (MQ=255)
                                                                                   ataaCCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGt                                                              <  4:800258/149‑1 (MQ=255)
                                                                                    taaCCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGt                                                              <  1:592535/148‑1 (MQ=255)
                                                                                        cTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGGTACTCCTATAATGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCATTGGCaa                                                         >  4:194890/1‑149 (MQ=255)
                                                                                         tGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTcc                                                                                                                 >  2:352417/1‑92 (MQ=255)
                                                                                         tGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTcc                                                                                                                 <  1:352417/92‑1 (MQ=255)
                                                                                              gTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTcc                                                    <  2:548018/148‑1 (MQ=255)
                                                                                                         tCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCg                                         >  1:524503/1‑148 (MQ=255)
                                                                                                              ataaataaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTg                                    <  2:490574/148‑1 (MQ=255)
                                                                                                              ataaataaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGc                                   <  4:601445/149‑1 (MQ=255)
                                                                                                               taaataaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGcc                                  <  3:194890/149‑1 (MQ=255)
                                                                                                                aaataaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGccc                                 >  3:110005/1‑149 (MQ=255)
                                                                                                                   taaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGt                                                              >  2:805448/1‑117 (MQ=255)
                                                                                                                   taaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGt                                                              <  1:805448/117‑1 (MQ=255)
                                                                                                                       tATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAatat                           <  1:424145/148‑1 (MQ=255)
                                                                                                                         ttCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtattt                         <  1:221966/148‑1 (MQ=255)
                                                                                                                          tCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtatttat                       <  3:714649/149‑1 (MQ=255)
                                                                                                                               tACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTAttcttc                  >  3:288627/1‑149 (MQ=255)
                                                                                                                               tACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTAttcttc                  <  4:324234/149‑1 (MQ=255)
                                                                                                                                aCTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTaa                              <  3:114780/136‑1 (MQ=255)
                                                                                                                                aCTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTaa                              >  4:114780/1‑136 (MQ=255)
                                                                                                                                    cTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATc             >  4:163690/1‑149 (MQ=255)
                                                                                                                                      ataAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtatttatt                      >  1:621345/1‑138 (MQ=255)
                                                                                                                                      ataAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtatttatt                      <  2:621345/138‑1 (MQ=255)
                                                                                                                                       taAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCAc           >  2:613970/1‑148 (MQ=255)
                                                                                                                                       taAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCAc           >  1:290085/1‑148 (MQ=255)
                                                                                                                                       aaaGGCGGCATCAAGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCAc           >  2:243847/2‑148 (MQ=255)
                                                                                                                                             ggcTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTaa    <  3:301188/149‑1 (MQ=255)
                                                                                                                                                ttGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATc  <  1:297059/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATC  >  NC_000913/1679254‑1679545

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: