Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,500,711 T→C pseudogene (81/1199 nt) yjgX ← KpLE2 phage‑like element; putative phosphoethanolamine transferase YjgX, N‑terminal fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,500,7110TC99.6% 1104.6 / ‑6.4 264pseudogene (81/1199 nt)yjgXKpLE2 phage‑like element; putative phosphoethanolamine transferase YjgX, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (122/141);  minor base G (0/1);  total (122/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 266 total aligned reads displayed.

CCCACAACTGTTCCTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGTGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAG  >  NC_000913/4500613‑4500808
                                                                                                  |                                                                                                 
    caACTGTTCCTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGc                                                                                             >  1:6672267/1‑101 (MQ=255)
    caACTGTTCCTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGc                                                                                             >  2:6672266/1‑101 (MQ=255)
     aaCTGTTCCTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGcc                                                                                            <  2:6671994/101‑1 (MQ=255)
     aaCTGTTCCTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGcc                                                                                            <  1:6671851/101‑1 (MQ=255)
      aCTGTTCCTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCa                                                                                           <  1:6671663/101‑1 (MQ=255)
      aCTGTTCCTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCa                                                                                           <  1:6671724/101‑1 (MQ=255)
      aCTGTTCCTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCa                                                                                           >  2:6672269/1‑101 (MQ=255)
         gTTCCTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGaa                                                                                        <  2:6672094/101‑1 (MQ=255)
          ttCCTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGaaa                                                                                       <  2:6671970/101‑1 (MQ=255)
           tCCTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGaaaa                                                                                      >  1:6672271/1‑101 (MQ=255)
            ccTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACAAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAAt                                                                                     <  1:6671825/101‑1 (MQ=255)
             cTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATc                                                                                    >  1:6672276/1‑101 (MQ=255)
               ttAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATcac                                                                                  >  2:6672278/1‑101 (MQ=255)
                 aaCCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCATAAAATcacca                                                                                >  1:6672281/1‑101 (MQ=255)
                    cATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGa                                                                             >  2:6672285/1‑101 (MQ=255)
                     aTCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGaa                                                                            <  2:6671801/101‑1 (MQ=255)
                     aTCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGaa                                                                            >  1:6672287/1‑101 (MQ=255)
                     aTCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGaa                                                                            <  1:6671820/101‑1 (MQ=255)
                      tCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAAt                                                                           <  2:6671960/101‑1 (MQ=255)
                      tCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAAt                                                                           >  2:6672292/1‑101 (MQ=255)
                      tCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAAt                                                                           <  1:6671939/101‑1 (MQ=255)
                         gCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATAcc                                                                        >  1:6672296/1‑101 (MQ=255)
                         gCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATAcc                                                                        <  2:6671912/101‑1 (MQ=255)
                           aaTGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAg                                                                      >  1:6672299/1‑101 (MQ=255)
                            aTGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGa                                                                     >  2:6672300/1‑101 (MQ=255)
                             tGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGaa                                                                    <  2:6671746/101‑1 (MQ=255)
                              ggCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAAt                                                                   >  1:6672302/1‑101 (MQ=255)
                               gCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTcc                                                                                               <  1:6671906/72‑1 (MQ=255)
                               gCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATc                                                                  <  2:6671848/101‑1 (MQ=255)
                               gCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATc                                                                  <  2:6672259/101‑1 (MQ=255)
                                ccAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATcg                                                                 <  2:6672022/101‑1 (MQ=255)
                                 cAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATcgc                                                                <  2:6671862/101‑1 (MQ=255)
                                 cAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATcgc                                                                >  2:6672305/1‑101 (MQ=255)
                                  aaCGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATcgcg                                                               >  2:6672307/1‑101 (MQ=255)
                                   aCGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCa                                                                                           <  2:6671817/72‑1 (MQ=255)
                                    cGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGc                                                             >  2:6672311/1‑101 (MQ=255)
                                    cGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGc                                                             <  2:6671940/101‑1 (MQ=255)
                                     ggTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCa                                                            >  2:6672315/1‑101 (MQ=255)
                                     ggTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCa                                                            <  1:6671952/101‑1 (MQ=255)
                                      gTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCAc                                                           >  2:6672317/1‑101 (MQ=255)
                                       tttACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACt                                                          >  1:6672318/1‑101 (MQ=255)
                                        ttACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTg                                                         <  1:6671933/101‑1 (MQ=255)
                                           cAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGAcc                                                      <  1:6671669/101‑1 (MQ=255)
                                            agagCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCa                                                     >  2:6672326/1‑101 (MQ=255)
                                            agagCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCa                                                     >  1:6672323/1‑101 (MQ=255)
                                             gagCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAg                                                    >  1:6672328/1‑101 (MQ=255)
                                               gCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGaa                                                  <  1:6672010/101‑1 (MQ=255)
                                                 ataataCGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATcaccac                                                                               <  1:6671786/70‑1 (MQ=255)
                                                 ataataCGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGaaaa                                                >  2:6672332/1‑101 (MQ=255)
                                                 ataataCGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGaaaa                                                >  1:6672334/1‑101 (MQ=255)
                                                  taataCGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAAt                                               >  2:6672338/1‑101 (MQ=255)
                                                  taataCGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAAt                                               <  1:6671845/101‑1 (MQ=255)
                                                   aataCGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATg                                              >  1:6672341/1‑101 (MQ=255)
                                                    ataCGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGa                                             <  2:6671870/101‑1 (MQ=255)
                                                     taCGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGAt                                            >  2:6672343/1‑101 (MQ=255)
                                                     taCGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGAt                                            <  1:6671981/101‑1 (MQ=255)
                                                        gAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGAc                                                       >  2:6672346/1‑87 (MQ=255)
                                                        gAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATaaa                                         <  1:6672027/101‑1 (MQ=255)
                                                         aaTAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAAc                                        <  2:6671890/101‑1 (MQ=255)
                                                         aaTAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAAc                                        <  1:6672213/101‑1 (MQ=255)
                                                         aaTAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAAc                                        <  1:6672052/101‑1 (MQ=255)
                                                            aGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAgg                                     >  1:6672349/1‑101 (MQ=255)
                                                             gCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACa                                       >  1:6672352/1‑98 (MQ=255)
                                                             gCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGa                                    <  2:6671918/101‑1 (MQ=255)
                                                              ccGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGaa                                   >  2:6672354/1‑101 (MQ=255)
                                                              ccGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGaa                                   <  1:6671908/101‑1 (MQ=255)
                                                               cGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGaaa                                  <  2:6671986/101‑1 (MQ=255)
                                                                 tAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAaca                                <  2:6671777/101‑1 (MQ=255)
                                                                  agagCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAacac                               >  2:6672358/1‑101 (MQ=255)
                                                                   gagCCAGCCAGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATcg                                                                 <  1:6671900/66‑1 (MQ=255)
                                                                     gccagccaGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGaa                                   <  1:6671805/94‑1 (MQ=255)
                                                                     gccagccaGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACggg                            <  1:6671974/101‑1 (MQ=255)
                                                                     gccagccaGTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACggg                            >  1:6672361/1‑101 (MQ=255)
                                                                       cagccagTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGc                          >  2:6672362/1‑101 (MQ=255)
                                                                       cagccagTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGc                          >  1:6672364/1‑101 (MQ=255)
                                                                       cagccagTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGc                          <  2:6672054/101‑1 (MQ=255)
                                                                        agccagTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCg                         >  1:6672367/1‑101 (MQ=255)
                                                                            agTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTc                     <  1:6671926/101‑1 (MQ=255)
                                                                            agTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTc                     >  1:6672370/1‑101 (MQ=255)
                                                                             gTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACgggg                           <  2:6671927/94‑1 (MQ=255)
                                                                             gTTACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCa                    <  2:6672038/101‑1 (MQ=255)
                                                                             gTTACGGTGGTCGCGATGGTGCGTCCACCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCa                    <  1:6672092/101‑1 (MQ=255)
                                                                              ttACGGTGGTCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCAc                   >  1:6672373/1‑101 (MQ=255)
                                                                                  ggtggtCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTg               <  2:6672006/101‑1 (MQ=255)
                                                                                  ggtggtCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTg               <  1:6672089/101‑1 (MQ=255)
                                                                                  ggtggtCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTg               >  2:6672376/1‑101 (MQ=255)
                                                                                    tggtCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCt             <  2:6672048/101‑1 (MQ=255)
                                                                                    tggtCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCt             >  1:6672378/1‑101 (MQ=255)
                                                                                      gtCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCttt           <  2:6672131/101‑1 (MQ=255)
                                                                                      gtCGCGATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCttt           <  1:6671868/101‑1 (MQ=255)
                                                                                        cgcgATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTtct         >  1:6672382/1‑101 (MQ=255)
                                                                                         gcgATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTtctc        <  1:6672000/101‑1 (MQ=255)
                                                                                         gcgATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTtctc        >  2:6672383/1‑101 (MQ=255)
                                                                                          cgATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCgg                                                              <  1:6671743/46‑1 (MQ=255)
                                                                                          cgATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGaaa                                  <  1:6671888/74‑1 (MQ=255)
                                                                                          cgATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTcc       >  1:6672385/1‑101 (MQ=255)
                                                                                           gATGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCt      <  2:6672067/101‑1 (MQ=255)
                                                                                             tGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATcgcg                                                               <  2:6672079/42‑1 (MQ=255)
                                                                                             tGGTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATaa                                          <  1:6671880/63‑1 (MQ=255)
                                                                                              ggTGCGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTcc       >  1:6672390/1‑97 (MQ=255)
                                                                                                  |                                                                                                 
CCCACAACTGTTCCTTTAACCATCCGCAATGGCCAACGGTTTACAGAGCATAATACGAATAGCCGTAGAGCCAGCCAGTTACGGTGGTCGCGATGGTGTGTCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAG  >  NC_000913/4500613‑4500808

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