Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	704,018	A→G	A25A (GCA→GCG)	nagE →	N‑acetylglucosamine specific PTS enzyme IIABC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	704,018	0	A	G	98.7%	266.5 / ‑4.8	78	A25A (GCA→GCG)	nagE	N‑acetylglucosamine specific PTS enzyme IIABC component
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (41/36);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TGCAGCAAAATAAAAATACGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGT  >  NC_000913/703883‑704153                                                                                                                                        |                                                                                                                                        tGCAGCAAAATAAAAATACGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGt                                                                                                                                     >  2:1116129/1‑140 (MQ=255)    agcaAAATAAAAATACGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTtgc                                                                                                                                  >  2:908438/1‑140 (MQ=255)    agcaAAATAAAAATACGGCTTGAAACGAGCCAAATAGGGGTCTCGCAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTtgc                                                                                                                                  >  2:782947/1‑140 (MQ=255)     gcaAAATAAAAATACGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTtgct                                                                                                                                 <  2:401898/140‑1 (MQ=255)        aaataaaaataCGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCg                                                                                                                              >  2:116859/1‑140 (MQ=255)         aataaaaataCGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGa                                                                                                                             >  2:1221034/1‑140 (MQ=255)             aaaataCGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCg                                                                                                                         <  2:430950/140‑1 (MQ=255)                  aCGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAg                                                                                                                    <  2:362371/140‑1 (MQ=255)                      cTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAg                                                                                                                >  1:173263/1‑140 (MQ=255)                        tGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGAt                                                                                                              <  1:1221034/140‑1 (MQ=255)                        tGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGAt                                                                                                              <  2:441833/140‑1 (MQ=255)                               aGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTa                                                                                                       <  2:366095/140‑1 (MQ=255)                                   aaaTAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGt                                                                                                   <  1:116859/140‑1 (MQ=255)                                    aaTAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGtt                                                                                                  <  2:525805/140‑1 (MQ=255)                                             ctcGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTACCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTAtt                                                                                         >  1:241080/1‑140 (MQ=255)                                              tcGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTg                                                                                        >  2:742444/1‑140 (MQ=255)                                              tcGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTg                                                                                        >  2:329853/1‑140 (MQ=255)                                              tcGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTg                                                                                        >  2:210041/1‑140 (MQ=255)                                              tcGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTg                                                                                        >  1:172341/1‑140 (MQ=255)                                              tcGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTg                                                                                        >  2:1067805/1‑140 (MQ=255)                                              tcGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTg                                                                                        <  2:173263/140‑1 (MQ=255)                                                   gggggAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAg                                                                                   >  2:1368530/1‑140 (MQ=255)                                                   gggggAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAg                                                                                   >  1:534944/1‑140 (MQ=255)                                                        aaTAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAggcggg                                                                              >  1:1107176/1‑140 (MQ=255)                                                        aaTAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAggcggg                                                                              >  1:731431/1‑140 (MQ=255)                                                        aaTAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAggcggg                                                                              >  2:1028830/1‑140 (MQ=255)                                                         aTAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAggcgggc                                                                             <  1:908438/140‑1 (MQ=255)                                                                  atatTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGAtt                                                                    >  1:580632/1‑140 (MQ=255)                                                                  atatTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGAtt                                                                    >  2:1405464/1‑140 (MQ=255)                                                                  atatTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGAtt                                                                    >  2:207615/1‑140 (MQ=255)                                                                  atatTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGAtt                                                                    >  2:1019849/1‑140 (MQ=255)                                                                  atatTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGAtt                                                                    >  2:1007282/1‑140 (MQ=255)                                                                   tatTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGAttt                                                                   >  1:999242/1‑140 (MQ=255)                                                                      tttAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTg                                                                >  1:599854/1‑140 (MQ=255)                                                                       ttAGGTTTTTTCCAGCGACTCGGTAGGGCGTTTCAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGa                                                               <  2:690808/140‑1 (MQ=255)                                                                        tAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGAt                                                              >  1:135177/1‑140 (MQ=255)                                                                           gTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAAc                                                           >  2:805521/1‑140 (MQ=255)                                                                           gTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAAc                                                           <  1:329853/140‑1 (MQ=255)                                                                            ttttttCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAAcc                                                          <  1:805521/140‑1 (MQ=255)                                                                            ttttttCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAAcc                                                          <  2:172341/140‑1 (MQ=255)                                                                            ttttttCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAAcc                                                          <  2:580632/140‑1 (MQ=255)                                                                            ttttttCCAGCGACTCGGTAGGGCGTTACAGCTCCCTACCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAAcc                                                          <  2:241080/140‑1 (MQ=255)                                                                             tttttCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCt                                                         <  1:1028830/140‑1 (MQ=255)                                                                             tttttCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCt                                                         <  1:1019849/140‑1 (MQ=255)                                                                               tttCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCg                                                       >  2:368650/1‑140 (MQ=255)                                                                                       gACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATc                                               <  1:1007282/140‑1 (MQ=255)                                                                                         cTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCtt                                             <  1:88817/140‑1 (MQ=255)                                                                                            ggTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGc                                          <  2:534944/140‑1 (MQ=255)                                                                                             gTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGcc                                         <  2:731431/140‑1 (MQ=255)                                                                                             gTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGcc                                         <  2:38117/140‑1 (MQ=255)                                                                                                gggCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATc                                      >  2:662117/1‑140 (MQ=255)                                                                                                gggCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATc                                      >  2:589041/1‑140 (MQ=255)                                                                                                   cGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGgt                                   >  2:231100/1‑140 (MQ=255)                                                                                                   cGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGgt                                   >  1:1221851/1‑140 (MQ=255)                                                                                                    gTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTGATCTTCGCCATCGgtg                                  <  2:135177/140‑1 (MQ=255)                                                                                                    gTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGgtg                                  >  1:1242973/1‑140 (MQ=255)                                                                                                    gTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGgtg                                  <  2:599854/140‑1 (MQ=255)                                                                                                      tACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGtgtg                                >  2:866796/1‑140 (MQ=255)                                                                                                      tACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGtgtg                                >  2:1077751/1‑140 (MQ=255)                                                                                                       aCAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTgg                               >  2:618750/1‑140 (MQ=255)                                                                                                                tATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCt                      <  1:742444/140‑1 (MQ=255)                                                                                                                tATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCt                      <  1:458158/140‑1 (MQ=255)                                                                                                                tATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCt                      <  1:231100/140‑1 (MQ=255)                                                                                                                tATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCt                      >  2:1036171/1‑140 (MQ=255)                                                                                                                tATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCt                      <  1:1167470/140‑1 (MQ=255)                                                                                                                tATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCt                      <  2:443001/140‑1 (MQ=255)                                                                                                                tATCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCt                      <  1:1077751/140‑1 (MQ=255)                                                                                                                  tCGCGGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTgg                    <  2:999242/140‑1 (MQ=255)                                                                                                                   cgcgGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGt                   <  1:207615/140‑1 (MQ=255)                                                                                                                   cgcgGTGCTGCCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGt                   <  1:1405464/140‑1 (MQ=255)                                                                                                                        tgctgcCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGaaa              >  1:1304623/1‑140 (MQ=255)                                                                                                                        tgctgcCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGaaa              >  1:101695/1‑140 (MQ=255)                                                                                                                        tgctgcCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCCGTGTGGCATCCAGCTGGTCGaaa              >  1:643553/1‑140 (MQ=255)                                                                                                                          ctgcCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGa            <  2:1172579/140‑1 (MQ=255)                                                                                                                          ctgcCGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGa            <  1:866796/140‑1 (MQ=255)                                                                                                                             ccGGTGGCGGCGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAg         <  2:643553/140‑1 (MQ=255)                                                                                                                                  ggcggcGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAAACAGCGCTg    >  1:343462/1‑140 (MQ=255)                                                                                                                                    cggcGCTGTTGCTGCGATTCGGTCAGCCAGATTTACTAAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGt  >  1:492580/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        TGCAGCAAAATAAAAATACGGCTTGAAACGAGCCAAATAGGGTTCTCGTAGGGGGAATAAGATGAATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGT  >  NC_000913/703883‑704153

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