Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,187,314 C→G 12.3% intergenic (+16/+66) rcnB → / ← yehA periplasmic modulator of Ni and Co efflux/putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,187,3140CG12.3% 204.7 / 19.1 65intergenic (+16/+66)rcnB/yehAperiplasmic modulator of Ni and Co efflux/putative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base C (25/32);  new base G (5/3);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

GTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACT  >  NC_000913/2187169‑2187460
                                                                                                                                                 |                                                                                                                                                  
gTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATg                                                                                                                                                 >  2:536141/1‑149 (MQ=255)
   gCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAgg                                                                                                                                              <  4:424619/149‑1 (MQ=255)
     gTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAgggg                                                                                                                                            >  3:134450/1‑149 (MQ=255)
     gTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAgggg                                                                                                                                            >  2:578215/1‑149 (MQ=255)
              ccGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATGATGAggggg                                                                                                                                           <  2:86719/141‑1 (MQ=255)
              ccGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATGATGAggggg                                                                                                                                           >  1:86719/1‑141 (MQ=255)
              ccGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAg                                                                                                                                   <  3:439178/149‑1 (MQ=255)
                  gcCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATCAAAGCCTACGACGGTGAGATTTCTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGacac                                                                                                                               <  2:204151/149‑1 (MQ=255)
                           cGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACAGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTAtt                                                                                                                      >  2:570982/1‑149 (MQ=255)
                                cGCACTGGACCTATATGGGTGGCGCGTTCGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTtta                                                                                                                 <  3:475297/149‑1 (MQ=255)
                                cGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTtta                                                                                                                 <  1:578215/149‑1 (MQ=255)
                                cGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTtta                                                                                                                 <  2:602849/149‑1 (MQ=255)
                                   aCTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAg                                                                                                                                               <  4:525588/116‑1 (MQ=255)
                                   aCTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAg                                                                                                                                               >  3:525588/1‑116 (MQ=255)
                                       gACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTttattatta                                                                                                           <  1:296985/148‑1 (MQ=255)
                                       gACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTttattatta                                                                                                           >  2:296985/1‑148 (MQ=255)
                                          cTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCa                                                                                                       <  2:577351/149‑1 (MQ=255)
                                          cTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCa                                                                                                       >  4:66594/1‑149 (MQ=255)
                                          cTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCa                                                                                                       >  4:23280/1‑149 (MQ=255)
                                           tataTGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCAc                                                                                                      <  3:331949/149‑1 (MQ=255)
                                           tataTGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCAc                                                                                                      <  1:389711/149‑1 (MQ=255)
                                            atatGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTAttttt                                                                                                                   >  3:591153/1‑135 (MQ=255)
                                            atatGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTAttttt                                                                                                                   <  4:591153/135‑1 (MQ=255)
                                             tatGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGc                                                                                                                                     >  2:463697/1‑116 (MQ=255)
                                             tatGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGc                                                                                                                                     <  1:463697/116‑1 (MQ=255)
                                               tGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGc                                                                                                  <  4:176044/149‑1 (MQ=255)
                                                    ggCGCGTACGTGTTGATCAGCGACACCGACAGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTc                                                                                             <  1:570982/149‑1 (MQ=255)
                                                            cGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGt                                                                                     <  3:391844/149‑1 (MQ=255)
                                                              tgtTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATGATGAgggggc                                                                                                                                          >  2:409071/1‑93 (MQ=255)
                                                              tgtTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATGATGAgggggc                                                                                                                                          <  1:409071/94‑2 (MQ=255)
                                                              tgtTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGtta                                                                                   >  3:567151/1‑149 (MQ=255)
                                                                 tGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGttatta                                                                                <  3:66594/149‑1 (MQ=255)
                                                                  gATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAg                                                                               >  2:22810/1‑149 (MQ=255)
                                                                   aTCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGt                                                                              >  3:675996/1‑149 (MQ=255)
                                                                   aTCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGt                                                                              <  1:465377/149‑1 (MQ=255)
                                                                    tCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTc                                                                             <  3:23280/149‑1 (MQ=255)
                                                                    tCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTc                                                                             >  3:492340/1‑149 (MQ=255)
                                                                       gCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTa                                                                          <  4:257233/149‑1 (MQ=255)
                                                                        cGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTAt                                                                         >  2:261091/1‑149 (MQ=255)
                                                                             ccGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCAc                                                                    <  3:46145/149‑1 (MQ=255)
                                                                              cGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGcc                                                                                                        <  3:330753/112‑1 (MQ=255)
                                                                              cGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGcc                                                                                                        >  4:330753/1‑112 (MQ=255)
                                                                              cGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACg                                                                   >  3:198745/1‑149 (MQ=255)
                                                                               gACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGt                                                                  >  3:347482/1‑149 (MQ=255)
                                                                                aCGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGtt                                                                 >  2:467966/1‑149 (MQ=255)
                                                                                    tAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGAtt                                                             >  4:135182/1‑149 (MQ=255)
                                                                                    tAAAATCATTAAAGCCTACGACGGGGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGAtt                                                             >  2:131428/1‑149 (MQ=255)
                                                                                          cATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCg                                                       >  3:397060/1‑149 (MQ=255)
                                                                                          cATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCg                                                       <  3:135182/149‑1 (MQ=255)
                                                                                              aaaGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGa                                                               >  1:412074/1‑137 (MQ=255)
                                                                                              aaaGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGa                                                               <  2:412074/137‑1 (MQ=255)
                                                                                                           gTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGcc                                      <  4:198745/149‑1 (MQ=255)
                                                                                                                   ttttATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTaa                              <  4:397060/149‑1 (MQ=255)
                                                                                                                       atcatcGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACttt                          <  4:567151/149‑1 (MQ=255)
                                                                                                                        tcatcGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACtttg                         <  1:601528/149‑1 (MQ=255)
                                                                                                                         catcGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACtttgt                        <  1:22810/149‑1 (MQ=255)
                                                                                                                          atcGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACtttgtt                       <  1:261091/149‑1 (MQ=255)
                                                                                                                           tcGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACtttgttt                      <  4:11513/149‑1 (MQ=255)
                                                                                                                           tcGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACtttgttt                      >  3:11513/1‑149 (MQ=255)
                                                                                                                                   aaaaGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTa              >  4:131068/1‑149 (MQ=255)
                                                                                                                                    atttcccccTCATGATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAg             >  1:549933/5‑149 (MQ=255)
                                                                                                                                     tttcccccTCATGATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGa            >  1:71065/4‑149 (MQ=255)
                                                                                                                                     aaGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGa            <  2:98964/149‑1 (MQ=255)
                                                                                                                                        cccccTCATGATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGc                                                  >  4:485872/1‑108 (MQ=255)
                                                                                                                                        cccccTCATGATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGc                                                  <  3:485872/108‑1 (MQ=255)
                                                                                                                                             tcatcaTGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTaa    <  3:486422/149‑1 (MQ=255)
                                                                                                                                               atcatGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACt  >  4:411392/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACT  >  NC_000913/2187169‑2187460

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