Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 478,533 C→T A29T (GCC→ACC)  ylaB ← putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913478,5330CT100.0% 253.3 / NA 77A29T (GCC→ACC) ylaBputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/37);  total (40/37)

AGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGCCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGT  >  NC_000913/478402‑478663
                                                                                                                                   |                                                                                                                                  
aGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGtcagccag                                                                                                                             >  2:973007/1‑139 (MQ=255)
  cgcCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGAt                                                                                                                           >  2:32765/1‑139 (MQ=255)
       tGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTa                                                                                                                      >  2:403482/1‑139 (MQ=255)
        gCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTaa                                                                                                                     >  2:225434/1‑139 (MQ=255)
               ttCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACa                                                                                                              <  2:553058/139‑1 (MQ=255)
               ttCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACa                                                                                                              <  1:934231/139‑1 (MQ=255)
                    gCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                         <  2:702634/139‑1 (MQ=255)
                        cTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAAt                                                                                                     >  1:236824/1‑139 (MQ=255)
                        cTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAAt                                                                                                     >  1:374392/1‑139 (MQ=255)
                          gCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATAc                                                                                                   >  1:216976/1‑139 (MQ=255)
                             gcgcATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGa                                                                                                >  1:123513/1‑139 (MQ=255)
                                cATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAg                                                                                             >  1:183703/1‑139 (MQ=255)
                                             tgtCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGa                                                                                                                >  1:154976/1‑107 (MQ=255)
                                             tgtCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGa                                                                                                                <  2:154976/107‑1 (MQ=255)
                                             tgtCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACa                                                                                                              <  2:242633/109‑1 (MQ=255)
                                             tgtCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACa                                                                                                              >  1:242633/1‑109 (MQ=255)
                                             tgtCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                                     >  1:715392/1‑134 (MQ=255)
                                             tgtCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTccc                                                                                >  1:1025/1‑139 (MQ=255)
                                             tgtCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACAATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                                     <  2:715392/134‑1 (MQ=255)
                                               tCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              <  2:677895/91‑1 (MQ=255)
                                               tCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              >  1:677895/1‑91 (MQ=255)
                                               tCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGa                                                                              >  1:1128836/1‑139 (MQ=255)
                                                 gCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGaaa                                                                            <  1:225434/139‑1 (MQ=255)
                                                  ccACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAAt                                                                           <  2:236824/139‑1 (MQ=255)
                                                               tCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTa                                                                                                                      >  1:1091586/1‑83 (MQ=255)
                                                               tCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTa                                                                                                                      <  2:1091586/83‑1 (MQ=255)
                                                                 gAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATg                                                            >  2:954961/1‑139 (MQ=255)
                                                                  aaTAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCAAAATAAGGCCGACCAGATGt                                                           >  2:309684/1‑139 (MQ=255)
                                                                       cGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTc                                                          >  1:597119/1‑135 (MQ=255)
                                                                       cGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTc                                                          <  2:597119/135‑1 (MQ=255)
                                                                         aCGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                         <  2:529573/86‑1 (MQ=255)
                                                                         aCGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                         >  1:529573/1‑86 (MQ=255)
                                                                             gggAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCAc                                                <  2:1025/139‑1 (MQ=255)
                                                                                             cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCaa                                                                                                      <  2:86482/69‑1 (MQ=255)
                                                                                             cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCaa                                                                                                      >  1:86482/1‑69 (MQ=255)
                                                                                             cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCa                                                 >  1:160090/1‑122 (MQ=255)
                                                                                             cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCa                                                 <  2:160090/122‑1 (MQ=255)
                                                                                             cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCttatta                                >  1:1094793/1‑139 (MQ=255)
                                                                                             cAGCTCTTCAATAAACGATGTTTCTACCCGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAGTACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCttatta                                >  1:704052/1‑139 (MQ=255)
                                                                                               gCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTcc                                                                                 <  1:79018/88‑1 (MQ=255)
                                                                                               gCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTcc                                                                                 >  2:79018/1‑88 (MQ=255)
                                                                                                ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             <  2:216976/139‑1 (MQ=255)
                                                                                                ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             <  1:32765/139‑1 (MQ=255)
                                                                                                ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             <  2:123513/139‑1 (MQ=255)
                                                                                                ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             >  2:249884/1‑139 (MQ=255)
                                                                                                ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             <  2:1094793/139‑1 (MQ=255)
                                                                                                ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             <  1:249884/139‑1 (MQ=255)
                                                                                                  ctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCata                                         <  1:283660/125‑1 (MQ=255)
                                                                                                  ctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCata                                         >  2:283660/1‑125 (MQ=255)
                                                                                                     cAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAgtg                        <  2:704052/139‑1 (MQ=255)
                                                                                                         aaaCGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAgg                    <  1:973007/139‑1 (MQ=255)
                                                                                                             gATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGaa                                                                                               <  1:893699/60‑1 (MQ=255)
                                                                                                             gATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGaa                                                                                               >  2:893699/1‑60 (MQ=255)
                                                                                                             gATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGc                                            <  2:39232/111‑1 (MQ=255)
                                                                                                             gATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGc                                            >  1:39232/1‑111 (MQ=255)
                                                                                                               tGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                         >  2:803775/1‑48 (MQ=255)
                                                                                                               tGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                         <  1:803775/48‑1 (MQ=255)
                                                                                                                gTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGta             <  2:714360/139‑1 (MQ=255)
                                                                                                                     tACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCt                                     >  2:1065152/1‑110 (MQ=255)
                                                                                                                     tACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCt                                     <  1:1065152/110‑1 (MQ=255)
                                                                                                                         tgctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCGCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATc                                       >  1:597422/1‑104 (MQ=255)
                                                                                                                         tgctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATc                                       <  2:597422/104‑1 (MQ=255)
                                                                                                                         tgctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGAtt    >  1:660226/1‑139 (MQ=255)
                                                                                                                          gctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGAcc                                                                 <  1:1111038/77‑1 (MQ=255)
                                                                                                                          gctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGAcc                                                                 >  2:1111038/1‑77 (MQ=255)
                                                                                                                           ctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                                     <  2:745001/56‑1 (MQ=255)
                                                                                                                           ctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                                     >  1:745001/1‑56 (MQ=255)
                                                                                                                           ctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATc                                       >  1:763051/1‑102 (MQ=255)
                                                                                                                           ctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATc                                       <  2:763051/102‑1 (MQ=255)
                                                                                                                           ctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAac                              >  1:846871/1‑111 (MQ=255)
                                                                                                                           ctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAac                              <  2:846871/111‑1 (MQ=255)
                                                                                                                           ctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATtgt  >  2:1083131/1‑139 (MQ=255)
                                                                                                                              aTGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGaa                                                                                               <  1:961462/43‑1 (MQ=255)
                                                                                                                              aTGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGaa                                                                                               >  2:961462/1‑43 (MQ=255)
                                                                                                                              aTGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGc                                                                      <  2:1010432/68‑1 (MQ=255)
                                                                                                                              aTGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGc                                                                      >  1:1010432/1‑68 (MQ=255)
                                                                                                                              aTGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             <  2:705672/109‑1 (MQ=255)
                                                                                                                              aTGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             >  1:705672/1‑109 (MQ=255)
                                                                                                                               tGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGaa                                                                                               <  1:864439/42‑1 (MQ=255)
                                                                                                                               tGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGaa                                                                                               >  2:864439/1‑42 (MQ=255)
                                                                                                                                   |                                                                                                                                  
AGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGCCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGT  >  NC_000913/478402‑478663

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: ‑