Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,341,631 T→C T847A (ACT→GCT)  yfaL ← adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,341,6310TC100.0% 240.6 / NA 61T847A (ACT→GCT) yfaLadhesin
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (40/21);  total (40/21)

ACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCATTTTGCAGGCTAACCGCACCGGA  >  NC_000913/2341486‑2341775
                                                                                                                                                 |                                                                                                                                                
aCGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTAc                                                                                                                                             >  2:6005/1‑151 (MQ=255)
 cGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTAcc                                                                                                                                            >  2:7551/1‑151 (MQ=255)
  gTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCg                                                                                                                                           >  2:564071/1‑151 (MQ=255)
   tAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGt                                                                                                                                          >  2:710910/1‑151 (MQ=255)
    aGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGtt                                                                                                                                         >  1:842206/1‑151 (MQ=255)
           gcctgccAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCAcc                                                                                                                                  >  1:672176/1‑151 (MQ=255)
               gccAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAAtt                                                                                                                              >  2:88996/1‑151 (MQ=255)
               gccAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAAt                                                                                                                               <  2:175283/150‑1 (MQ=255)
                  aGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGAt                                                                                                                           <  2:55512/151‑1 (MQ=255)
                  aGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGAt                                                                                                                           <  2:1012623/151‑1 (MQ=255)
                       gAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCtt                                                                                                                      <  1:646352/151‑1 (MQ=255)
                                aTTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTc                                                                                                              >  2:354058/1‑150 (MQ=255)
                                  tGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATc                                                                                                           >  1:539782/1‑151 (MQ=255)
                                   gTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCg                                                                                                          <  2:166665/151‑1 (MQ=255)
                                             tGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACCGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAACTGATCGCTTACCGAGTCATCGCCGTTTAAtt                                                                                                >  2:1156390/1‑151 (MQ=255)
                                                gTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGc                                                                                             >  1:513188/1‑151 (MQ=255)
                                                  gggATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCt                                                                                            <  2:546615/150‑1 (MQ=255)
                                                      tCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCg                                                                                        >  2:90339/1‑150 (MQ=255)
                                                          cTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGc                                                                                     >  2:654665/1‑149 (MQ=255)
                                                           tGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGt                                                                                  >  1:1011365/1‑151 (MQ=255)
                                                              gAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGc                                                                               >  1:383401/1‑151 (MQ=255)
                                                                  tCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACATGTGTGTTGCCAGCATCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCGATCGAGCAGTAGCGGAc                                                                           >  1:522379/1‑151 (MQ=255)
                                                                  tCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTAc                                                                           >  2:29056/1‑151 (MQ=255)
                                                                    aaCCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTAccg                                                                         >  1:802886/1‑151 (MQ=255)
                                                                    aaCCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTAccg                                                                         <  1:923877/151‑1 (MQ=255)
                                                                     aCCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTAccgc                                                                        <  1:283827/151‑1 (MQ=255)
                                                                         cTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTAccgccgcc                                                                    >  1:1192932/1‑151 (MQ=255)
                                                                         cTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTAccgccgc                                                                     <  1:670138/150‑1 (MQ=255)
                                                                           ttAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTAccgccgcc                                                                    <  2:676433/149‑1 (MQ=255)
                                                                            tAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTAccgccgccga                                                                  <  1:1078728/150‑2 (MQ=255)
                                                                             aaTGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGt                                                                 <  1:939213/150‑1 (MQ=255)
                                                                              aTGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGtat                                                               >  1:496608/1‑151 (MQ=255)
                                                                                        aTGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGtt                                                     >  1:247248/1‑151 (MQ=255)
                                                                                         tGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGttt                                                    >  1:939854/1‑151 (MQ=255)
                                                                                         tGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCACTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGttt                                                    >  2:1123306/1‑151 (MQ=255)
                                                                                          gTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGttt                                                    >  1:615907/1‑150 (MQ=255)
                                                                                              gCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGt                                               >  2:215071/1‑151 (MQ=255)
                                                                                              gCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGt                                               >  2:443798/1‑151 (MQ=255)
                                                                                                  aCCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGc                                           <  2:625164/151‑1 (MQ=255)
                                                                                                    cAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCg                                          >  1:759059/1‑150 (MQ=255)
                                                                                                    cAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGt                                         <  2:437743/151‑1 (MQ=255)
                                                                                                     aaTCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGt                                         >  2:574865/1‑150 (MQ=255)
                                                                                                            gTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGcc                                   >  2:950574/1‑149 (MQ=255)
                                                                                                                  gAGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGc                           >  2:305783/1‑151 (MQ=255)
                                                                                                                   aGTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGcc                          >  2:928948/1‑151 (MQ=255)
                                                                                                                    gTTAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCACCACCAATTTATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCTCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCa                         <  2:509116/151‑1 (MQ=255)
                                                                                                                     ttAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCAt                        >  2:26839/1‑151 (MQ=255)
                                                                                                                     ttAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCAt                        >  2:1121000/1‑151 (MQ=255)
                                                                                                                     ttAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCAt                        >  2:1133091/1‑151 (MQ=255)
                                                                                                                      tAACCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCAtt                       >  1:880465/1‑151 (MQ=255)
                                                                                                                       aaCCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCAttt                      >  1:606487/1‑151 (MQ=255)
                                                                                                                       aaCCACCACAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCAttt                      >  2:1113206/1‑151 (MQ=255)
                                                                                                                          caccacAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCATTTTGc                   >  2:970898/1‑151 (MQ=255)
                                                                                                                            ccacAGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCATTTTGCAg                 <  1:504031/151‑1 (MQ=255)
                                                                                                                                aGTTGTGTTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCATTTTGCAGGCTa             <  1:6005/151‑1 (MQ=255)
                                                                                                                                    gtgtTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCATTTTGCAGGCTAACCg         <  2:853234/151‑1 (MQ=255)
                                                                                                                                     tgtTGCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCATTTTGCAGGCTAACCGc        >  2:1022677/1‑151 (MQ=255)
                                                                                                                                         gCCAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCATTTTGCAGGCTAACCGCACCg    <  1:272039/151‑1 (MQ=255)
                                                                                                                                          ccAGCAGCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTTGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGGATCGCCAGCGACGCCATTATGCAGGCTACCCGCACCg    <  1:582283/150‑1 (MQ=255)
                                                                                                                                           cagcagCATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTCCCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCATTTTGCAGGCTAACCGCACCGGa  <  2:288280/151‑1 (MQ=255)
                                                                                                                                           caacagCATTACCGTTCATCACCAATTGAGCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCATTTTGCAGGCTAACCGCACCgg   >  2:164989/4‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCATCGCCGTTTAATTCGCTATCGAGCAGTAGCGTACCGCCGCCGGTATAATCACCGTTTACCGTCAGCGTATCGCCAGCGACGCCATTTTGCAGGCTAACCGCACCGGA  >  NC_000913/2341486‑2341775

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