Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,569,309	G→A	pseudogene (312/921 nt)	yjiP →	pseudogene, transposase_31 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,569,309	0	G	A	100.0%	220.0 / NA	68	pseudogene (312/921 nt)	yjiP	pseudogene, transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/32);  total (36/32)

CGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACGGAAG  >  NC_000913/4569173‑4569444                                                                                                                                         |                                                                                                                                        cgcTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAacg                                                                                                                                       <  1:167260/139‑1 (MQ=255)   cTACAGTCCGCTATTTTAGAGTTGGTCAAGAGCCGCGAAGGTGATGGCTATATCTATGTGGTGATTGAAGATCAGAGCCGCGATGGCATTCATACGGCCTTTCGCCTGATGCGGTATTCCATGGCCGTGATGTAACGcc                                                                                                                                     <  2:636496/139‑1 (MQ=255)         tCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATag                                                                                                                               <  1:298874/139‑1 (MQ=255)             atatTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAacg                                                                                                                                       >  1:993607/1‑127 (MQ=255)             atatTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCa                                                                                                                           <  1:331680/139‑1 (MQ=255)                   tATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATaa                                                                                                                     >  1:840624/1‑139 (MQ=255)                     tGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAAc                                                                                                                   >  1:172685/1‑139 (MQ=255)                       gTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAacg                                                                                                                                       >  1:215814/1‑117 (MQ=255)                       gTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAacg                                                                                                                                       <  2:215814/117‑1 (MQ=255)                         cGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCa                                                                                                               <  1:145310/139‑1 (MQ=255)                         cGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGAAACGCCATATAGAGCATGATAAACGCCa                                                                                                               >  1:633455/1‑139 (MQ=255)                          ggTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAg                                                                                                              >  2:828467/1‑139 (MQ=255)                             aaaGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                           <  1:1049392/139‑1 (MQ=255)                                    cgcgAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGt                                                                                                    >  1:475637/1‑139 (MQ=255)                                             gATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                           <  2:992613/139‑1 (MQ=255)                                                           tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCAtata                                                                                                                                >  2:533653/1‑88 (MQ=255)                                                           tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCAtata                                                                                                                                <  1:533653/88‑1 (MQ=255)                                                           tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCa                                                                                                               >  1:312789/1‑105 (MQ=255)                                                           tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCa                                                                                                               <  2:312789/105‑1 (MQ=255)                                                           tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGAt                                                                                       <  2:154107/129‑1 (MQ=255)                                                           tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGAt                                                                                       >  1:154107/1‑129 (MQ=255)                                                            gtggtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGa                                                                                        <  1:692754/127‑1 (MQ=255)                                                            gtggtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGa                                                                                        >  2:692754/1‑127 (MQ=255)                                                            gtggtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATc                                                                            >  1:335770/1‑139 (MQ=255)                                                               gtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACg                                                                         >  2:834128/1‑139 (MQ=255)                                                                tgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACgg                                                                        >  2:813248/1‑139 (MQ=255)                                                                 gATTGAACATCAGAGCCGCGAGTACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGt                                                                       >  1:483289/1‑139 (MQ=255)                                                                 gATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGt                                                                       >  1:385362/1‑139 (MQ=255)                                                                  aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCg                                                                                         >  1:319049/1‑120 (MQ=255)                                                                  aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCg                                                                                         <  2:319049/120‑1 (MQ=255)                                                                  aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTa                                                                      >  1:660121/1‑139 (MQ=255)                                                                  aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTa                                                                      >  1:75838/1‑139 (MQ=255)                                                                  aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTa                                                                      >  2:505101/1‑139 (MQ=255)                                                                  aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTa                                                                      >  1:928671/1‑139 (MQ=255)                                                                      aaCATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCg                                                                  <  1:828467/139‑1 (MQ=255)                                                                             gagCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                           <  2:385362/139‑1 (MQ=255)                                                                              agCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                             >  2:53398/1‑120 (MQ=255)                                                                              agCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                             <  1:53398/120‑1 (MQ=255)                                                                              agCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCtt                                                          >  2:338206/1‑139 (MQ=255)                                                                                ccGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTAc                                                        >  1:584282/1‑139 (MQ=255)                                                                                  gcgAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCg                                                                  >  2:586689/1‑127 (MQ=255)                                                                                  gcgAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCg                                                                  <  1:586689/127‑1 (MQ=255)                                                                                   cgAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTa                                                                >  1:1027823/1‑128 (MQ=255)                                                                                   cgAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTa                                                                <  2:1027823/128‑1 (MQ=255)                                                                                         cATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTcc                                               <  1:338206/139‑1 (MQ=255)                                                                                         cATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTcc                                               <  1:505101/139‑1 (MQ=255)                                                                                          aTTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTccc                                              >  1:517654/1‑139 (MQ=255)                                                                                                 tGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTg                                                   <  2:145520/127‑1 (MQ=255)                                                                                                 tGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTg                                                   >  1:145520/1‑127 (MQ=255)                                                                                                     cTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGcc                                                                                                                                     <  2:751883/41‑1 (MQ=255)                                                                                                     cTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGcc                                                                                                                                     >  1:751883/1‑41 (MQ=255)                                                                                                      tttCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCa                                                                                                               <  2:629321/62‑1 (MQ=255)                                                                                                      tttCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCa                                                                                                               >  1:629321/1‑62 (MQ=255)                                                                                                       ttCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                           <  2:563713/113‑1 (MQ=255)                                                                                                       ttCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                           >  1:563713/1‑113 (MQ=255)                                                                                                        tCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                           <  2:870109/64‑1 (MQ=255)                                                                                                        tCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                           >  1:870109/1‑64 (MQ=255)                                                                                                                 gcgcTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGcc                       <  1:813248/139‑1 (MQ=255)                                                                                                                            tGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTAc            <  2:584282/139‑1 (MQ=255)                                                                                                                            tGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTAc            >  2:841852/1‑139 (MQ=255)                                                                                                                               cGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTa                                                                <  2:929959/84‑1 (MQ=255)                                                                                                                               cGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTa                                                                >  1:929959/1‑84 (MQ=255)                                                                                                                                 gTGATGTTACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAc       >  1:918243/1‑139 (MQ=255)                                                                                                                                 gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCa                                                                                              <  1:729608/52‑1 (MQ=255)                                                                                                                                 gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCa                                                                                              >  2:729608/1‑52 (MQ=255)                                                                                                                                 gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAc       >  2:273882/1‑139 (MQ=255)                                                                                                                                 gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAc       >  2:367603/1‑139 (MQ=255)                                                                                                                                    aTGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACGGa    <  1:841852/139‑1 (MQ=255)                                                                                                                                    aTGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACGGa    <  2:413830/139‑1 (MQ=255)                                                                                                                                    aTGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACGGa    <  1:3562/139‑1 (MQ=255)                                                                                                                                      gTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACGGAAg  <  1:273882/139‑1 (MQ=255)                                                                                                                                      gTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACGGAAg  <  1:936836/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        CGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACGGAAG  >  NC_000913/4569173‑4569444

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