Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,157,712 (C)6→7 coding (63/927 nt) yqhG → DUF3828 domain‑containing protein YqhG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,157,7061.C100.0% 1007.4 / NA 238Q19H (CAG→CAC) yqhGDUF3828 domain‑containing protein YqhG
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (136/102);  total (136/102)

Only 100 of 268 total aligned reads displayed.

CTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑CCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCG  >  NC_000913/3157606‑3157801
                                                                                                     |                                                                                               
 tCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑c                                                                                                <  1:4616010/101‑1 (MQ=255)
 tCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑c                                                                                                <  1:4615706/101‑1 (MQ=255)
   aTCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑ccc                                                                                              >  2:4616135/1‑101 (MQ=255)
     cTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑ccccc                                                                                            >  2:4616139/1‑101 (MQ=255)
     cTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑c                                                                                                >  1:4616140/1‑97 (MQ=255)
       atatCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑ccccccc                                                                                          >  1:4616145/1‑100 (MQ=255)
        tatCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑ccccccct                                                                                         >  2:4616146/1‑99 (MQ=255)
         atCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑ccccccctc                                                                                        <  1:4615711/101‑4 (MQ=255)
           ccGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑ccccccctcac                                                                                      <  1:4615681/101‑6 (MQ=255)
           ccGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑ccccccctcac                                                                                      <  2:4615563/101‑6 (MQ=255)
              ggATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGa                                                                                    >  2:4616150/1‑101 (MQ=255)
                 tAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCg                                                                                 >  1:4616152/1‑101 (MQ=255)
                  aaaCGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGt                                                                                >  2:4616157/1‑101 (MQ=255)
                  aaaCGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGt                                                                                >  1:4616156/1‑101 (MQ=255)
                   aaCGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTa                                                                               <  1:4615902/101‑1 (MQ=255)
                     cGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGa                                                                             >  2:4616161/1‑101 (MQ=255)
                     cGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGa                                                                             <  2:4615873/101‑1 (MQ=255)
                      gCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGaa                                                                            <  1:4615601/101‑1 (MQ=255)
                      gCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGaa                                                                            >  1:4616163/1‑101 (MQ=255)
                        aaCGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGaaca                                                                          >  1:4616165/1‑101 (MQ=255)
                        aaCGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGaaca                                                                          >  2:4616169/1‑101 (MQ=255)
                         aCGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGaacaa                                                                         >  1:4616171/1‑101 (MQ=255)
                          cGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGaacaaa                                                                        <  1:4615878/101‑1 (MQ=255)
                            gagaGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGaacaaaca                                                                      <  1:4615835/101‑1 (MQ=255)
                             agagGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAg                                                                     >  2:4616173/1‑101 (MQ=255)
                               agGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTc                                                                   <  1:4615858/101‑1 (MQ=255)
                                ggTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTcc                                                                  >  2:4616176/1‑101 (MQ=255)
                                 gTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCg                                                                 <  1:4615761/101‑1 (MQ=255)
                                   gAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGc                                                               <  1:4615941/101‑1 (MQ=255)
                                   gAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGc                                                               <  1:4615933/101‑1 (MQ=255)
                                   gAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGc                                                               >  1:4616179/1‑101 (MQ=255)
                                     aaaaCAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCat                                                             <  2:4615764/101‑1 (MQ=255)
                                     aaaaCAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCat                                                             >  2:4616182/1‑101 (MQ=255)
                                       aaCAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCatat                                                           >  1:4616185/1‑101 (MQ=255)
                                       aaCAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCatat                                                           >  2:4616187/1‑101 (MQ=255)
                                        aCAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATAtt                                                          >  1:4616189/1‑101 (MQ=255)
                                         cAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATAttt                                                         <  2:4615812/101‑1 (MQ=255)
                                          aGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTa                                                        >  1:4616195/1‑101 (MQ=255)
                                          aGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTa                                                        >  2:4616194/1‑101 (MQ=255)
                                          aGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTa                                                        <  2:4615778/101‑1 (MQ=255)
                                            aTGAAAATAATTCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATc                                                      <  2:4615895/101‑1 (MQ=255)
                                            aTGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATc                                                      >  1:4616197/1‑101 (MQ=255)
                                             tGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCa                                                     >  2:4616199/1‑101 (MQ=255)
                                              gAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAg                                                    <  1:4615800/101‑1 (MQ=255)
                                               aaaaTAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGa                                                   <  1:4615409/101‑1 (MQ=255)
                                               aaaaTAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGa                                                   >  2:4616203/1‑101 (MQ=255)
                                                aaaTAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGaa                                                  >  1:4616206/1‑101 (MQ=255)
                                                 aataatACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAAc                                                 <  2:4615801/101‑1 (MQ=255)
                                                  ataataCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACt                                                >  2:4616209/1‑101 (MQ=255)
                                                   taataCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACta                                               >  2:4616211/1‑101 (MQ=255)
                                                    aataCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACtat                                              <  1:4616018/101‑1 (MQ=255)
                                                     ataCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACtata                                             >  2:4616214/1‑101 (MQ=255)
                                                     ataCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACtata                                             >  1:4616216/1‑101 (MQ=255)
                                                      taCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATaa                                            <  2:4615671/101‑1 (MQ=255)
                                                       aCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATaaa                                           <  1:4615949/101‑1 (MQ=255)
                                                        cTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAAt                                          >  1:4616219/1‑101 (MQ=255)
                                                         ttCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACAGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATc                                         >  2:4616222/1‑101 (MQ=255)
                                                          tCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCa                                        >  2:4616224/1‑101 (MQ=255)
                                                           cTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAg                                       <  2:4615644/101‑1 (MQ=255)
                                                             gTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTcc                                                                  >  2:4616229/1‑72 (MQ=255)
                                                             gTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGAt                                     <  2:4616045/101‑1 (MQ=255)
                                                             gTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGAt                                     >  1:4616228/1‑101 (MQ=255)
                                                              tttttAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATg                                    <  2:4615926/101‑1 (MQ=255)
                                                               ttttAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGc                                   >  1:4616232/1‑101 (MQ=255)
                                                               ttttAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGc                                   <  2:4615750/101‑1 (MQ=255)
                                                                tttAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGcc                                  >  1:4616234/1‑101 (MQ=255)
                                                                tttAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGcc                                  <  1:4615746/101‑1 (MQ=255)
                                                                  tAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCAc                                <  2:4615942/101‑1 (MQ=255)
                                                                  tAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCAc                                >  2:4616239/1‑101 (MQ=255)
                                                                   agcagcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACt                               >  1:4616240/1‑101 (MQ=255)
                                                                    gcagcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTg                              <  2:4615973/101‑1 (MQ=255)
                                                                    gcagcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTg                              <  1:4615796/101‑1 (MQ=255)
                                                                    gcagcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTg                              <  1:4615837/101‑1 (MQ=255)
                                                                     cagcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGc                             >  2:4616244/1‑101 (MQ=255)
                                                                     cagcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGc                             <  2:4615930/101‑1 (MQ=255)
                                                                     cagcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGc                             >  1:4616245/1‑101 (MQ=255)
                                                                      agcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAAt                                          >  2:4616249/1‑87 (MQ=255)
                                                                       gcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGccc                           >  1:4616250/1‑101 (MQ=255)
                                                                         ccTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGCAg                                                                              >  2:4616255/1‑48 (MQ=255)
                                                                          cTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCtt                        <  2:4615842/101‑1 (MQ=255)
                                                                          cTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCtt                        >  2:4616257/1‑101 (MQ=255)
                                                                           tGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTa                       >  2:4616260/1‑101 (MQ=255)
                                                                             gCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTAtt                     <  1:4615911/101‑1 (MQ=255)
                                                                               aaGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTAtttt                   >  1:4616265/1‑101 (MQ=255)
                                                                               aaGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTAtttt                   >  2:4616263/1‑101 (MQ=255)
                                                                               aaGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTAtttt                   >  2:4616266/1‑101 (MQ=255)
                                                                                aGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTg                  <  1:4615818/101‑1 (MQ=255)
                                                                                 gTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGGCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTgg                 <  1:4615965/101‑1 (MQ=255)
                                                                                  ttttACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGt                <  2:4615867/101‑1 (MQ=255)
                                                                                   tttACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTg               <  1:4615822/101‑1 (MQ=255)
                                                                                    ttACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGa              <  1:4616000/101‑1 (MQ=255)
                                                                                     tACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGaa             >  2:4616271/1‑101 (MQ=255)
                                                                                         gTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCgg                                                                >  1:4616272/1‑46 (MQ=255)
                                                                                             acGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGgaga     >  1:4616277/1‑101 (MQ=255)
                                                                                             acGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGgaga     <  1:4615782/101‑1 (MQ=255)
                                                                                              cGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGagag    >  1:4616279/1‑101 (MQ=255)
                                                                                               gCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGc   <  1:4615890/101‑1 (MQ=255)
                                                                                               gCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGc   <  1:4615865/101‑1 (MQ=255)
                                                                                                cacaGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCg  >  1:4616285/1‑101 (MQ=255)
                                                                                                cacaGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCg  >  2:4616282/1‑101 (MQ=255)
                                                                                                     |                                                                                               
CTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑CCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCG  >  NC_000913/3157606‑3157801

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