Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,379,266	(C)7→8	coding (732/1128 nt)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,259	1	.	C	86.9%	185.4 / 24.6	61	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand):  ref base . (1/7);  new base C (25/28);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

CAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTT  >  NC_000913/3379130‑3379403                                                                                                                                    |                                                                                                                                                 cAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAAc                                                                                                                               >  2:733909/1‑150 (MQ=255) cAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAAc                                                                                                                               >  1:620095/1‑150 (MQ=255)   atatAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCTTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCAt                                                                                                                            <  2:86468/151‑1 (MQ=255)   atatAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCAt                                                                                                                            >  2:838884/1‑151 (MQ=255)   atatAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCAt                                                                                                                            >  2:601548/1‑151 (MQ=255)      tattCATGTTGCCTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGCCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGt                                                                                                                         <  2:764118/148‑1 (MQ=255)      taGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGCCCGCCAGCGTCTGGTTCTCGACGCCCTTTTTTGCTAATGCCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGt                                                                                                                         <  1:695577/151‑1 (MQ=255)      taGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCCTAGt                                                                                                                         >  2:503616/1‑151 (MQ=255)      taGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGt                                                                                                                         >  2:608523/1‑151 (MQ=255)      taGTCAGGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGCCCGCCAGCGTCTGGTTCTCGACGCCCATTTTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGt                                                                                                                         <  2:53196/151‑1 (MQ=255)        gTCATGCTGCCTGCGGGCGTCGCCCCACAGCGTAGTAAAAGAGCCCGCAAGCGTCTGGTTCTCGACGCCCATTTTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGtt                                                                                                                        <  1:1002002/150‑1 (MQ=255)        gTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGttt                                                                                                                       >  2:960426/1‑151 (MQ=255)         tctggCTGACTGCGGGCGTCGACGCACAGGGTAGTAAAGGAGGCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGtt                                                                                                                        <  2:624223/145‑1 (MQ=255)         tCATGCTGACTGCGGGCTTCGACGCACAGGGTGGTAAAAGAGACCGCCAGCGTTTGGTTCTCGCCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGtt                                                                                                                        <  2:1044692/149‑1 (MQ=255)         tCAGGCTGACTGCGGGCGCCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGGCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGtt                                                                                                                        <  2:139725/149‑1 (MQ=255)             gCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCCCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGCTTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCa                                                                                                                  >  1:273193/1‑151 (MQ=255)              gcaacaatgggCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCAt                                                                                                                 >  2:425421/9‑151 (MQ=255)                 aCTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGCCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                                                                              <  1:521546/151‑1 (MQ=255)                      gggCTTCGACCCCCAGTTTAGTAAAAGAGACCCCCAGCGTCTGGTTCTCGACGCCCTTTTTTGCTAATGGCCGATGGTTGCTTTCTACCGTCGGTGAATTTTCTCTTT‑TCCCCCCCCGCCAGCGCCAACCATTGTTTATCCATTTGCGcc                                                                                                          <  2:371416/150‑1 (MQ=255)                       ggCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                                                                        >  1:483021/1‑151 (MQ=255)                                     gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGaa                                                                                          >  2:509652/1‑151 (MQ=255)                                     gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGGGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCCTCGTGaa                                                                                          >  2:126450/1‑151 (MQ=255)                                        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCCGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                                                                       >  2:768062/1‑151 (MQ=255)                                        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                                                                       >  1:511970/1‑151 (MQ=255)                                         agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                                                                      >  2:8588/1‑151 (MQ=255)                                          gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                                                                      >  1:171124/1‑150 (MQ=255)                                          gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                                                                      >  2:807786/1‑150 (MQ=255)                                          gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                                                                      >  2:589631/1‑150 (MQ=255)                                                   cccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTGGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                                                                              <  1:874401/148‑1 (MQ=255)                                                          gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                                                                     >  2:995897/1‑151 (MQ=255)                                                           catCTGGTTCTCGACGCCCTTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGTTGAATTTTCTCGTTCC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                                                                       <  1:960569/146‑1 (MQ=255)                                                                ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag                                                               >  1:817402/1‑151 (MQ=255)                                                                 gTTCTCGACGCCCTTTTTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag                                                               <  2:770904/150‑1 (MQ=255)                                                                    ctcGACGCCCATTTTTGCTAATGGCCGATGTTTGATTTCTAACGTCGGTAAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag                                                           <  1:273700/151‑1 (MQ=255)                                                                      cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTAAATTTTCTCGTTCC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg                                                         <  1:312727/151‑1 (MQ=255)                                                                          gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg                                                     >  1:574635/1‑150 (MQ=255)                                                                           cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc                                                    <  1:640413/151‑1 (MQ=255)                                                                             cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAg                                                  >  2:564081/1‑151 (MQ=255)                                                                              ttttttgCAAATGGCCGATGGTTGATTTCTAACGTCGTTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGa                                                 <  2:238236/147‑1 (MQ=255)                                                                              tttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGa                                                 <  1:291959/150‑1 (MQ=255)                                                                                     caaaTGGCCGAGGTTTGATTTCTAACGTGGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAAt                                          <  2:419510/149‑1 (MQ=255)                                                                                         tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCaa                                       >  2:998354/1‑150 (MQ=255)                                                                                                      aTTTCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATACATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACg                         >  1:44205/1‑151 (MQ=255)                                                                                                      aTTTCTAACGTCGGGGAATTTTCGCGTTCC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACg                         <  1:575228/151‑1 (MQ=255)                                                                                                       tttCTAACGTCGGGGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGt                        >  1:211543/1‑151 (MQ=255)                                                                                                         tCTAACGTCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTc                      <  1:760442/151‑1 (MQ=255)                                                                                                             acttCGGTGATTTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTa                  <  2:48232/148‑1 (MQ=255)                                                                                                              cgtcgtGGATTTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTAc                 <  1:712285/151‑1 (MQ=255)                                                                                                                tCGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGCTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACAt               >  1:875324/1‑151 (MQ=255)                                                                                                                 cGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATc              <  1:425266/151‑1 (MQ=255)                                                                                                                 cGGTGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATc              <  1:141496/151‑1 (MQ=255)                                                                                                                    taaaTTTTCTCGTTCC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATcac            <  2:620079/148‑1 (MQ=255)                                                                                                                    tGAATTTTCTCGTTCC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATcac            <  1:601301/150‑1 (MQ=255)                                                                                                                    tGAATTTTCTCGTTCC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATcac            <  2:581586/150‑1 (MQ=255)                                                                                                                    tGAATTTTCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATcac            <  2:906700/150‑1 (MQ=255)                                                                                                                        ttttCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGc       <  1:376219/151‑1 (MQ=255)                                                                                                                        ttttCTCGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGc       <  2:584097/151‑1 (MQ=255)                                                                                                                          tttccGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACTCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTg     <  2:214345/147‑1 (MQ=255)                                                                                                                          ttCTCGTTCC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTg     <  1:995541/151‑1 (MQ=255)                                                                                                                             tcGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGttt  <  2:234039/151‑1 (MQ=255)                                                                                                                             ccGTT‑TCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGttt  <  1:780239/150‑1 (MQ=255)                                                                                                                                    |                                                                                                                                                 CAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCACAATGCTGTTT  >  NC_000913/3379130‑3379403

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