Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,574,259	A→G	pseudogene (341/2937 nt)	yjiV →	pseudogene, conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,574,259	0	A	G	100.0%	247.0 / NA	66	pseudogene (341/2937 nt)	yjiV	pseudogene, conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/31);  total (35/31)

TTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGG  >  NC_000913/4574127‑4574394                                                                                                                                      |                                                                                                                                        ttGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATaa                                                                                                                                   >  2:841605/1‑139 (MQ=255)            tatgtgAACCAGCGACAAAATTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATg                                                                                                                        >  2:523614/1‑139 (MQ=255)              tgtgAACCAGCGACAAAATTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGaa                                                                                                                      >  1:144656/1‑139 (MQ=255)              tgtgAACCAGCGACAAAATTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGaa                                                                                                                      >  2:371047/1‑139 (MQ=255)              tgtgAACCAGCGACAAAATTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGaa                                                                                                                      <  1:249711/139‑1 (MQ=255)              tgtgAACCAGCGACAAAATTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGaa                                                                                                                      >  1:654049/1‑139 (MQ=255)                   aCCAGCGACAAAATTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATtctc                                                                                                                 >  2:468245/1‑139 (MQ=255)                    ccAGCGACAAAATTGTGACTGTGGAAGCCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATtctc                                                                                                                 >  2:855229/1‑139 (MQ=255)                    ccAGCGACAAAATTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTctct                                                                                                                >  2:384095/1‑139 (MQ=255)                    ccAGCGACAAAATTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTctct                                                                                                                >  1:699546/1‑139 (MQ=255)                     cAGCGACAAAATTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCtt                                                                                                               >  2:287640/1‑139 (MQ=255)                         gACAAAATTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCa                                                                                                           <  1:1019723/139‑1 (MQ=255)                            aaaaTTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGa                                                                                                        <  1:655335/139‑1 (MQ=255)                              tattgtgattgtgGAAG‑CCCTGTATACGAAATTGCATACTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAATGTGGCGATGAATTCTCTTTGCAGGATg                                                                                                      <  2:468914/138‑1 (MQ=255)                                    gACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTAc                                                                                                <  2:699546/139‑1 (MQ=255)                                    gACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTAc                                                                                                <  1:893643/139‑1 (MQ=255)                                    gACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTAc                                                                                                <  1:217976/139‑1 (MQ=255)                                        gtgGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAACGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGt                                                                                            >  1:689199/1‑139 (MQ=255)                                         tgGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGtt                                                                                           <  2:1001553/139‑1 (MQ=255)                                          ggAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTg                                                                                          <  1:515036/139‑1 (MQ=255)                                          ggAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTg                                                                                          <  1:502896/139‑1 (MQ=255)                                                cccTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACAt                                                                                     <  1:384095/139‑1 (MQ=255)                                                       taCGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAATAGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGGAGTACCTGTTGAACATGACGCTa                                                                              <  2:755117/139‑1 (MQ=255)                                                               tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCa                                                                                                           <  1:488120/102‑1 (MQ=255)                                                               tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCa                                                                                                           >  2:488120/1‑102 (MQ=255)                                                               tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGtt                                                                                           <  2:219737/118‑1 (MQ=255)                                                               tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGtt                                                                                           >  1:219737/1‑118 (MQ=255)                                                               tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTa                                                                              <  2:111941/131‑1 (MQ=255)                                                               tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTa                                                                              >  1:111941/1‑131 (MQ=255)                                                                gCATTCTGTAATGATTGTAATGAGCCTCATCTTCTCGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTCGCGATGAATTCCCTTTGCAGGATGAAGTACCTGTTGAACATGCCGCTACAGAAGaaa                                                                     >  2:61986/1‑139 (MQ=255)                                                                        tAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGaa                                                             <  1:523614/139‑1 (MQ=255)                                                                              ttGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAagag                                                       <  2:343123/139‑1 (MQ=255)                                                                                gTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTGTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTAc                                                                             >  2:802484/1‑115 (MQ=255)                                                                                gTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTAc                                                                             <  1:802484/115‑1 (MQ=255)                                                                                  aaTGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACa                                                   >  2:705309/1‑139 (MQ=255)                                                                                  aaTGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACa                                                   >  1:576814/1‑139 (MQ=255)                                                                                        ccgcATCTTCTCGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAACGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTc                                             >  2:651340/4‑139 (MQ=255)                                                                                        ccTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGaa                                                                                        <  2:595348/96‑1 (MQ=255)                                                                                        ccTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGaa                                                                                        >  1:595348/1‑96 (MQ=255)                                                                                        ccTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTc                                             >  2:289982/1‑139 (MQ=255)                                                                                        ccTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTc                                             >  1:499423/1‑139 (MQ=255)                                                                                        ccTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTc                                             >  1:1115928/1‑139 (MQ=255)                                                                                           cATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCg                                                               >  2:581656/1‑118 (MQ=255)                                                                                           cATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCg                                                               <  1:581656/118‑1 (MQ=255)                                                                                            aTCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGcc                                         >  2:831432/1‑139 (MQ=255)                                                                                             tcttctGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCt                                        <  2:116865/139‑1 (MQ=255)                                                                                                  tGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGtt                                   >  2:10925/1‑139 (MQ=255)                                                                                                       cGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTa                                                                                                 <  2:820773/72‑1 (MQ=255)                                                                                                       cGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTa                                                                                                 >  1:820773/1‑72 (MQ=255)                                                                                                         ggACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACAt                                                                                     <  1:1111323/82‑1 (MQ=255)                                                                                                         ggACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACAt                                                                                     >  2:1111323/1‑82 (MQ=255)                                                                                                          gACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAgg                                                                                                         <  1:848178/61‑1 (MQ=255)                                                                                                          gACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAgg                                                                                                         >  2:848178/1‑61 (MQ=255)                                                                                                          gACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGcc                                         >  1:705565/1‑125 (MQ=255)                                                                                                          gACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGcc                                         <  2:705565/125‑1 (MQ=255)                                                                                                             aaaaaGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCg                        <  1:855229/139‑1 (MQ=255)                                                                                                                    gCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGgaga                 >  2:937824/1‑139 (MQ=255)                                                                                                                    gCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGgaga                 >  2:677294/1‑139 (MQ=255)                                                                                                                    gCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGgaga                 >  2:102060/1‑139 (MQ=255)                                                                                                                          tAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAg           <  1:287640/139‑1 (MQ=255)                                                                                                                            gTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCg         <  1:289982/139‑1 (MQ=255)                                                                                                                             tCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGa        >  1:854677/1‑139 (MQ=255)                                                                                                                               cAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAgg      <  1:532755/139‑1 (MQ=255)                                                                                                                               cAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAgg      <  2:854677/139‑1 (MQ=255)                                                                                                                                aGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGc     <  1:468245/139‑1 (MQ=255)                                                                                                                                  tGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAg   >  2:35399/1‑139 (MQ=255)                                                                                                                                   ggggAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGATCATGACGCTACATAAGATAAAGACGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAgg  <  2:3460/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        TTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAG‑CCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGG  >  NC_000913/4574127‑4574394

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