Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,618,281	Δ1 bp	coding (628/666 nt)	marC ←	UPF0056 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,618,281	0	T	.	100.0%	316.2 / NA	65	coding (628/666 nt)	marC	UPF0056 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (37/28);  total (37/28)

GGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGA  >  NC_000913/1618137‑1618421                                                                                                                                                 |                                                                                                                                             ggTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAAc                                                                                                                                       <  2:25943/151‑1 (MQ=255)   tGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTg                                                                                                                                     >  1:914711/1‑151 (MQ=255)     ccTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGc                                                                                                                                    >  2:63940/1‑150 (MQ=255)        gCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGcac                                                                                                                                  <  2:729606/149‑1 (MQ=255)         cTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACAcc                                                                                                                               >  2:537105/1‑151 (MQ=255)          ttttGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACAccc                                                                                                                              >  1:794922/1‑151 (MQ=255)           tttGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACAccc                                                                                                                              <  2:721720/150‑1 (MQ=255)              gCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcatac                                                                                                                          >  1:420685/1‑151 (MQ=255)              gCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcata                                                                                                                           <  1:23805/150‑1 (MQ=255)                cgcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcatacat                                                                                                                        >  1:306927/1‑151 (MQ=255)                cgcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcatacat                                                                                                                        <  2:787353/151‑1 (MQ=255)                 gcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcatacat                                                                                                                        >  2:433683/1‑150 (MQ=255)                     tAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACca                                                                                                                    <  2:95255/150‑1 (MQ=255)                      aaTTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACcagc                                                                                                                  <  1:544528/151‑1 (MQ=255)                      aaTTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACcag                                                                                                                   <  2:706069/150‑1 (MQ=255)                           ggTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGa                                                                                                              >  2:439093/1‑150 (MQ=255)                           ggTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGa                                                                                                              >  2:638930/1‑150 (MQ=255)                             tCAATCATTATATTTCGCCGCAGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGCACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGa                                                                                                              >  2:444721/1‑148 (MQ=255)                              cAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATcc                                                                                                          >  2:76679/1‑151 (MQ=255)                                  cATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATc                                                                                                      <  1:732185/151‑1 (MQ=255)                                     tataTTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATc                                                                                                      >  1:301326/1‑148 (MQ=255)                                       tatTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACg                                                                                                 >  1:790063/1‑151 (MQ=255)                                       tatTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACg                                                                                                 >  2:516015/1‑151 (MQ=255)                                        atTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACg                                                                                                 <  2:212142/150‑1 (MQ=255)                                            cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGgag                                                                                             <  1:301764/150‑1 (MQ=255)                                              ccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACggg                                                                                               >  1:683224/1‑146 (MQ=255)                                               cgcTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCg                                                                                         <  1:456472/151‑1 (MQ=255)                                                gcTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGc                                                                                        >  2:176930/1‑151 (MQ=255)                                                   ggCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCt                                                                                       <  1:302735/149‑1 (MQ=255)                                                    gCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTc                                                                                     >  1:882604/1‑150 (MQ=255)                                                      cAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAAt                                                                                  <  1:49378/151‑1 (MQ=255)                                                      cAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAAt                                                                                  >  2:258658/1‑151 (MQ=255)                                                       aGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAAt                                                                                  >  1:492051/1‑150 (MQ=255)                                                       aGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATg                                                                                 >  1:841946/1‑151 (MQ=255)                                                           acacTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCAc                                                                             <  1:810448/151‑1 (MQ=255)                                                            cacTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCAcc                                                                            >  2:45318/1‑151 (MQ=255)                                                                     aaCAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCccac                                                                   >  2:107927/1‑151 (MQ=255)                                                                       cAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCccacc                                                                  >  1:813535/1‑150 (MQ=255)                                                                        aGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAg                                                                <  2:558275/151‑1 (MQ=255)                                                                           cGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCg                                                             <  1:807987/151‑1 (MQ=255)                                                                              aaTAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCAt                                                          >  2:214456/1‑151 (MQ=255)                                                                                        aGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCg                                                <  1:748419/151‑1 (MQ=255)                                                                                         ggCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCgg                                               <  1:193579/151‑1 (MQ=255)                                                                                         ggCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCgg                                               <  1:278141/151‑1 (MQ=255)                                                                                              aTTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTAGCCCACCAGCCGCATAATTGCGCCGGAGCt                                           >  1:229911/1‑150 (MQ=255)                                                                                                 cGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGt                                       >  1:839898/1‑151 (MQ=255)                                                                                                  gggCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTa                                      <  2:569943/151‑1 (MQ=255)                                                                                                   ggCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCAt                                                          >  2:52855/1‑130 (MQ=255)                                                                                                   ggCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCAt                                                          <  1:52860/130‑1 (MQ=255)                                                                                                        tttttAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTc                                >  1:782615/1‑151 (MQ=255)                                                                                                          tttAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTccc                              >  2:728846/1‑151 (MQ=255)                                                                                                            tAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTcccc                             >  1:493541/1‑150 (MQ=255)                                                                                                             aaTGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCcac                           >  1:850203/1‑151 (MQ=255)                                                                                                              aTGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCcaca                          <  1:737708/151‑1 (MQ=255)                                                                                                                 gTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACaaaa                       <  2:310668/151‑1 (MQ=255)                                                                                                                  tACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACaaaa                       >  2:118245/1‑150 (MQ=255)                                                                                                                       tttAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAAcc                  <  2:19250/150‑1 (MQ=255)                                                                                                                        ttAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGc                >  2:120740/1‑151 (MQ=255)                                                                                                                          aaTGATTTCCAGGATGCCGTTAA‑AATAAACTGCACCCCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGcca              <  1:225169/151‑1 (MQ=255)                                                                                                                          aaTGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGcc               >  2:97007/1‑150 (MQ=255)                                                                                                                             gATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGccacc            >  2:747882/1‑150 (MQ=255)                                                                                                                              aTTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGccacc            >  1:535553/1‑149 (MQ=255)                                                                                                                               tttCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAgg         >  1:343692/1‑151 (MQ=255)                                                                                                                                  ccAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGaa       <  1:336107/150‑1 (MQ=255)                                                                                                                                      gATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGa  <  1:76685/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                             GGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGA  >  NC_000913/1618137‑1618421

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