Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,878,894	T→A	5.2%	D409V (GAC→GTC)	gyrB ←	DNA gyrase, subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,878,894	0	T	A	5.2%	319.3 / 5.0	94	D409V (GAC→GTC)	gyrB	DNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base T (48/41);  new base A (3/2);  total (51/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCA  >  NC_000913/3878746‑3879038                                                                                                                                                     |                                                                                                                                                 gAAGCGCGCTTTCTCGACGTCGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGt                                                                                                                                                  >  2:70542/1‑149 (MQ=255)  aaGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGt                                                                                                                                                  <  4:824631/148‑1 (MQ=255)   aGCGCGCTTTCTCTACGTTGAAGATTTTCCCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCt                                                                                                                                                <  2:840494/149‑1 (MQ=255)   aGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCt                                                                                                                                                <  2:400847/149‑1 (MQ=255)   aGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCt                                                                                                                                                <  2:540983/149‑1 (MQ=255)     cgcgcTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGc                                                                                                                                              <  1:666369/149‑1 (MQ=255)       cgcTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGcc                                                                                                                                             >  3:441238/1‑148 (MQ=255)         cTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGACTGCCAgg                                                                                                                                          <  2:607294/149‑2 (MQ=255)          tttCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGtt                                                                                                                                         <  1:890428/149‑1 (MQ=255)          tttCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGACTGCCAgga                                                                                                                                         <  1:555686/149‑3 (MQ=255)           ttCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAgga                                                                                                                                         <  4:543858/148‑3 (MQ=255)           ttCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGttt                                                                                                                                        <  1:792442/149‑1 (MQ=255)            tctcGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGttt                                                                                                                                        <  3:488372/148‑1 (MQ=255)             ctcGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTg                                                                                                                                       <  3:204434/148‑1 (MQ=255)                    ttGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCa                                                                                                                                >  3:505247/1‑148 (MQ=255)                        ggATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGcc                                                                                                                            >  3:453561/1‑148 (MQ=255)                          aTTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCg                                                                                                                          <  3:819665/148‑1 (MQ=255)                            tttACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCt                                                                                                                        >  3:91078/1‑148 (MQ=255)                             ttACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTaa                                                                                                                      <  1:55349/149‑1 (MQ=255)                                  cTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCgt                                                                                                                 >  1:463320/1‑148 (MQ=255)                                  cTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCg                                                                                                                  <  4:269327/148‑1 (MQ=255)                                  cTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCg                                                                                                                  >  3:572864/1‑148 (MQ=255)                                  cTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCg                                                                                                                  <  3:539247/148‑1 (MQ=255)                                   ttCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGa                                                                                                                 >  3:811449/1‑148 (MQ=255)                                   ttCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAg                                                                                                                >  2:299787/1‑149 (MQ=255)                                       gCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAgcgc                                                                                                             <  4:970146/148‑1 (MQ=255)                                         ggCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGc                                                                                                                             <  1:1060114/130‑1 (MQ=255)                                         ggCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGc                                                                                                                             >  2:1060114/1‑130 (MQ=255)                                           cAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCt                                                                                                         >  3:76804/1‑148 (MQ=255)                                           cAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCAGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCtt                                                                                                        >  1:620470/1‑149 (MQ=255)                                              aaTCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAc                                                                                                     <  1:1071183/149‑1 (MQ=255)                                               aTCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAcg                                                                                                    <  1:70542/149‑1 (MQ=255)                                                   ccTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAcggcg                                                                                                 <  4:255887/148‑1 (MQ=255)                                                   ccTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAcggcg                                                                                                 <  3:591137/148‑1 (MQ=255)                                                    cTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGCGCGCACCTTTACGGCggg                                                                                               >  2:981160/1‑149 (MQ=255)                                                     tGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCggg                                                                                               >  4:487305/1‑148 (MQ=255)                                                       gTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCa                                                                                            >  1:583293/1‑149 (MQ=255)                                                        ttCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCAt                                                                                           >  2:198317/1‑149 (MQ=255)                                                          cTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCAttt                                                                                         <  1:92912/149‑1 (MQ=255)                                                           ttGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTc                                                                                        >  2:291782/1‑149 (MQ=255)                                                           ttGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTc                                                                                        >  1:219543/1‑149 (MQ=255)                                                                gTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCAcgc                                                                                    >  3:226041/1‑148 (MQ=255)                                                                 ttACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAcggc                                                                                                  <  3:305996/133‑1 (MQ=255)                                                                 ttACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAcggc                                                                                                  >  4:305996/1‑133 (MQ=255)                                                                 ttACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCAcgcgc                                                                                  <  1:198317/149‑1 (MQ=255)                                                                    cGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCAcgcgcgc                                                                                <  4:811449/148‑1 (MQ=255)                                                                      ccccTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGAc                                                                             <  2:662039/149‑1 (MQ=255)                                                                        ccTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGAcgc                                                                           >  1:379828/1‑149 (MQ=255)                                                                        ccTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGAcg                                                                            >  3:144463/1‑148 (MQ=255)                                                                          tGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGAcgcg                                                                          >  4:578296/1‑148 (MQ=255)                                                                          tGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGAcgcg                                                                          >  4:887233/1‑148 (MQ=255)                                                                          tGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACgcgc                                                                         >  1:965001/1‑149 (MQ=255)                                                                           gCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGcgcc                                                                        >  2:895343/1‑149 (MQ=255)                                                                            cTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGcgcc                                                                        >  3:792990/1‑148 (MQ=255)                                                                             ttCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAcggcgg                                                                                                <  2:954198/123‑1 (MQ=255)                                                                             ttCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAcggcgg                                                                                                >  1:954198/1‑123 (MQ=255)                                                                                  agagCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTca                                                                  <  4:111981/148‑1 (MQ=255)                                                                                      ccgcccgcGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTcacggg                                                              >  4:993119/1‑148 (MQ=255)                                                                                       cgcccgcGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTcacgggc                                                             >  3:958508/1‑148 (MQ=255)                                                                                        gcccgcGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCt                                                                                                         >  2:742210/1‑103 (MQ=255)                                                                                        gcccgcGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCt                                                                                                         <  1:742210/103‑1 (MQ=255)                                                                                          ccgcGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTcacgggcacg                                                          <  4:76804/148‑1 (MQ=255)                                                                                                     ccTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACTGGCACGGGCAGCATCGa                                               >  4:554593/1‑148 (MQ=255)                                                                                                     ccTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGa                                               >  4:573833/1‑148 (MQ=255)                                                                                                     ccTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGAt                                              >  1:1025919/1‑149 (MQ=255)                                                                                                     ccTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGAt                                              <  2:379828/149‑1 (MQ=255)                                                                                                      cTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATa                                             >  2:41044/1‑149 (MQ=255)                                                                                                       ttCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGGAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGGCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCATCATCGATa                                             >  3:73843/1‑148 (MQ=255)                                                                                                           accaGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAAtttt                                        <  2:219543/149‑1 (MQ=255)                                                                                                           accaGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAAttt                                         <  4:453561/148‑1 (MQ=255)                                                                                                              cggTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGc                                      >  4:726305/2‑148 (MQ=255)                                                                                                              aGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGc                                      <  3:887233/148‑1 (MQ=255)                                                                                                                gTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCaa                                   >  1:148690/1‑149 (MQ=255)                                                                                                                   cAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCa                                <  1:749807/149‑1 (MQ=255)                                                                                                                     gTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAAc                                  >  1:717626/1‑145 (MQ=255)                                                                                                                     gTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAAc                                  <  2:717626/145‑1 (MQ=255)                                                                                                                     gTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACg                              >  2:808126/1‑149 (MQ=255)                                                                                                                      ttCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGa                             <  1:41044/149‑1 (MQ=255)                                                                                                                      ttCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGa                             <  2:936405/149‑1 (MQ=255)                                                                                                                        cGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGAt                            <  3:573833/148‑1 (MQ=255)                                                                                                                             aGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGGTAATTTTGCCAACCACGATTTTcg                       >  4:1022608/1‑148 (MQ=255)                                                                                                                                   ggATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCTACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGt                >  2:434646/1‑149 (MQ=255)                                                                                                                                         cgTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGt           >  3:209471/1‑148 (MQ=255)                                                                                                                                          gTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGttt         <  1:469217/149‑1 (MQ=255)                                                                                                                                            tCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTcg                       <  4:130361/133‑1 (MQ=255)                                                                                                                                            tCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTcg                       >  3:130361/1‑133 (MQ=255)                                                                                                                                            tCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGtttt        <  4:226041/148‑1 (MQ=255)                                                                                                                                            aacTGGCAGACTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTc       >  2:388319/3‑149 (MQ=255)                                                                                                                                            aacTGGCAGACTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTc       >  2:918987/3‑149 (MQ=255)                                                                                                                                             ccTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTc       >  4:328345/1‑148 (MQ=255)                                                                                                                                               tGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcag    >  1:728224/1‑149 (MQ=255)                                                                                                                                               tGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcag    >  1:607332/1‑149 (MQ=255)                                                                                                                                               tGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcag    >  2:1000629/1‑149 (MQ=255)                                                                                                                                               tGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCca     >  3:825900/1‑148 (MQ=255)                                                                                                                                                 gCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcagca  <  1:895343/149‑1 (MQ=255)                                                                                                                                                 gCAGACTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTTATTTCACGCGCGCGACGCGCCGCTTCACTGGCACGGGCCGCATCGATAATTTTGTCAACCACGATTTTCGCTTCGGTTGGGGTTTCcagca  >  1:530392/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCA  >  NC_000913/3878746‑3879038

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