Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,681,932 T→C P33P (CCA→CCG ydgC ← GlpM family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,681,9320TC100.0% 319.7 / NA 87P33P (CCA→CCGydgCGlpM family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (41/46);  total (41/46)

CGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCCC  >  NC_000913/1681836‑1682028
                                                                                                |                                                                                                
cGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGa                                                                                               >  2:373907/1‑100 (MQ=255)
cGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGa                                                                                               <  2:1650320/100‑1 (MQ=255)
 gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                              <  1:1919550/100‑1 (MQ=255)
 gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGAAAAGGTCGGGAAAAGCGGaa                                                                                              <  2:341881/100‑1 (MQ=255)
  aCCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAAt                                                                                             >  1:648137/1‑100 (MQ=255)
  aCCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAAt                                                                                             >  2:267976/1‑100 (MQ=255)
    cacaTACTAAAAATGATTGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCa                                                                                           <  2:1486945/100‑1 (MQ=255)
    cacaTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCa                                                                                           <  2:1075941/100‑1 (MQ=255)
     acaTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAg                                                                                          >  1:2149653/1‑100 (MQ=255)
     acaTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAg                                                                                          <  2:877649/100‑1 (MQ=255)
     acaTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAg                                                                                          >  1:1476438/1‑100 (MQ=255)
     acaTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCCGGAACAGCGGAATCAg                                                                                          >  2:1395855/1‑100 (MQ=255)
      caTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGc                                                                                         >  1:2192590/1‑100 (MQ=255)
       aTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGcc                                                                                        <  2:1723218/100‑1 (MQ=255)
        tACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGccc                                                                                       >  2:655308/1‑100 (MQ=255)
         aCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCg                                                                                      >  2:459125/1‑100 (MQ=255)
          cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                     <  1:2279831/100‑1 (MQ=255)
            aaaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCg                                                                                   <  2:1187243/100‑1 (MQ=255)
            aaaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCg                                                                                   >  1:1045090/1‑100 (MQ=255)
             aaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGa                                                                                  >  1:1115512/1‑100 (MQ=255)
             aaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGa                                                                                  >  2:1280555/1‑100 (MQ=255)
              aaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGat                                                                                 <  1:1026261/100‑1 (MQ=255)
               aaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGata                                                                                <  2:81523/100‑1 (MQ=255)
                atgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGatat                                                                               <  2:1374387/100‑1 (MQ=255)
                      gTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATaataat                                                                         >  1:1081150/1‑100 (MQ=255)
                        tGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttt                                                                       >  2:741045/1‑100 (MQ=255)
                         gcgcGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAtttt                                                                      <  1:594431/100‑1 (MQ=255)
                          cgcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttttt                                                                     >  1:1065592/1‑100 (MQ=255)
                          cgcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttttt                                                                     >  2:1969984/1‑100 (MQ=255)
                           gcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAATGGCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTc                                                                    >  2:777110/1‑100 (MQ=255)
                           gcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTc                                                                    <  1:1280558/100‑1 (MQ=255)
                            cgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCg                                                                   >  1:74125/1‑100 (MQ=255)
                                aGGCTTTAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGtttt                                                               <  1:2333981/100‑1 (MQ=255)
                                 ggCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttttt                                                              <  2:1166606/100‑1 (MQ=255)
                                  gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                             <  2:648136/100‑1 (MQ=255)
                                  gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                             >  1:582787/1‑100 (MQ=255)
                                  gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                             <  1:2002357/100‑1 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                            <  1:486854/100‑1 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                            <  2:92661/100‑1 (MQ=255)
                                    ttCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCt                                                           <  2:1777316/100‑1 (MQ=255)
                                     tCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTa                                                          >  2:2295562/1‑100 (MQ=255)
                                      cAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaa                                                         <  2:696139/100‑1 (MQ=255)
                                        aTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCGAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaaa                                                       >  2:1281378/1‑100 (MQ=255)
                                         tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaa                                                        >  1:1848413/1‑98 (MQ=255)
                                          gCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAaca                                                     >  2:1916778/1‑100 (MQ=255)
                                          gCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAaca                                                     >  1:613872/1‑100 (MQ=255)
                                            cgcgTTCTCTGGCAACAATATAATGCGCAATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACAcc                                                   <  1:1643165/100‑1 (MQ=255)
                                            cgcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGTCCGGCGATATAATAATTTTTCGTTCTTGCTAAAACAcc                                                   <  2:1442030/100‑1 (MQ=255)
                                            cgcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACAcc                                                   <  1:1930578/100‑1 (MQ=255)
                                            cgcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACAcc                                                   >  1:277140/1‑100 (MQ=255)
                                              cgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCaa                                                 <  1:373909/100‑1 (MQ=255)
                                                ttCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATc                                               >  1:1044554/1‑100 (MQ=255)
                                                ttCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATc                                               <  1:665910/100‑1 (MQ=255)
                                                gatcttGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATc                                               <  2:391874/95‑1 (MQ=255)
                                                    cTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                           <  2:1044552/100‑1 (MQ=255)
                                                    cTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                           <  2:420174/100‑1 (MQ=255)
                                                     tGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCagatc                                          >  1:391876/1‑96 (MQ=255)
                                                     tGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAg                                          <  2:881156/100‑1 (MQ=255)
                                                       gCAACAATATAATGCGCGATAAGCGCAAAGGTGGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCa                                        >  1:2214785/1‑100 (MQ=255)
                                                       gCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCa                                        >  1:135118/1‑100 (MQ=255)
                                                       gCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCa                                        >  1:45240/1‑100 (MQ=255)
                                                        caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc                                       >  2:865169/1‑100 (MQ=255)
                                                        caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc                                       >  2:1996857/1‑100 (MQ=255)
                                                           caaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTAc                                    >  2:1180801/1‑100 (MQ=255)
                                                              tataATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAg                                 <  2:135116/100‑1 (MQ=255)
                                                               ataATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgc                                <  2:1211325/100‑1 (MQ=255)
                                                                   tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                            <  1:622900/100‑1 (MQ=255)
                                                                      gcgATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGa                         >  2:638154/1‑100 (MQ=255)
                                                                       cgATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAg                        >  1:997264/1‑100 (MQ=255)
                                                                          tAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAgcgc                     >  1:1276229/1‑100 (MQ=255)
                                                                          tAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAgcgc                     >  2:814854/1‑100 (MQ=255)
                                                                           aaGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTACAACACCAATCAACAGCACTACCAGCGCACCGAGcgcg                    <  2:1249873/100‑1 (MQ=255)
                                                                            aGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCgg                   >  1:222132/1‑100 (MQ=255)
                                                                            aGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCgg                   >  1:1630086/1‑100 (MQ=255)
                                                                             gcgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGc                  >  1:1974193/1‑100 (MQ=255)
                                                                               gcAAATGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCtt                <  2:624469/100‑1 (MQ=255)
                                                                                cAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCttt               <  1:814856/100‑1 (MQ=255)
                                                                                cAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCttt               >  1:2311410/1‑100 (MQ=255)
                                                                                 aaaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTg              <  2:277138/100‑1 (MQ=255)
                                                                                   aGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAt            <  2:1044336/100‑1 (MQ=255)
                                                                                   aGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAt            <  1:1281381/100‑1 (MQ=255)
                                                                                    ggTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAtt           <  1:2295571/100‑1 (MQ=255)
                                                                                    ggTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAtt           <  2:1045088/100‑1 (MQ=255)
                                                                                       cGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAcc        <  1:638155/100‑1 (MQ=255)
                                                                                       cGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAcc        <  1:625916/100‑1 (MQ=255)
                                                                                            aaaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGccc   <  1:759684/100‑1 (MQ=255)
                                                                                             aaGCGGAATCAGCCCGGCGATATAATAATTGTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGcccc  <  2:222130/100‑1 (MQ=255)
                                                                                                |                                                                                                
CGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCCC  >  NC_000913/1681836‑1682028

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: