Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,664,599:1	+T	9.2%	coding (94/615 nt)	dmsD →	twin‑argninine leader‑binding protein for DmsA and TorA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,664,599	1	.	T	9.2%	257.6 / 21.6	65	coding (94/615 nt)	dmsD	twin‑argninine leader‑binding protein for DmsA and TorA
Reads supporting (aligned to +/- strand):  ref base . (33/26);  new base T (3/3);  total (36/29)
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.88e-01

GTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGC  >  NC_000913/1664459‑1664747                                                                                                                                              |                                                                                                                                                     gTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCActctct                                                                                                                                                  >  1:1524136/1‑140 (MQ=255)          ttGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCgg                                                                                                                                          <  2:1643171/144‑1 (MQ=255)          ttGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCgg                                                                                                                                          >  1:1643171/1‑144 (MQ=255)                cGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGc                                                                                                                                            >  1:685642/1‑136 (MQ=255)                cGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGc                                                                                                                                            <  2:685642/136‑1 (MQ=255)                         ccGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCACTCTCTgga                                                                                                                                               <  3:184090/124‑2 (MQ=255)                         ccGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCACTCTCTgga                                                                                                                                               >  4:184090/1‑123 (MQ=255)                            cacaTTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCACTCTCTgga                                                                                                                                               <  3:352476/121‑2 (MQ=255)                            cacaTTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAG‑GTGCGGCTTCCGCACTCTCTgga                                                                                                                                               >  4:352476/1‑120 (MQ=255)                                            aGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTAc                                                                                                   >  4:548376/1‑149 (MQ=255)                                              gAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCa                                                                                                 >  4:20780/1‑149 (MQ=255)                                                    cccATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGt                                                                                                            <  1:152908/132‑1 (MQ=255)                                                    cccATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGt                                                                                                            >  2:152908/1‑132 (MQ=255)                                                     ccATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGt                                                                                          >  3:1554352/1‑149 (MQ=255)                                                         ttttCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCAtt                                                                                      <  3:548376/149‑1 (MQ=255)                                                           ttCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAg                                                                                    <  4:1554352/149‑1 (MQ=255)                                                                    aGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCAt                                                                                       <  1:1089503/137‑1 (MQ=255)                                                                    aGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCAt                                                                                       >  2:1089503/1‑137 (MQ=255)                                                                     gATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGt                                                                          <  3:20780/149‑1 (MQ=255)                                                                               cTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCtttt                                                                >  3:1507179/1‑149 (MQ=255)                                                                                              tgtTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGc                                                                                            <  2:1387855/106‑1 (MQ=255)                                                                                              tgtTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGc                                                                                            >  1:1387855/1‑106 (MQ=255)                                                                                                ttGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGAc                                                                                                                            >  3:19363/1‑72 (MQ=255)                                                                                                ttGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGAc                                                                                                                            <  4:19363/72‑1 (MQ=255)                                                                                                  gggCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAg                                                                                                             >  2:1540282/1‑85 (MQ=255)                                                                                                  gggCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAg                                                                                                             <  1:1540282/85‑1 (MQ=255)                                                                                                    gcgcgcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCa                                                                                                                                 <  3:175138/63‑1 (MQ=255)                                                                                                    gcgcgcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCa                                                                                                                                 >  4:175138/1‑63 (MQ=255)                                                                                                      gcgcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCa                                         >  3:1894094/1‑149 (MQ=255)                                                                                                        gcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACg                                       >  1:1549931/1‑149 (MQ=255)                                                                                                        gcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACg                                       >  1:76170/1‑149 (MQ=255)                                                                                                           gTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGccc                                    >  2:665110/1‑149 (MQ=255)                                                                                                           gTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGccc                                    >  4:252659/1‑149 (MQ=255)                                                                                                           gTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGccc                                    >  4:564414/1‑149 (MQ=255)                                                                                                                 ttaCGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTg                                                                                                           <  1:1617916/72‑1 (MQ=255)                                                                                                                 ttaCGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTg                                                                                                           >  2:1617916/1‑72 (MQ=255)                                                                                                                  taCGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTg                             <  4:1507179/149‑1 (MQ=255)                                                                                                                      aaccAGAG‑GTGCGGCTTCCGCACTCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGa                                                  <  2:955469/122‑1 (MQ=38)                                                                                                                      aaccAGAG‑GTGCGGCTTCCGCACTCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGa                                                  >  1:955469/3‑124 (MQ=38)                                                                                                                       tCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGc                        <  3:564414/149‑1 (MQ=255)                                                                                                                       tCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGc                        >  1:1018841/1‑149 (MQ=255)                                                                                                                       tCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGc                        >  1:1517468/1‑149 (MQ=255)                                                                                                                       tCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGc                        <  1:665110/149‑1 (MQ=255)                                                                                                                         cAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGa                                                                        <  3:942896/99‑1 (MQ=255)                                                                                                                         cAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGa                                                                        >  4:942896/1‑99 (MQ=255)                                                                                                                          agagagTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGa                                                            <  1:530828/110‑1 (MQ=255)                                                                                                                          agagagTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGa                                                            >  2:530828/1‑110 (MQ=255)                                                                                                                             gagTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgt                  <  3:252659/149‑1 (MQ=255)                                                                                                                              agTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgtt                 >  3:1487747/1‑149 (MQ=255)                                                                                                                              agTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgtt                 <  4:1894094/149‑1 (MQ=255)                                                                                                                              agTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgtt                 <  2:1549931/149‑1 (MQ=255)                                                                                                                               gTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgttt                <  2:1517468/149‑1 (MQ=255)                                                                                                                                tGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGtttgtttg               >  1:194566/1‑149 (MQ=255)                                                                                                                                 gCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTttgtttgt              <  2:76170/149‑1 (MQ=255)                                                                                                                                  cGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTc             <  2:194566/149‑1 (MQ=255)                                                                                                                                  cGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTc             >  3:1579536/1‑149 (MQ=255)                                                                                                                                   ggAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCtttt                                                                >  2:812168/1‑97 (MQ=255)                                                                                                                                   ggAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCtgtt                                                                <  1:812168/97‑1 (MQ=255)                                                                                                                                    gAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgt                  <  4:1062035/142‑1 (MQ=255)                                                                                                                                    gAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgt                  >  3:1062035/1‑142 (MQ=255)                                                                                                                                       gCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGccc        <  2:1018841/149‑1 (MQ=255)                                                                                                                                         cGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCCGCGTTTGTTTGTCGGCCGGt      >  4:255833/1‑149 (MQ=255)                                                                                                                                          gCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTTGGCCCGTg     >  2:1532774/1‑149 (MQ=255)                                                                                                                                          gCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCTGCCCGTg     >  1:640012/1‑149 (MQ=255)                                                                                                                                          gCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTg     >  3:1644151/1‑149 (MQ=255)                                                                                                                                          gCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTg     >  3:1527900/1‑149 (MQ=255)                                                                                                                                             c‑cTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGAAACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGGCGGCCAGTAGGc  >  4:1464132/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                     GTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCAC‑CTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGC  >  NC_000913/1664459‑1664747

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