Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,993,436	G→A	6.3%	E762K (GAA→AAA)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,993,436	0	G	A	6.3%	163.6 / 2.5	63	E762K (GAA→AAA)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (21/38);  new base A (2/2);  total (23/40)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

AACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATC  >  NC_000913/3993289‑3993580                                                                                                                                                    |                                                                                                                                                 aaCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcgacga                                                                                                                                                 >  2:718087/1‑149 (MQ=255)   cAAACTGCACGGCCTGTCAGTGCCGGTTGAAACCAATCACGTCAAATTACAGGAGGTGGTGGACGGCTATGCCAGCGAAGTGATCATCCAGTTATTATTCGAAGAAACGCAAGCCGAGAATGGCTTGAATCTCTACAGACGCGACGcaa                                                                                                                                               >  2:698346/1‑149 (MQ=255)   cAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGGGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaaa                                                                                                                                               <  1:501616/149‑1 (MQ=255)   cAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaaa                                                                                                                                               <  1:387599/149‑1 (MQ=255)    aaaCAGCACGGCCTGTAAGTGCAGGTTGAAACCAATCAAGACAAATTACCGGCGGTGATGGACGCCATTGCCAGCGAAGGGATCATCCAGATCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCGCGACGaacg                                                                                                                                              >  2:189769/1‑147 (MQ=255)       cTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCaa                                                                                                                                           >  1:463558/1‑149 (MQ=255)         gCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAAcc                                                                                                                                         >  2:939834/1‑149 (MQ=255)            cGGCCTGGCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCATTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCggg                                                                                                                                      <  2:808967/149‑1 (MQ=255)            cGGCATGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAAGGGCTTTAATATCTACATTCTCGACGAAAGCAACCggg                                                                                                                                      <  1:703780/149‑1 (MQ=255)                 tGTTAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGTCTTTGCCAGCGATGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAAGCGGGTTGAg                                                                                                                                 <  2:825326/149‑1 (MQ=255)                  gTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAgg                                                                                                                                <  1:727048/149‑1 (MQ=255)                      gttcAGGTTGAAACCAATCACGTCAAATTACCGGCGGGGGGGGACGGCTTTGCCAGCGACGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGTCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGtata                                                                                                                            <  1:617449/146‑1 (MQ=255)                      gTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACAAAAGCAACCGGGTTGAGGtata                                                                                                                            >  2:48170/1‑149 (MQ=255)                      gTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATAATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATTGCTTTAATATCTACATTCTCGACAAAAGCAACCGGGTTGAGGtata                                                                                                                            >  2:958375/1‑149 (MQ=255)                          aGGTTGATACCATTCACGTCAAATTTACGGCGTTGGTGGACGTCTTTGCGAGCGAAGGTATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATTTCTACTTTCTCGACGAAAGCAACCGGGTTGAGGTATATcac                                                                                                                        <  2:904741/149‑1 (MQ=255)                            gTTGAAACCAATCACGTCAAATTACCGGCGGGGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCACGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTAGATcacca                                                                                                                      <  2:195395/149‑1 (MQ=255)                            gTTGAAACCAATCACGTCAAATAACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACTAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATcacca                                                                                                                      <  2:197306/149‑1 (MQ=255)                                 aaCCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGCAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATTTCACCACTGCg                                                                                                                 >  2:898652/1‑149 (MQ=255)                                   ccAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTGTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGaa                                                                                                               >  2:186749/1‑149 (MQ=255)                                       tCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGgca                                                                                                           >  1:59772/1‑149 (MQ=255)                                        cACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGgcag                                                                                                          <  1:939834/149‑1 (MQ=255)                                        cACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGgcag                                                                                                          <  2:697273/149‑1 (MQ=255)                                              aaaTTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAaga                                                                                                    <  1:227987/149‑1 (MQ=255)                                                     cggcggTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTg                                                                                             >  2:409461/1‑149 (MQ=255)                                                       gcggTGGTGGACGGCTTTGCCAGCGAAGTGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGt                                                                                           <  1:958655/149‑1 (MQ=255)                                                            ggtggACGGCGTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGACACGCAAGCCGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTg                                                                                      <  1:750729/149‑1 (MQ=255)                                                                gACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGt                                                                                  <  1:9307/149‑1 (MQ=255)                                                                gACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGt                                                                                  <  1:17716/149‑1 (MQ=255)                                                                   ggCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAg                                                                               >  2:456624/1‑149 (MQ=255)                                                                    gCTTTGCCAGCGAAGTGATCATCTAGTTCTTTGTCGAAGAAACGCAAGACGAGAATGGGTTTAATATCTAGTTTCTCGACGAATGCAACGGGGTTGAGGTATATCACCACTGCGAAGGCAGAAAAGAGGAGCTGGTACGTGACGTCAGt                                                                              <  1:567790/149‑1 (MQ=255)                                                                     cTTTGCCAGCGAAGGGATCATCCATTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGGTGAGGTAGATCACCACTGCGAAGGCAGCAAAGTGGAGCTGGTCCGTGACGTCAGTc                                                                             >  2:859192/1‑149 (MQ=255)                                                                       ttGCCAGCGAAGGTATCAACCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCATCAAAGAGGAGCTGGTACGTGACGTCAGTCGc                                                                           <  1:905392/149‑1 (MQ=255)                                                                       ttGCCAGCGAAGGGATCATCCACTACTCTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATCTCTACATTCACGACGAAAGCCACCGGGTTGAGGTATGGCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGc                                                                           <  2:962210/149‑1 (MQ=255)                                                                         gCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCtt                                                                         >  1:725900/1‑149 (MQ=255)                                                                         gCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCtt                                                                         <  1:882550/149‑1 (MQ=255)                                                                          ccAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTc                                                                        <  1:967437/149‑1 (MQ=255)                                                                          ccAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTc                                                                        <  1:6066/149‑1 (MQ=255)                                                                          ccAGCGAAGGGATCATCCAGTTATTATTCGAAGAAACGCAAGAGGAGAATCGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTc                                                                        <  1:970625/149‑1 (MQ=255)                                                                             gCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTAc                                                                     >  2:49880/1‑149 (MQ=255)                                                                               gAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGGTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTc                                                                   <  1:641938/149‑1 (MQ=255)                                                                                   ggATCATCCAGTTCTTATTCCAAGAAACGCAAGACGAGAAAGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCa                                                               <  2:442295/149‑1 (MQ=255)                                                                                    gATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGCCGAAAGCAACCGGGGGGAGGTATATCACCACTCCGAAGGCAGCAAAGAGGATCTGGAACGAGACGTCGCTCTCTGCTACTAGTCAt                                                              >  1:694565/1‑149 (MQ=255)                                                                                     atcatcCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATc                                                             >  1:187876/1‑149 (MQ=255)                                                                                     atcatcCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATc                                                             <  1:2498/149‑1 (MQ=255)                                                                                        atcCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCa                                                          >  2:30015/1‑149 (MQ=255)                                                                                         tcCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAt                                                         >  1:130082/1‑149 (MQ=255)                                                                                          ccAGTTCTTTTTCGATGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACAGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGTGCTGGTACGTGACGTCAGTCGCTGGTACTCGTCATCGCATg                                                        >  2:684912/1‑149 (MQ=255)                                                                                                cTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACTTCAGTCGCTTCTACTCGTCATCGCATGACCGtt                                                  >  1:957570/1‑149 (MQ=255)                                                                                                  ttttCGAAGAAACGCAAGACTAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGAAGCAAAGAGGAGCTGGTTCGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACAGtttt                                                <  2:502546/149‑1 (MQ=255)                                                                                                    ttCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTAc                                              >  2:116034/1‑149 (MQ=255)                                                                                                      cGAAGAAACGCAAGACGAGACTCGCTTTAATATCAAGCTTCTCGACGAACGCAACCGGGGTGAGGTATATCACCCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCt                                            <  2:678493/149‑1 (MQ=255)                                                                                                      cGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGGTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCt                                            <  2:463558/149‑1 (MQ=255)                                                                                                        aagaaACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTAc                                          >  1:55511/1‑149 (MQ=255)                                                                                                                 aaGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGc                                 <  1:520347/149‑1 (MQ=255)                                                                                                                   taCGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGGTTTACCTACGGCTCAAGCtt                               <  1:208628/148‑1 (MQ=255)                                                                                                                    aCGAGAATGGCTTTAATATCTACTTTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtc                              <  1:519340/149‑1 (MQ=255)                                                                                                                       cgacggccTTTAATATCTACATCCTCGACGAAAGCCACCGGGTTGAGGTATATCACCCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtcatc                           <  1:708273/142‑1 (MQ=255)                                                                                                                           tGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACt                       <  2:113820/149‑1 (MQ=255)                                                                                                                               tttAATATCTACTTTCGCGACGAAAGCCTCCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTTGTAAGTGACGTCATTCGCTTCTACTCGTCATCGCATGACCGTTTTACGTACGGCTCAAGCTTCATCAACTTCac                   <  1:189769/149‑2 (MQ=255)                                                                                                                                        tACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCa          >  1:810077/1‑149 (MQ=255)                                                                                                                                          cATTCTCGACAAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCACGCTTCATCAACTTCAACCTGCCGCAGt        <  1:48170/149‑1 (MQ=255)                                                                                                                                          aaTTCTCGACAAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGt        <  1:958375/148‑1 (MQ=255)                                                                                                                                               tcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTAt   <  2:342782/149‑1 (MQ=255)                                                                                                                                               tcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTAt   <  2:158453/149‑1 (MQ=255)                                                                                                                                                cgacgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATc  <  1:30015/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 AACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATC  >  NC_000913/3993289‑3993580

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