Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,060,264	A→C	6.0%	intergenic (+34/‑183)	bipA → / → yihL	ribosome‑dependent GTPase, ribosome assembly factor/putative transcriptional regulator YihL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,060,264	0	A	C	6.0%	213.7 / 5.0	89	intergenic (+34/‑183)	bipA/yihL	ribosome‑dependent GTPase, ribosome assembly factor/putative transcriptional regulator YihL
Reads supporting (aligned to +/- strand):  ref base A (37/44);  new base C (3/3);  total (40/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.76e-01

CGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAA  >  NC_000913/4060122‑4060404                                                                                                                                               |                                                                                                                                              cGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGa                                                                                                                                               <  1:573767/143‑1 (MQ=255)     gatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGc                                                                                                                                           <  1:331142/143‑1 (MQ=255)     gatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                           >  2:68558/1‑143 (MQ=255)      atgaCGAACTGGTAGAAGCGCCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                                                          <  1:37011/143‑1 (MQ=255)       tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                                                                         <  2:52188/143‑1 (MQ=255)       tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATGCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAg                                                                                                                                         <  1:430223/143‑1 (MQ=255)        gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAACACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTCCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                                                                        >  2:699460/1‑143 (MQ=255)        gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGCCACCTGACGGGAAACGATCGTCGCCGCGCCAACCGCGCACCGAACGACGATTAATTCCGTTTTTTAGTCATAAAAAACCTCCCAGCGCTGGCAgg                                                                                                                                        >  2:175075/1‑143 (MQ=255)            aaCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACATGCCAGCGATGGCAGG‑tttt                                                                                                                                   >  2:725270/1‑143 (MQ=255)             aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTACACGTCACCTGACGGAAAACGCTCGTCGCCGGGCCCACCGCGCACCGAAAGACGCTTAATCCCGTTTTTTGGTCATGTAACACCTGCCAGCGATGGCAGG‑ttttt                                                                                                                                  >  2:314157/1‑143 (MQ=255)               tGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttt                                                                                                                                >  1:731387/1‑143 (MQ=255)               tGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttt                                                                                                                                >  1:759781/1‑143 (MQ=255)                gggataaGTTACTCCGACCTCTATCCGTAGTCGTAAACGTCACCTGACGGGAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttat                                                                                                                               <  1:560547/138‑3 (MQ=255)                 gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAGCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAACACACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                                                              >  1:622962/1‑143 (MQ=255)                  tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTg                                                                                                                             <  2:731387/143‑1 (MQ=255)                     aaGTGACTCCGCCCTCTATCCGTATGCTCAAACGTCACATGACGGAAAGCGATCGTCGCCGCGCCAACCGCTAACCGAAAGACGATTAATTTCTTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCt                                                                                                                          <  2:444628/143‑1 (MQ=255)                        tgccTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTc                                                                                                                       <  1:549359/140‑1 (MQ=255)                          aCTCCGCCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAACCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAACACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGt                                                                                                                     <  1:38453/143‑1 (MQ=255)                            tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAACAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGa                                                                                                                   >  2:224852/1‑143 (MQ=255)                            tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGCTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TCTTTTTTTGCCTGTGGTGa                                                                                                                   >  2:452300/1‑143 (MQ=255)                               gACCTCTAGCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAgcg                                                                                                                <  1:461254/143‑1 (MQ=255)                                   tctATCCGTATTCGTAAACGTCACCTGACGGAACACGATCGTAGCCGCGCCAACCGCGCACCGAAAGACGCTTACTTTCGTTTTTTAGTCATAAAAAACATGCCAGCGATGGCAGG‑CTTATTTTTGCCTCTCGTGAGCGgacg                                                                                                            >  2:99806/1‑142 (MQ=255)                                      aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGGCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCa                                                                                                         <  1:652429/143‑1 (MQ=255)                                         cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                                                                      >  2:17930/1‑143 (MQ=255)                                         cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATCAATCTCGTTTTTTCGTCATAACAAACCTGCCAGCGAAGGCCGG‑TATCTTTTTGCCTGTCGTGAGCGCACCCCAGTg                                                                                                      >  2:168549/1‑143 (MQ=255)                                            ggtgggaaaCGTCACCTGAGGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTCCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                   <  2:353496/137‑1 (MQ=255)                                            atccGTAACCGACACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAGAGATGATTATTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                   <  1:223716/140‑1 (MQ=255)                                            aTTCGTACACGTCACCCGACTGAAAACGATCGTCTCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                   <  2:546891/143‑1 (MQ=255)                                            aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCTCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                   <  2:797882/143‑1 (MQ=255)                                            aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                   >  1:591470/1‑143 (MQ=255)                                            aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                   <  1:215391/143‑1 (MQ=255)                                             tgcGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                                                                  <  2:512518/141‑1 (MQ=255)                                             tccGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                                                                  <  2:645961/141‑1 (MQ=255)                                                gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                                                               >  1:670569/1‑143 (MQ=255)                                                 tAAACGTCACCTGACGGTAAACGATCGTCGCCGCCCCAACCGCCCCCCGAAAGACGCTTACTTTCGGTTTTTCGTCAGCAAAACCCTGCCAGCGATGGCAGG‑TTTTGTTTTGCGTGGCGTGACCCCACACCCGTGCgcgcacc                                                                                              >  2:545971/1‑137 (MQ=255)                                                  aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCACCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                                                             <  1:682298/143‑1 (MQ=255)                                                  aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTGATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACCCACATGAGAGTAtct                                                                                             >  2:277858/1‑143 (MQ=255)                                                     cggcACCTGGCGGAAAACGATCGTCGCAGCACCAACAGCGCATCGAAAGAAGATAAATATAGTATTTCAGTCAGAAAAAACCGGCCAGCGATGGCAGG‑TATTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                                                          <  2:139454/140‑1 (MQ=255)                                                          ccTGACGGAAAACGAGCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTACGTTTTTTAGTCATAAAAAACCTGCCCGCGATGGCAGG‑TTTTTTTTTGCCTGTCGCGAGCGCACACCAATGCGAGCATCTCCCTtatc                                                                                     >  2:142377/1‑142 (MQ=255)                                                           cTGACGGAATACGATTGTCGCCGCGCCGACCGCTCACCGAAAGACGATTTATTTCGTTTTTTAGTCATAAAAAGCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCGTATaa                                                                                    <  2:507793/143‑1 (MQ=255)                                                              aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                                                                 >  1:252719/1‑143 (MQ=255)                                                               cGGAAATCGATCGTCGCCGCGCCAACCGCGCACCGACAGACGATTAATTTCGTTTTTTACTCATAAACAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                                                                <  2:62692/143‑1 (MQ=255)                                                                ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTTTTGTTAGCGCACAACAATACGAGTATCGCCCTTTCAAGtata                                                                               >  2:721557/1‑140 (MQ=255)                                                                ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTGTTTGTCGGTCGTGAGGGCACACCAATGCGAGTATCTCCCTTATAAGTCtg                                                                               >  2:617444/1‑143 (MQ=255)                                                                  aaaaCGATCGTCGCCGCGCCAACCTCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                             <  1:173214/143‑1 (MQ=255)                                                                  aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                             >  1:254265/1‑143 (MQ=255)                                                                  aaaaCGATCGTCCCCGGGTCCACCGCGCACCGAAAGACGATTAGTTGCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                             <  2:559946/143‑1 (MQ=255)                                                                   acaCGCTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAGACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                             <  2:706562/141‑1 (MQ=255)                                                                   aaaCGATCGGCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTGTAGTCATAAAATAGCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGGGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                            <  2:531952/143‑1 (MQ=255)                                                                    accGATCGTCGCCGCGCCAGCCGCGCACCGTGAGACTATTAATTTCGTTTTTTAGTCATAAAATACCTGCCAGCGGTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGTGAGTATCTCCCTTATAAGTCTGTGGt                                                                           <  2:629467/141‑1 (MQ=255)                                                                    aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGCGTATATCCCTTATAAGTCTGTGGt                                                                           >  1:510191/1‑143 (MQ=255)                                                                           gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAAATTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑CTTTTTTTTGCCTGTCGTGCGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                                                                    >  1:672127/1‑143 (MQ=255)                                                                            tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATACAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCTTGTCGTGAGCGCACACAAATGAGAGTATCTCCCTTATAAGTCTGTGGTTACGTcca                                                                   >  1:99401/1‑141 (MQ=255)                                                                             ggcggcgcCAACCGCGCATCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCGGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                                                                  <  2:76201/139‑1 (MQ=255)                                                                             cgtcgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAATTACCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                                                                  <  2:40841/140‑1 (MQ=255)                                                                               ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                                                                >  1:360910/1‑143 (MQ=255)                                                                                 gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTGTAGTCATAAAAAACCTGCAAGCGATGGCAGG‑TTTTTTTTTGACTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTAAAAGTCTGTGGTTACGTCAATagag                                                              <  1:44770/143‑1 (MQ=255)                                                                                    ccAACCGCGCACCGAAAGACGATTACTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGCGCGCACACCAATGCGAGTCTCTCCCTTATACGTCCGTGGTTCCGCCCATACAgcgc                                                           >  2:155173/1‑143 (MQ=255)                                                                                       aCCGCGCACCGAAAGACGATTAATTACGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtct                                                        >  1:361275/1‑141 (MQ=255)                                                                                         cgcgCACCGAAAGACGATTACTTTCGTTTTTTAGTCATAAAAAGCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTCTAATTCTATGGTTACGTCAATAGAGAGCTTTGt                                                      <  2:253130/143‑1 (MQ=255)                                                                                          gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTGTGTc                                                     >  1:102590/1‑143 (MQ=255)                                                                                           cgcACCGAGAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca                                                    <  2:689294/143‑1 (MQ=255)                                                                                                cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt                                               <  1:231702/143‑1 (MQ=255)                                                                                                cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt                                               >  1:453701/1‑143 (MQ=255)                                                                                                  aaaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTAt                                             <  1:805624/143‑1 (MQ=255)                                                                                                    aGACGCTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTGTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCCCCATTATCt                                           <  2:629478/143‑1 (MQ=255)                                                                                                      aCGATTAATTTCTTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGc                                         <  2:750951/143‑1 (MQ=255)                                                                                                      aCGATTAATTTCTTTGTTTATTCACAAGAAACCTGCCAGCGATCGCAGG‑TTTTTGTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGc                                         <  2:666183/143‑1 (MQ=255)                                                                                                      aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGc                                         <  1:104581/143‑1 (MQ=255)                                                                                                      aCGATTAATTACAATTCTAAGTCATGCAGAACCTGCCAGCGATCGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGc                                         <  1:553419/143‑1 (MQ=255)                                                                                                               ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATa                                <  1:343007/143‑1 (MQ=255)                                                                                                                  gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACATTAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGc                             >  2:220672/1‑143 (MQ=255)                                                                                                                   ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCTTGTCGTGAGCGCACACCAATGAGAGTATCTCCCTTATACGTCTGTGGTTACGTCAGTAGAGCGCTTTGTCACAATTATCTGCAAAGTCATACGcc                            >  1:786726/1‑143 (MQ=255)                                                                                                                   ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTCTCACCATTATCTGCAAATTCATACGcag                           >  1:105354/1‑141 (MQ=255)                                                                                                                   ttttttAGTCATAAAAAACCTGCCAGAGATGGCAGG‑TTTTTTTTTGCCTGTAGTGAGCGCACACCAATGAGAGTATCTTCCGTATAAGTCTGTGGTTAAGTAAATAGAGATATTTGTCACACTTATCTGCTAAGTCCTaatcc                            >  2:349972/1‑139 (MQ=255)                                                                                                                     gtttAGTCAGAAAGAGCGTGCCAGAGAGGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGt                          <  1:720444/142‑1 (MQ=255)                                                                                                                        tAGTCATAACACCCCTGCCCTCGTTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTaa                       <  1:731206/143‑1 (MQ=255)                                                                                                                           tCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAAtt                     >  1:505766/1‑143 (MQ=255)                                                                                                                             aaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCAAACCCATGCGAGTATCGTCCTTATAAGTCTGTGGTTACTTCAATAGAGAGTCTTGTCCCCATTAGCTGCAAAGCTCTACGCCGTTTATTGcg                  >  2:309492/3‑142 (MQ=255)                                                                                                                              tAAAAAACCTTCCAGCGATGGGAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctt                 <  1:68558/143‑1 (MQ=255)                                                                                                                                aaagaCCTGCCAGCGATGGCAGG‑TTCCTTTTGGCCTGTCGTGAGCGCACACCAATACGAGTATCTCGCAGATAAGTCTGTGGGTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttc               <  2:14883/143‑1 (MQ=255)                                                                                                                                    gccTGCCAGCGAGTGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttt           <  2:104686/142‑1 (MQ=255)                                                                                                                                    aCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGACTTTCTCCCTTATAAGTCTGTGGTTACGTCGATAGAAAGCTTAGTAACAATTATCTGGAAACTCATACGACGTACGTTGgtttctcta           >  2:509452/1‑140 (MQ=255)                                                                                                                                       tGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCCCAAACCTGAATCGCTCTCTGTGATAAGTCTGTGGTGACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg        <  1:186447/143‑1 (MQ=255)                                                                                                                                       tGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg        >  1:501807/1‑143 (MQ=255)                                                                                                                                       tGCCAGCGATGGCAGG‑TTTTTTCTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg        <  1:17930/143‑1 (MQ=255)                                                                                                                                         ccAGCGATGGCGGG‑TTTTGTCTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc      <  2:520132/143‑1 (MQ=255)                                                                                                                                         ccAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTACTAAGTCAGGGGTTACGTCAATAGAGAGCATTGTAACAATTATCTGCAAAGACATACGCCGTTAATTGCTTTCTTTTTTGgc      >  1:366774/1‑143 (MQ=255)                                                                                                                                             cGATGGCAGG‑TTTTTTTTTTTTTGTCGTGAGCGCACACCAATGCGAGTATCTCCGTTATTAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaa  >  2:41330/1‑143 (MQ=255)                                                                                                                                             cGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaa  <  1:62823/143‑1 (MQ=255)                                                                                                                                               |                                                                                                                                              CGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAA  >  NC_000913/4060122‑4060404

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