Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,509,873 A→G pseudogene (70/330 nt) yjhV → KpLE2 phage‑like element; putative uncharacterized protein YjhV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,509,8730AG100.0% 1067.9 / NA 258pseudogene (70/330 nt)yjhVKpLE2 phage‑like element; putative uncharacterized protein YjhV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (130/128);  total (130/128)

Only 100 of 261 total aligned reads displayed.

GGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCACGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCACACCATCCTATTGTTG  >  NC_000913/4509775‑4509970
                                                                                                  |                                                                                                 
ggAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTcgcg                                                                                                 <  2:5969444/101‑1 (MQ=255)
 gAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGa                                                                                                <  2:5969452/101‑1 (MQ=255)
  aGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGCCGCGAt                                                                                               >  2:5969968/1‑101 (MQ=255)
    gCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGa                                                                                             >  1:5969972/1‑101 (MQ=255)
     ccAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGAc                                                                                            >  2:5969975/1‑101 (MQ=255)
     ccAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGAc                                                                                            >  2:5969973/1‑101 (MQ=255)
      cAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGAct                                                                                           >  1:5969979/1‑101 (MQ=255)
       aGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGActc                                                                                          >  2:5969982/1‑101 (MQ=255)
       aGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGActc                                                                                          >  1:5969983/1‑101 (MQ=255)
        gTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGActct                                                                                         <  1:5969542/101‑1 (MQ=255)
        gTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGActct                                                                                         <  1:5969588/101‑1 (MQ=255)
         ttttAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTg                                                                                        >  2:5969985/1‑101 (MQ=255)
          tttAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTgg                                                                                       <  2:5969537/101‑1 (MQ=255)
          tttAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTgg                                                                                       <  1:5969803/101‑1 (MQ=255)
           ttAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGt                                                                                      >  2:5969989/1‑101 (MQ=255)
            tAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGtt                                                                                     <  2:5969497/101‑1 (MQ=255)
             aGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTa                                                                                    >  2:5969992/1‑101 (MQ=255)
             aGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTa                                                                                    <  1:5969662/101‑1 (MQ=255)
                  aCCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGAt                                                                               <  2:5969504/101‑1 (MQ=255)
                   ccACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATa                                                                              <  2:5969525/101‑1 (MQ=255)
                    cacaACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATAc                                                                             >  2:5969997/1‑101 (MQ=255)
                    cacaACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATAc                                                                             >  1:5969996/1‑101 (MQ=255)
                      caacaGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCg                                                                           <  1:5969563/101‑1 (MQ=255)
                       aacaGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCgg                                                                          <  2:5969512/101‑1 (MQ=255)
                        acaGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGc                                                                         <  1:5969486/101‑1 (MQ=255)
                         caGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCt                                                                        <  2:5969560/101‑1 (MQ=255)
                          aGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCtt                                                                       >  2:5970003/1‑101 (MQ=255)
                          aGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCtt                                                                       >  1:5970005/1‑101 (MQ=255)
                          aGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCtt                                                                       <  2:5969770/101‑1 (MQ=255)
                           gACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCttt                                                                      <  1:5969932/101‑1 (MQ=255)
                            aCTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTg                                                                     >  1:5970007/1‑101 (MQ=255)
                             cTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGa                                                                    >  2:5970010/1‑101 (MQ=255)
                              tACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGAc                                                                   >  1:5970012/1‑101 (MQ=255)
                              tACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGAc                                                                   >  1:5970014/1‑101 (MQ=255)
                                ccTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGtt                                                                                     <  2:5969301/81‑1 (MQ=255)
                                ccTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGAccc                                                                 >  2:5970016/1‑101 (MQ=255)
                                   gAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGt                                                              >  1:5970021/1‑101 (MQ=255)
                                     aGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTAt                                                            <  2:5969620/101‑1 (MQ=255)
                                      gggAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTAtt                                                           <  1:5969198/101‑1 (MQ=255)
                                       ggAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTAttt                                                          <  2:5969582/101‑1 (MQ=255)
                                       ggAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTAttt                                                          >  1:5970024/1‑101 (MQ=255)
                                         aGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTc                                                        <  1:5969426/101‑1 (MQ=255)
                                          gCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTcc                                                       >  2:5970026/1‑101 (MQ=255)
                                          gCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTcc                                                       <  2:5969488/101‑1 (MQ=255)
                                           ccGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCg                                                      >  2:5970028/1‑101 (MQ=255)
                                             gcgGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGc                                                    >  2:5970032/1‑101 (MQ=255)
                                             gcgGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGc                                                    <  1:5969636/101‑1 (MQ=255)
                                             gcgGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGc                                                    >  1:5970030/1‑101 (MQ=255)
                                              cgGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGct                                                   <  2:5969575/101‑1 (MQ=255)
                                               ggTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGctc                                                  <  2:5969631/101‑1 (MQ=255)
                                                gTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGctct                                                 <  2:5969664/101‑1 (MQ=255)
                                                 tCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTa                                                <  2:5969698/101‑1 (MQ=255)
                                                   gCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATc                                              <  2:5969370/101‑1 (MQ=255)
                                                   gCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATc                                              >  1:5970036/1‑101 (MQ=255)
                                                    ccTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCg                                             >  2:5970039/1‑101 (MQ=255)
                                                     cTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGc                                            >  2:5970041/1‑101 (MQ=255)
                                                       ggCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTg                                          <  2:5969607/101‑1 (MQ=255)
                                                       ggCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTg                                          >  1:5970042/1‑101 (MQ=255)
                                                        gCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGa                                         >  2:5970045/1‑101 (MQ=255)
                                                         cAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATGCCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTg                                          <  1:5969514/99‑1 (MQ=255)
                                                         cAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGaa                                        <  1:5969568/101‑1 (MQ=255)
                                                            ttGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCa                                     >  1:5970049/1‑101 (MQ=255)
                                                            ttGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCa                                     >  2:5970050/1‑101 (MQ=255)
                                                            ttGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCa                                     <  1:5969551/101‑1 (MQ=255)
                                                               cAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATgg                                  >  1:5970055/1‑101 (MQ=255)
                                                                aGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGc                                 <  1:5969725/101‑1 (MQ=255)
                                                                 gTAGCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCg                                >  2:5970057/1‑101 (MQ=255)
                                                                    gCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTca                             <  1:5969547/101‑1 (MQ=255)
                                                                    gCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTca                             <  1:5969611/101‑1 (MQ=255)
                                                                    gCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTca                             <  2:5969775/101‑1 (MQ=255)
                                                                    gCAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTca                             >  1:5970058/1‑101 (MQ=255)
                                                                     cAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCAt                                    >  2:5970062/1‑93 (MQ=255)
                                                                     cAGGAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTcag                            <  1:5969510/101‑1 (MQ=255)
                                                                       ggAGATCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTcagca                          <  2:5969673/101‑1 (MQ=255)
                                                                        gagaTCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTcagcag                         >  1:5970064/1‑101 (MQ=255)
                                                                         agaTCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCagcagc                        <  1:5969349/101‑1 (MQ=255)
                                                                         agaTCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCagcagc                        <  1:5969764/101‑1 (MQ=255)
                                                                          gaTCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGcc                       >  2:5970066/1‑101 (MQ=255)
                                                                          gaTCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGcc                       <  1:5969459/101‑1 (MQ=255)
                                                                          gaTCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGcc                       <  1:5969697/101‑1 (MQ=255)
                                                                          gaTCATCCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGcc                       <  2:5969587/101‑1 (MQ=255)
                                                                            tcatcCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGt                     <  2:5969736/101‑1 (MQ=255)
                                                                             catcCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTg                    >  2:5970070/1‑101 (MQ=255)
                                                                              atcCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGatc                                                                  >  1:5970073/1‑52 (MQ=39)
                                                                               tcCAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGca                  <  2:5969779/101‑1 (MQ=255)
                                                                                ccAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGcac                 >  1:5970076/1‑101 (MQ=255)
                                                                                ccAGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGcac                 >  2:5970075/1‑101 (MQ=255)
                                                                                  aGCCGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCacac               >  1:5970079/1‑101 (MQ=255)
                                                                                     cGCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTTAATCATGGCGTCAGCAGCCGTGCACACCAt            <  2:5969605/101‑1 (MQ=255)
                                                                                      gCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCACACCATc           >  1:5970084/1‑101 (MQ=255)
                                                                                      gCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCACACCATc           <  2:5969628/101‑1 (MQ=255)
                                                                                      gCAGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCACACCATc           >  2:5970083/1‑101 (MQ=255)
                                                                                        aGTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCACACCATCCt         >  2:5970088/1‑101 (MQ=255)
                                                                                         gTACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCACACCATCCTa        <  1:5969676/101‑1 (MQ=255)
                                                                                          tACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCACACCATCCTAt       >  1:5970093/1‑101 (MQ=255)
                                                                                          tACAGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCACACCATCCTAt       >  2:5970091/1‑101 (MQ=255)
                                                                                             aGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCACACCATCCTAttgt    >  2:5970096/1‑101 (MQ=255)
                                                                                             aGGTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCACACCATCCTAttgt    <  1:5969592/101‑1 (MQ=255)
                                                                                              ggTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCACACCATCCTAttgtt   <  2:5969584/101‑1 (MQ=255)
                                                                                               gTCGCGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCACACCATCCTAttgttg  <  1:5969642/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
GGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCACGATGACTCTGGTTAATGATACCGGCTTTGACCCCGTATTTTCCGGCTCTATCGCTGAATCATGGCGTCAGCAGCCGTGCACACCATCCTATTGTTG  >  NC_000913/4509775‑4509970

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