Predicted mutation
evidence seq id position mutation annotation gene description
RA 4 23,302 (G)5→4 intergenic (‑53/+122) unknown ← / ← unknown Ferric siderophore transport system, biopolymer transport protein ExbB/Putative membrane protein YfjD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*423,2980G.100.0% 158.9 / NA 39intergenic (‑49/+126)unknown/unknownFerric siderophore transport system, biopolymer transport protein ExbB/Putative membrane protein YfjD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (20/19);  total (20/19)

CACAGTAGAAACTATGGTGAGGTTACGGTTTAAATCAATTTCAAGTGCATAAATGGTGGAATATTGACTCACTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGAGGGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCATATTCTCACTCACTTCCAACACGGTAACCTGCAAGCCATTAAGTTCACAGACTGTGCCTTCGTGTGGGATTTCTTCCAAATG  >  4/23052‑23537
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
cacaGTAGAAACTATGGTGAGGTTACGGTTTAAATCAATTTCAAGTGCATAAATGGTGGAATATTGACTCACTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGAggggc                                                                                                                                                                                                                                             <  1:218906/251‑2 (MQ=255)
  caGTAGAAACTATGGTGAGGTTACGGTTTAAATCAATTTCAAGTGCATAAATGGTGGAATATTGACTCACTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGAggggcag                                                                                                                                                                                                                                           >  1:78552/1‑248 (MQ=255)
    gTAGAAACTATGGTGAGGTTACGGTTTAAATCAATTTCAAGTGCATAAATGGTGGAATATTGACTCACTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGAggggc                                                                                                                                                                                                                                             <  1:691392/247‑2 (MQ=255)
                gTGAGGTTACGGTTTAAATCAATTTCAAGTGCATAAATGGTGGAATATTGACTCACTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTc                                                                                                                                                                                                                             <  1:237099/250‑1 (MQ=255)
                    ggTTACGGTTTAAATCAATTTCAAGTGCATAAATGGTGGAATATTGACTCACTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATAcc                                                                                                                                                                                                                                  >  1:214205/1‑241 (MQ=255)
                      ttACGGTTTAAATCAATTTCAAGTGCATAAATGGTGGAATATTGACTCACTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAAc                                                                                                                                                                                                                       <  1:756445/250‑1 (MQ=255)
                                           cgccattcaactgggctggatTGACTCACTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTg                                                                                                                                                                                                 >  1:623987/20‑251 (MQ=255)
                                            gTGCATAAATGGTGGAATATTGACTCACTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGt                                                                                                                                                                                                >  1:868911/1‑251 (MQ=255)
                                            gTGCATAAATGGTGGAATATTGACTCACTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAAAGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTAt                                                                                                                                                                                                   >  1:741227/1‑248 (MQ=255)
                                              gCATAAATGGTGGAATATTGACTCACTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGAggggc                                                                                                                                                                                                                                             <  1:890964/205‑2 (MQ=255)
                                                      ggtggAATATTGACTCACTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAAc                                                                                                                                                                                                                       >  1:571832/1‑218 (MQ=255)
                                                                       cTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATc                                                                                                                                                                                       <  1:146336/233‑1 (MQ=255)
                                                                                   aTAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATAcc                                                                                                                                                                                                                                  <  1:877949/178‑1 (MQ=255)
                                                                                                aCGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGt                                                                                                                                            <  1:203106/251‑1 (MQ=255)
                                                                                                   ttCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGttt                                                                                                                                          >  1:898357/1‑250 (MQ=255)
                                                                                                                    gCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATAccc                                                                                                                                                                                                                                 <  1:215197/146‑1 (MQ=255)
                                                                                                                     cATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATAccc                                                                                                                                                                                                                                 <  1:268385/145‑1 (MQ=255)
                                                                                                                       ttAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAAtt                                                                                                                     >  1:646879/1‑251 (MQ=255)
                                                                                                                               cTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATaa                                                                                                              <  1:95465/250‑1 (MQ=255)
                                                                                                                               cTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACATGATCGTTTCTGATTACGTAAACACAAATTc                                                                                                                    >  1:185241/1‑244 (MQ=255)
                                                                                                                                       cAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAAc                                                                                                                                                                                                                       <  1:272275/137‑1 (MQ=255)
                                                                                                                                                  cGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAAtt                                                                                                            <  1:793926/233‑1 (MQ=255)
                                                                                                                                                     taataaTGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGAggggcag                                                                                                                                                                                                                                           <  1:492400/104‑4 (MQ=255)
                                                                                                                                                       ataatGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAAt                                                                                      >  1:543910/1‑250 (MQ=255)
                                                                                                                                                         aatGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCat                                                                                   <  1:692819/251‑1 (MQ=255)
                                                                                                                                                                   cTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTc                                                                                                                                         >  1:99410/1‑187 (MQ=255)
                                                                                                                                                                   cTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGCGATCTGTTTAAATAAATAGCAGAAAGA‑GGTGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCatat                                                                                 >  1:591211/1‑243 (MQ=255)
                                                                                                                                                                                 gAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGAtt                                                                                                                                    <  1:548192/178‑1 (MQ=255)
                                                                                                                                                                                   aTTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCATATTCTCACTCACTTCCAACACGGTaa                                                         >  1:284816/1‑251 (MQ=255)
                                                                                                                                                                                             ttGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCATATTCTCACTCACTTCCAACACGGTAACCTGCAAGcc                                               <  1:408668/251‑1 (MQ=255)
                                                                                                                                                                                                        gCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCAGCTTAATCATATTCTCACTCACtt                                                                     >  1:311455/1‑218 (MQ=255)
                                                                                                                                                                                                           tttCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCATATTCTCACTCACTTCCAACACGGT‑ACCTGCAAGCCATTAAGTTCACAGAc                                <  1:317779/251‑1 (MQ=255)
                                                                                                                                                                                                            ttCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCATATTCTCACTCACTTCCAACACGGTAACCTGCAAGCCATTAAGTTCACAGAc                                <  1:807028/251‑1 (MQ=255)
                                                                                                                                                                                                            ttCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCATATTCTCACTCACTTCCAACACGGTAACCTGCAAGCCATTAAGTTCACAGAc                                <  1:56913/251‑1 (MQ=255)
                                                                                                                                                                                                                      aaTGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATAcc                                                                                                                                                                                                                                  >  1:412828/1‑47 (MQ=255)
                                                                                                                                                                                                                      aaTGTGATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTaaaaaaa                                                                                                                                                                                         >  1:881212/1‑88 (MQ=255)
                                                                                                                                                                                                                         gtgATCTGTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAacac                                                                                                                                                    <  1:370118/122‑1 (MQ=255)
                                                                                                                                                                                                                                gTTTAAATAAATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCATATTCTCACTCACTTCCAACACGGTAACCTGCAAGCCATTAAGTTCACAGACTGTGCCTTCGTGTGGGAttt            >  1:848202/1‑251 (MQ=255)
                                                                                                                                                                                                                                gTTTAAAT‑AATAGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCATATTCTCACTCACTTCCAACACGGTAACCTGCAAGCCATTAAGTTCACAGACTGTGCCTTCGTGTGGGAttt            >  1:352056/1‑250 (MQ=255)
                                                                                                                                                                                                                                        aaataGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAAcc                                                                                                    <  1:384348/155‑1 (MQ=255)
                                                                                                                                                                                                                                          ataGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCATATTCTCACTCACTTCCAACACGGTAACCTGc                                                    >  1:585788/1‑201 (MQ=255)
                                                                                                                                                                                                                                          ataGTAGAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCATATTCTCACTCACTTCCAACACGGTAACCTGCAAGCCATTAAGTTCACAGACTGTGCCTTCGTGTGGGATTTCTTCCAAATg  >  1:405443/1‑251 (MQ=255)
                                                                                                                                                                                                                                           tagtagAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCATCTTCTCACTCACTTCCAACACGGTAACCTGCAAGCCATTAAGTTCACAGACtgtgggagcaaacgcgccctatgtggg     >  1:454847/1‑224 (MQ=255)
                                                                                                                                                                                                                                             gtagAAAGA‑GGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCATATtct                                                                              <  1:719882/172‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
CACAGTAGAAACTATGGTGAGGTTACGGTTTAAATCAATTTCAAGTGCATAAATGGTGGAATATTGACTCACTTTAACACGGGATAAAAAGAAGAAACGTTCGATGACCATGGTGAGCATTAATACGCTCATAAACAGTAATAAACCGATAATAATGTAGTCACTATAATTTTGTAAGAATTCAAAAATTTGTGGCATTTGCTTTTCCTTTTTAAATGTGATCTGTTTAAATAAATAGTAGAAAGAGGGGGCAGATCATACCCCTTCCATAACTGATAATCAATCTCATTTATTGTAAAAAAATCTGCGCTTCTTTTTTTAAAAAAAGCGATAAAAACACTTGATCGTTTCTGATTACGTAAACACAAATTCTTATAATTTAATAACCTTCGCCCGCTTAATCATATTCTCACTCACTTCCAACACGGTAACCTGCAAGCCATTAAGTTCACAGACTGTGCCTTCGTGTGGGATTTCTTCCAAATG  >  4/23052‑23537

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