Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 562,986 T→C G217G (GGT→GGC)  sfmH → FimA homolog, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913562,9860TC100.0% 270.9 / NA 76G217G (GGT→GGC) sfmHFimA homolog, function unknown
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (37/39);  total (37/39)

GTTACACCTGTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAAT  >  NC_000913/562852‑563118
                                                                                                                                      |                                                                                                                                    
gTTACACCTGTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATc                                                                                                                                  >  1:1199006/1‑139 (MQ=255)
 ttACACCTGTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCg                                                                                                                                 <  2:306984/139‑1 (MQ=255)
   acacCTGTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAc                                                                                                                               >  1:611017/1‑139 (MQ=255)
   acacCTGTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAc                                                                                                                               <  2:1199006/139‑1 (MQ=255)
       cTGTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAc                                                                                                                               >  1:1244076/1‑135 (MQ=255)
       cTGTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAc                                                                                                                               <  2:1244076/135‑1 (MQ=255)
          tttATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAcc                                                                                                                              <  1:300844/133‑1 (MQ=255)
          tttATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAcc                                                                                                                              >  2:300844/1‑133 (MQ=255)
          tttATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGTTgg                                                                                                                        >  2:5067/1‑139 (MQ=255)
          tttATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTgg                                                                                                                        >  2:1177609/1‑139 (MQ=255)
               aCCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCa                                                                                                                   <  1:1177609/139‑1 (MQ=255)
                   ttAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACctgc                                                                                                                           <  2:866308/127‑1 (MQ=255)
                   ttAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACctgc                                                                                                                           >  1:866308/1‑127 (MQ=255)
                   ttAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGcc                                                                                                               >  2:876958/1‑139 (MQ=255)
                    tAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAc                                                                                                                               <  1:260558/122‑1 (MQ=255)
                    tAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAc                                                                                                                               >  2:260558/1‑122 (MQ=255)
                    tAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCg                                                                                                              >  1:178867/1‑139 (MQ=255)
                    tAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCg                                                                                                              >  2:94664/1‑139 (MQ=255)
                           ggtggCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGGTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCa                                                                                                       >  1:211997/1‑139 (MQ=255)
                            gtggCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAc                                                                                                                               >  2:1228124/1‑114 (MQ=255)
                            gtggCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAc                                                                                                                               <  1:1228124/114‑1 (MQ=255)
                                 cGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAc                                                                                                                               >  2:761479/1‑109 (MQ=255)
                                 cGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAc                                                                                                                               <  1:761479/109‑1 (MQ=255)
                                              cGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTc                                                                                         <  1:534715/134‑1 (MQ=255)
                                              cGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTc                                                                                         >  2:534715/1‑134 (MQ=255)
                                              cGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATg                                                                                    >  2:87780/1‑139 (MQ=255)
                                               gCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGa                                                                                                                                >  2:140708/1‑94 (MQ=255)
                                               gCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGa                                                                                                                                <  1:140708/94‑1 (MQ=255)
                                               gCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGt                                                                                   <  2:178867/139‑1 (MQ=255)
                                               gCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGt                                                                                   <  2:838415/139‑1 (MQ=255)
                                               gCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGt                                                                                   <  2:495020/139‑1 (MQ=255)
                                                cAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAAc                                                                                                         <  2:178137/116‑1 (MQ=255)
                                                cAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAAc                                                                                                         >  1:178137/1‑116 (MQ=255)
                                                   aaCTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGa                                                                               >  1:1337788/1‑139 (MQ=255)
                                                    aCTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGaa                                                                              <  1:1429391/139‑1 (MQ=255)
                                                       gCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCAt                                                                                                >  1:775916/1‑118 (MQ=255)
                                                       gCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCAt                                                                                                <  2:775916/118‑1 (MQ=255)
                                                           aTTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTtgct                                                                       >  2:290228/1‑139 (MQ=255)
                                                              aaaTGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTgcg                                                                    <  2:93043/139‑1 (MQ=255)
                                                               aaTGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAc                                                                                                                               >  1:552946/1‑79 (MQ=255)
                                                               aaTGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGAc                                                                                                                               <  2:552946/79‑1 (MQ=255)
                                                                 tGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGTTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAg                                                                 <  1:5067/139‑1 (MQ=255)
                                                                 tGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAg                                                                 <  1:87780/139‑1 (MQ=255)
                                                                 tGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAg                                                                 <  1:876958/139‑1 (MQ=255)
                                                                 tGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAg                                                                 <  1:744499/139‑1 (MQ=255)
                                                                      ggCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTc                                                                                         <  1:1491279/110‑1 (MQ=255)
                                                                      ggCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTc                                                                                         >  2:1491279/1‑110 (MQ=255)
                                                                                 gAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGc                                                 <  1:1431752/139‑1 (MQ=255)
                                                                                     ttGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAAt                                                                                                                                   >  2:1010523/1‑53 (MQ=255)
                                                                                     ttGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAAt                                                                                                                                   <  1:1010523/53‑1 (MQ=255)
                                                                                     ttGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTAATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCt                                             >  1:955949/1‑139 (MQ=255)
                                                                                         ttttGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGaa                                                                              >  1:741469/1‑102 (MQ=255)
                                                                                         ttttGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGaa                                                                              <  2:741469/102‑1 (MQ=255)
                                                                                            tGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTc                                                                                                                    <  2:1026442/61‑1 (MQ=255)
                                                                                            tGGTGATATCGGCGCATCGTTATTTAGTGAGGCAGGGCCGGGCAATCGACCTGCTGGTGTc                                                                                                                    >  1:1026442/1‑61 (MQ=255)
                                                                                              gTGATATCGGCGCATCGTTTTTTAGTGCGGCAGGGCCGGGCAATCGACctgc                                                                                                                           <  1:496492/52‑1 (MQ=255)
                                                                                              gTGATATCGGCGCATCGTTTTTTAGTGCGGCAGGGCCGGGCAATCGACctgc                                                                                                                           >  2:496492/1‑52 (MQ=255)
                                                                                                  tatCGGCGCATCGTTATTTAGTGCGACAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGcc                                >  1:808126/1‑139 (MQ=255)
                                                                                                       gcgcATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTg                                            >  1:1023329/1‑122 (MQ=255)
                                                                                                       gcgcATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTg                                            <  2:1023329/122‑1 (MQ=255)
                                                                                                        cgcATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTaa                                                       >  2:1264939/1‑110 (MQ=255)
                                                                                                        cgcATCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTaa                                                       <  1:1264939/110‑1 (MQ=255)
                                                                                                           aTCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTtg                                                                         <  1:1352048/89‑1 (MQ=255)
                                                                                                           aTCGTTATTTAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTtg                                                                         >  2:1352048/1‑89 (MQ=255)
                                                                                                                  tttAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTc                                                                                         >  1:439862/1‑66 (MQ=255)
                                                                                                                  tttAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTc                                                                                         <  2:439862/66‑1 (MQ=255)
                                                                                                                    tAGTGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATg              <  1:94664/139‑1 (MQ=255)
                                                                                                                       tGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGc                                                                                                             <  1:842681/41‑1 (MQ=255)
                                                                                                                       tGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGc                                                                                                             >  2:842681/1‑41 (MQ=255)
                                                                                                                       tGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTc                                                                                         >  2:584391/1‑61 (MQ=255)
                                                                                                                       tGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTc                                                                                         <  1:584391/61‑1 (MQ=255)
                                                                                                                       tGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTgg                        >  2:59829/1‑126 (MQ=255)
                                                                                                                       tGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTgg                        <  1:59829/126‑1 (MQ=255)
                                                                                                                       tGCGGCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGgtg           >  2:82640/1‑139 (MQ=255)
                                                                                                                          ggCAGGGCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCg        <  2:655742/139‑1 (MQ=255)
                                                                                                                                gCCGGGCAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAAt  <  1:1473350/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
GTTACACCTGTTTATACCATTAGTTTTGGTGGCCGGGTTGAAGTACCGCAAAACTGCGAATTAAATGCCGGGCAGATTGTTGAATTTGATTTTGGTGATATCGGCGCATCGTTATTTAGTGCGGCAGGGCCGGGTAATCGACCTGCTGGTGTCATGCCGCAAACCAAGAGCATTGCGGTCAAATGTACGAATGTTGCTGCGCAGGCTTATTTAACAATGCGTCTGGAAGCCAGTGCCGTTTCTGGTCAGGCGATGGTGTCGGACAAT  >  NC_000913/562852‑563118

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