Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,810,065	T→G	8.2%	intergenic (+17/‑146)	yniC → / → ydjM	2‑deoxyglucose‑6‑P phosphatase/inner membrane protein regulated by LexA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,810,065	0	T	G	8.2%	202.4 / 10.5	61	intergenic (+17/‑146)	yniC/ydjM	2‑deoxyglucose‑6‑P phosphatase/inner membrane protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base G (2/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATC  >  NC_000913/1809917‑1810210                                                                                                                                                     |                                                                                                                                                  gCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAt                                                                                                                                                   <  2:605100/149‑1 (MQ=255)  ccTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAt                                                                                                                                                   >  4:133901/1‑148 (MQ=255)    tctAAAGCAGCCCGCATGCGTTCCATCATGGTTCCTGCGCCAGAAGCGCATAATGATCCACGTTTTGTATTAGCAGACGTCATACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtctc                                                                                                                                                >  2:701412/1‑149 (MQ=255)     ctAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtctc                                                                                                                                                >  3:53703/1‑148 (MQ=255)       aaaGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGc                                                                                                                                              >  3:148089/1‑148 (MQ=255)        aaGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCt                                                                                                                                             >  4:446996/1‑148 (MQ=255)        aaGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTg                                                                                                                                            >  1:362655/1‑149 (MQ=255)           cagcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGccc                                                                                                                                         <  2:156286/149‑1 (MQ=255)            agcCCGGATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGcccc                                                                                                                                        <  2:616272/149‑1 (MQ=255)             gcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGcccc                                                                                                                                        <  4:517889/148‑1 (MQ=255)                       gTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                                                                                             <  2:548100/149‑1 (MQ=255)                        ttCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                                                                                             >  4:46344/1‑148 (MQ=255)                            aTCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGa                                                                                                                         >  4:498682/1‑148 (MQ=255)                            aTCGTCGTTCCGGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaa                                                                                                                        >  1:797660/1‑149 (MQ=255)                             tcgtcgTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaaa                                                                                                                       >  2:892019/1‑149 (MQ=255)                                 cgtTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGc                                                                                                                   <  1:372362/149‑1 (MQ=255)                                  gtTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGc                                                                                                                   <  4:436945/148‑1 (MQ=255)                                   ttCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGccc                                                                                                                                         <  3:517269/125‑1 (MQ=255)                                    tCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAt                                                                                                                 <  4:54452/148‑1 (MQ=255)                                     tttGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAtt                                                                                                                >  4:673533/3‑148 (MQ=255)                                      cTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTc                                                                                                               >  3:966868/1‑148 (MQ=255)                                       tGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCa                                                                                                              <  4:466804/148‑1 (MQ=255)                                          gcCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTg                                                                                                          >  2:122345/1‑149 (MQ=255)                                                      aaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                                                                                               >  4:972439/1‑148 (MQ=255)                                                          gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                                                                           <  3:164901/148‑1 (MQ=255)                                                                 gTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                                                                   <  1:614528/149‑1 (MQ=255)                                                                       tATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                                                                             <  1:892019/149‑1 (MQ=255)                                                                            gCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTac                                                                        >  2:96902/1‑149 (MQ=255)                                                                                  gTCTAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTc                                                                  >  1:926642/1‑149 (MQ=255)                                                                                   tCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                                                                                            <  4:856527/106‑1 (MQ=255)                                                                                   tCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                                                                                            >  3:856527/1‑106 (MQ=255)                                                                                    cAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                                                                             <  4:682766/136‑1 (MQ=255)                                                                                    cAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                                                                             >  3:682766/1‑136 (MQ=255)                                                                                           tCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTcc                                                         <  2:362655/149‑1 (MQ=255)                                                                                           tCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTc                                                          <  4:148089/148‑1 (MQ=255)                                                                                             aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGt                                                       >  1:676235/1‑149 (MQ=255)                                                                                              tCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGt                                                                                             >  1:796352/1‑110 (MQ=255)                                                                                              tCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGt                                                                                             <  2:796352/110‑1 (MQ=255)                                                                                                 cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCa                                                    <  4:53703/148‑1 (MQ=255)                                                                                                       gAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGtat                                              >  3:423519/1‑148 (MQ=255)                                                                                                            cACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTATGAGAAAGAGTTCCACCGTCACGACTTCCGTACACTGTATaaaaa                                         >  4:20161/1‑148 (MQ=255)                                                                                                            cACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATaaaaa                                         >  3:317865/1‑148 (MQ=255)                                                                                                            cACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATaaaaa                                         >  4:768480/1‑148 (MQ=255)                                                                                                            cACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATaaaaa                                         >  4:849112/1‑148 (MQ=255)                                                                                                                gCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCt                                     <  3:46344/148‑1 (MQ=255)                                                                                                                gCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCt                                     <  4:423519/148‑1 (MQ=255)                                                                                                                  aaaaGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCtat                                   <  4:317865/148‑1 (MQ=255)                                                                                                                   aaaGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATAc                                 <  1:959795/149‑1 (MQ=255)                                                                                                                      gACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGt                              >  2:646284/1‑149 (MQ=255)                                                                                                                        ccTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTa                             <  3:972439/148‑1 (MQ=255)                                                                                                                           tctcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGt                          >  4:605798/1‑148 (MQ=255)                                                                                                                           tctcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTa                         >  2:687163/1‑149 (MQ=255)                                                                                                                             tcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTAt                        <  3:446996/148‑1 (MQ=255)                                                                                                                                gTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGa                     >  4:987086/1‑148 (MQ=255)                                                                                                                                 ttAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGAc                    >  3:141147/1‑148 (MQ=255)                                                                                                                                     tGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGtt                <  4:966868/148‑1 (MQ=255)                                                                                                                                     gaaCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGttt               >  1:705690/3‑149 (MQ=255)                                                                                                                                      gACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTa              <  1:174216/149‑1 (MQ=255)                                                                                                                                        ccAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGt            >  2:88052/1‑149 (MQ=255)                                                                                                                                           ggggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGa          >  4:270448/1‑148 (MQ=255)                                                                                                                                            gggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGt        >  1:898936/1‑149 (MQ=255)                                                                                                                                              gCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGtt       <  3:634872/148‑1 (MQ=255)                                                                                                                                                aGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGtttt     <  3:768480/148‑1 (MQ=255)                                                                                                                                                 gTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTAt   <  1:687163/149‑1 (MQ=255)                                                                                                                                                  tGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATc  <  1:701412/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GCCTCTAAAGCAGCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTATC  >  NC_000913/1809917‑1810210

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