Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	858,673	A→C	5.1%	L130R (CTC→CGC) ‡	ybiU ←	DUF1479 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	858,673	0	A	C	5.1%	229.3 / 4.4	79	L130R (CTC→CGC) ‡	ybiU	DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base C (2/2);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

ATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAT  >  NC_000913/858525‑858818                                                                                                                                                     |                                                                                                                                                  aTCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGa                                                                                                                                                   <  7:16733/149‑1 (MQ=255)    cGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCAGAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTTCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGc                                                                                                                                                 <  3:207027/148‑1 (MQ=255)    cGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGc                                                                                                                                                 <  2:302082/148‑1 (MQ=255)        ttAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTccc                                                                                                                                            <  5:302555/149‑1 (MQ=255)             ccATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCgaag                                                                                                                                        >  6:212381/1‑148 (MQ=255)              cATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCgaaga                                                                                                                                       >  1:418732/1‑148 (MQ=255)              cATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCgaaga                                                                                                                                       >  5:207462/1‑148 (MQ=255)                    tttCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATc                                                                                                                                >  8:146660/1‑149 (MQ=255)                     ttCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCga                                                                                                                                          >  7:170988/1‑138 (MQ=255)                     ttCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCga                                                                                                                                          <  8:170988/138‑1 (MQ=255)                        cATCACTTGAAAATGTCCACAGCCGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGTTGTGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCg                                                                                                                            <  6:67412/149‑1 (MQ=255)                         aTCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCgg                                                                                                                           >  3:131052/1‑149 (MQ=255)                         aTCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCgg                                                                                                                           >  8:439916/1‑149 (MQ=255)                              tttCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCt                                                                                                                      >  3:215731/1‑149 (MQ=255)                              tttCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCt                                                                                                                      <  6:55163/149‑1 (MQ=255)                              tttCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCt                                                                                                                      >  3:36984/1‑149 (MQ=255)                              tttCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCgagat                                                                                                                                       <  1:501502/132‑4 (MQ=255)                              tttCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCgagat                                                                                                                                       >  2:501502/1‑129 (MQ=255)                                tCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCt                                                                                                                      <  7:146660/147‑1 (MQ=255)                                  aaaTGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTg                                                                                                                   >  1:71716/1‑148 (MQ=255)                                        ccACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCt                                                                                                            >  4:413042/1‑149 (MQ=255)                                          acaGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCg                                                                                                          >  6:481859/1‑149 (MQ=255)                                          acaGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCg                                                                                                          >  7:487952/1‑149 (MQ=255)                                           caGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGt                                                                                                         <  4:215731/149‑1 (MQ=255)                                            aGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTCGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTc                                                                                                        <  2:227467/149‑1 (MQ=255)                                            aGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGt                                                                                                         >  8:32326/1‑148 (MQ=255)                                             gACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTc                                                                                                        >  7:381157/1‑148 (MQ=255)                                                     agaAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGtt                                                                                               <  5:481859/149‑1 (MQ=255)                                                     agaAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGtt                                                                                               <  7:439916/149‑1 (MQ=255)                                                      gaAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTg                                                                                              >  1:129683/1‑149 (MQ=255)                                                           gACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTc                                                                                         <  5:212381/149‑1 (MQ=255)                                                            aCTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCTACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTcc                                                                                        <  4:36984/149‑1 (MQ=255)                                                            aCTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTcc                                                                                        <  8:381157/149‑1 (MQ=255)                                                            aCTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTc                                                                                         >  7:247120/1‑148 (MQ=255)                                                                cgcATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGa                                                                                     <  3:413042/148‑1 (MQ=255)                                                                 gcATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGAGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGa                                                                                     >  6:413216/1‑147 (MQ=255)                                                                       gCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATa                                                                                   <  2:185493/143‑1 (MQ=255)                                                                       gCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATa                                                                                   >  1:185493/1‑143 (MQ=255)                                                                        ccATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAg                                                                             >  1:260294/1‑148 (MQ=255)                                                                          aTTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCAt                                                                          <  6:207462/149‑1 (MQ=255)                                                                          aTTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCAGGGTCATTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCGAGCAt                                                                          >  5:235179/1‑149 (MQ=255)                                                                             tcttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAATTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATAGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCAAGCATCGa                                                                       >  1:113757/1‑149 (MQ=255)                                                                              cttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGAc                                                                      >  7:377089/1‑149 (MQ=255)                                                                               ttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACt                                                                     >  5:58862/1‑149 (MQ=255)                                                                                 cACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGAc                                                                      <  4:131052/146‑1 (MQ=255)                                                                                   ctctGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATc                                                                 >  5:106749/1‑149 (MQ=255)                                                                                        gCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCa                                                              <  2:71716/147‑1 (MQ=255)                                                                                            gCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGAc                                                                      <  1:417705/135‑1 (MQ=255)                                                                                            gCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGAc                                                                      >  2:417705/1‑135 (MQ=255)                                                                                              cTGCATTTGCGCCTGCGCCCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgcg                                                      >  7:243815/1‑149 (MQ=255)                                                                                              cTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGc                                                                                             >  8:187408/1‑110 (MQ=255)                                                                                              cTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGc                                                                                             <  7:187408/110‑1 (MQ=255)                                                                                              cTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgcg                                                      >  3:87777/1‑149 (MQ=255)                                                                                                   tttGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGt                                                 >  7:292618/1‑149 (MQ=255)                                                                                                         ccTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCgg                                           <  7:32326/149‑1 (MQ=255)                                                                                                          cTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCggg                                          >  7:460266/1‑149 (MQ=255)                                                                                                             cGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCggggg                                        <  2:418732/148‑1 (MQ=255)                                                                                                                 cAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATg                                   <  8:487952/149‑1 (MQ=255)                                                                                                                   gTAGATGGGGTAAATCTCGTGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGACGGGGCCTTTTTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGAATGATCCCAGCCTAGCGCTTGTACGCGGGGGCACTGGc                                 >  6:476506/1‑149 (MQ=255)                                                                                                                   gTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGCCCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTa                                                         >  2:118886/1‑125 (MQ=255)                                                                                                                   gTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTa                                                         <  1:118886/125‑1 (MQ=255)                                                                                                                      gATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCt                               <  8:247120/148‑1 (MQ=255)                                                                                                                       aTGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGATATTATCGCCGTGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCTACTGATCCCAGCCTAGCGCTTGt                                                 >  1:37701/1‑129 (MQ=255)                                                                                                                       aTGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGt                                                 <  2:37701/129‑1 (MQ=255)                                                                                                                            gTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGCCTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCa                        >  3:36522/1‑149 (MQ=255)                                                                                                                             tAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCAc                       <  5:413216/149‑1 (MQ=255)                                                                                                                               aaTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCg                     <  2:260294/149‑1 (MQ=255)                                                                                                                                aTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCg                     >  7:14712/1‑148 (MQ=255)                                                                                                                                   tcGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAAt                                    >  7:75094/1‑130 (MQ=255)                                                                                                                                   tcGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAAt                                    <  8:75094/130‑1 (MQ=255)                                                                                                                                     gggACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAAc               >  7:396284/1‑149 (MQ=255)                                                                                                                                             aaGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGa                                                                                     >  4:9202/1‑71 (MQ=255)                                                                                                                                             aaGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGa                                                                                     <  3:9202/71‑1 (MQ=255)                                                                                                                                             aaGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGt       >  4:178452/1‑149 (MQ=255)                                                                                                                                               gcgcTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGtt      >  1:182091/1‑148 (MQ=255)                                                                                                                                                cgcTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTa    <  6:106749/149‑1 (MQ=255)                                                                                                                                                 gcTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTaa   >  5:60475/1‑149 (MQ=255)                                                                                                                                                  cggAGCGGCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCCGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAt  <  5:476506/147‑1 (MQ=255)                                                                                                                                                  cTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAt  <  8:292618/149‑1 (MQ=255)                                                                                                                                                  cTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAt  <  8:377089/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  ATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAT  >  NC_000913/858525‑858818

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