Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,992,624	T→G	7.9%	L491R (CTG→CGG)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,992,624	0	T	G	7.9%	203.7 / 7.3	61	L491R (CTG→CGG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (24/32);  new base G (2/3);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCG  >  NC_000913/3992479‑3992772                                                                                                                                                  |                                                                                                                                                     gATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACtggt                                                                                                                                                   <  1:193403/149‑1 (MQ=255)  aTCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACtggtg                                                                                                                                                  <  2:603414/149‑1 (MQ=255)  aTCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACtggtg                                                                                                                                                  >  1:606910/1‑149 (MQ=255)      ctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGt                                                                                                                                              <  2:606910/149‑1 (MQ=255)       tcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTg                                                                                                                                             >  3:613225/1‑149 (MQ=255)        cGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTgg                                                                                                                                            <  4:329791/149‑1 (MQ=255)        cGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTgg                                                                                                                                            <  1:658568/149‑1 (MQ=255)            aCCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCAt                                                                                                                                        <  4:30555/149‑1 (MQ=255)                aaTCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGtt                                                                                                                                    >  2:550623/1‑149 (MQ=255)                    tGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAAc                                                                                                                                <  1:136494/149‑1 (MQ=255)                     gACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACg                                                                                                                               >  2:139317/1‑149 (MQ=255)                       ccTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGc                                                                                                                             >  4:292066/1‑149 (MQ=255)                        cTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGcc                                                                                                                            >  3:570150/1‑149 (MQ=255)                         tttatttaTGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGc                                                                                                                             >  3:50138/1‑147 (MQ=255)                         tttatttaTGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGc                                                                                                                             <  4:50138/147‑1 (MQ=255)                         tttatttaTGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCct                                                                                                                           >  1:266082/1‑149 (MQ=255)                         tttatttaTGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCct                                                                                                                           >  1:240077/1‑149 (MQ=255)                                atGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGAcc                                                                                                                    <  1:158711/149‑1 (MQ=255)                                   tgCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAGCGGCCTGCTGACCTcg                                                                                                                 >  4:296947/1‑149 (MQ=255)                                    gccgccGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcgc                                                                                                                <  4:38455/149‑1 (MQ=255)                                        ccgggccgggCTAACCGTTCAGGTTAGTATCTGTTTAACCGCGCGCCAAATATTGCGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAAGGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCAc                                                                                                            <  4:374168/149‑1 (MQ=255)                                              cgggcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGttt                                                                                                      <  1:627297/149‑1 (MQ=255)                                                ggcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGt                                                                                                    <  1:83523/149‑1 (MQ=255)                                                 gcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGta                                                                                                   <  1:10008/149‑1 (MQ=255)                                                        gTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                                                            <  4:570150/149‑1 (MQ=255)                                                        gTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                                                            <  2:540246/149‑1 (MQ=255)                                                        gTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCAGCAGCCGCTGGAATATAACCGTTACCTGAATAAACTCGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                                                            <  4:119670/149‑1 (MQ=255)                                                         ttCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAg                                                                                           <  2:389590/149‑1 (MQ=255)                                                             ggTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTaa                                                                                       <  4:629991/149‑1 (MQ=255)                                                                  gTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAAc                                                                                                                                <  1:650807/103‑1 (MQ=255)                                                                  gTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAAc                                                                                                                                >  2:650807/1‑103 (MQ=255)                                                                        gTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCg                                                                            <  1:25439/149‑1 (MQ=255)                                                                        gTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCg                                                                            >  2:90478/1‑149 (MQ=255)                                                                           taACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGAtt                                                                         <  2:127/149‑1 (MQ=255)                                                                            aaCCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGAttt                                                                        <  1:608513/149‑1 (MQ=255)                                                                            aaCCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGAttt                                                                        <  2:96128/149‑1 (MQ=255)                                                                                gcgggcCAAATATTGAGTCGATCATCAGCAATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCt                                                                    <  2:247133/149‑1 (MQ=255)                                                                                      cAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGttt                                                                                                      <  1:563490/109‑1 (MQ=255)                                                                                      cAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGttt                                                                                                      >  2:563490/1‑109 (MQ=255)                                                                                        aaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCa                                                            <  2:428474/149‑1 (MQ=255)                                                                                        aaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCa                                                            <  1:139317/149‑1 (MQ=255)                                                                                         atatTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAg                                                           <  3:321571/149‑1 (MQ=255)                                                                                           atTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGga                                                         <  2:266082/149‑1 (MQ=255)                                                                                               aGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATg                                                     >  2:656133/1‑149 (MQ=255)                                                                                                 tCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGt                                                   >  4:231077/1‑149 (MQ=255)                                                                                                 tCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGt                                                   >  2:384147/1‑149 (MQ=255)                                                                                                   gATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCg                                                 >  3:397770/1‑149 (MQ=255)                                                                                                    atcatcAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGc                                                <  3:300549/149‑1 (MQ=255)                                                                                                     tcatcaGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGta                                                                                                   <  3:531099/97‑1 (MQ=255)                                                                                                     tcatcaGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGta                                                                                                   >  4:531099/1‑97 (MQ=255)                                                                                                       atcagccatcagccGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGa                                             >  4:155027/1‑149 (MQ=255)                                                                                                        tcagccatcagccGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGAc                                            >  4:127591/1‑149 (MQ=255)                                                                                                         cagccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                                                            <  3:609382/100‑1 (MQ=255)                                                                                                         cagccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                                                            >  4:609382/1‑100 (MQ=255)                                                                                                           gccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTAGCGTGGGCATGGATTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTCGCCTAAGTTGCAGGAGCTGGTCGCCGACgtg                                         >  3:422038/1‑149 (MQ=255)                                                                                                               tcagccGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGc                                     >  3:450903/1‑149 (MQ=255)                                                                                                               tcagccGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGc                                     <  1:550623/149‑1 (MQ=255)                                                                                                                        ggAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTcc                            >  2:628883/1‑149 (MQ=255)                                                                                                                        ggAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTcc                            <  4:313087/149‑1 (MQ=255)                                                                                                                           atatAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGcc                                                                     <  2:305242/105‑1 (MQ=255)                                                                                                                           atatAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGcc                                                                     >  1:305242/1‑105 (MQ=255)                                                                                                                                aCCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTgcgc                    <  2:489290/149‑1 (MQ=255)                                                                                                                                    ttACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTAc                >  1:661399/1‑149 (MQ=255)                                                                                                                                        cTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGc            >  2:265669/1‑149 (MQ=255)                                                                                                                                                  tggtggCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgacaccg  >  4:13379/1‑149 (MQ=255)                                                                                                                                                  tggtggCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgacaccg  >  4:634421/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCG  >  NC_000913/3992479‑3992772

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