Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,504,799	0	A	G	42.4%	81.4 / 32.0	68	G361G (GGA→GGG)	insQ	IS609 transposase B
Reads supporting (aligned to +/- strand):  ref base A (10/28);  new base G (28/0);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 1.47e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAAC  >  NC_000913/1504660‑1504948                                                                                                                                             |                                                                                                                                                         cgTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACG‑CCGGTCTTGGCCGTGGGGAGATgg                                                                                                                                                   >  1:370583/1‑147 (MQ=255)  gTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACG‑CCCGTCTTGCCTTTGGGGGGAAGGTg                                                                                                                                                 >  2:248681/1‑148 (MQ=255)  gTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCGCC‑CCCGTCCTTCCTGTGGGGGGATgg                                                                                                                                                   >  2:683506/1‑146 (MQ=255)    aCTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCCTTCTTTCCCGTGgggaggagggggagtc                                                                                                                                            >  1:108339/1‑140 (MQ=255)    aCTGGGGTAATACAGCAAAAGAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTc                                                                                                                                            <  1:828923/151‑1 (MQ=255)      tgtgGTCATACAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTc                                                                                                                                            <  2:999149/149‑1 (MQ=255)             aTACAGAGAAAGAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGgccg‑c                                                                                                                                    <  1:899196/149‑1 (MQ=255)              tACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACG‑CCGGTCTTTCCTGTGGGGGGATGGGGCGGGCAGGCCG‑ccc                                                                                                                                  >  2:338413/1‑150 (MQ=255)               aCAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACG‑CCGTTGTTGCCGGTGGGGAGAGGGGGCAGTCAGGCCG‑ccc                                                                                                                                  >  1:339583/1‑149 (MQ=255)                cAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACG‑CCCGTCTTTCCCGTGGGGAGAGGGGGccggaagggccccc                                                                                                                                   >  2:399567/1‑135 (MQ=255)                cAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGAAG‑CCGGTCTTTCCCGTGGGGGGAAGGGGCAGTGAGGGCC‑CCCgct                                                                                                                               >  1:963727/1‑149 (MQ=255)                    gAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACG‑CCGGTCTTTCCCGTGGGGGGAAGGGGGAGGCAGGGCG‑CCCgct                                                                                                                               >  2:735568/1‑145 (MQ=255)                     aaaGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGGCG‑CCGTTCTTTTCTGTGGGGAAATTGGGCTGTCAGGCCG‑CCCGTTgga                                                                                                                            >  2:443823/1‑147 (MQ=255)                      aaGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGGTCTTTCCTGTGGAGAGATGGGGCAGTCAGGCCG‑cccccttaag                                                                                                                           >  1:1166539/1‑142 (MQ=255)                        gagaACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCCGTCTTTCCCGTGGGGAGAAGGGGGAGGCAGGGCG‑CCCGCTG‑GAGCA‑gg                                                                                                                     >  2:1084406/1‑151 (MQ=255)                           aaCCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGaa                                                                                                                   <  2:779961/150‑1 (MQ=255)                             ccGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGGTCTTTCCTGTGGGGAGAAGGTGCAGTCAAGCCG‑CCCGCTG‑AAGGA‑GGAAccc                                                                                                                >  1:59700/1‑151 (MQ=255)                                   gTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGaa                                                                                                          <  2:270111/151‑1 (MQ=255)                                     cacaAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACC‑CCCGTCTTTCCTGTGGGGGGGTGGGGGAGTCAGGCCGCCCCGTTG‑‑AGCA‑GGGACCCACCCaaa                                                                                                         >  2:870850/1‑150 (MQ=255)                                     cacaAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCGCC‑CCGGTCTTTCCCGGGGGGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCC‑GGAACCCCCCgca                                                                                                          >  1:1042808/1‑147 (MQ=255)                                       caAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCGCG‑CCCTGCCTTCCTGTGGGGGAAAGGGGCAGCCCAGCCG‑GCCACTG‑GAACA‑AGAACCCAcccca                                                                                                          >  2:897554/1‑144 (MQ=255)                                          aGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTc                                                                                                   <  1:425933/151‑1 (MQ=255)                                           gTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGTGCCAG‑CCGGTTTTTTCTCTGGTGGGATAGGGCGGTCAGGTCG‑CCCCCTG‑GAGAC‑AGAAACCACCGAAAAGAtt                                                                                                    >  1:615628/1‑149 (MQ=255)                                           gTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGCTCTTTCCTGTGGGGGGAAGGTGCCGTTAGGGCG‑CCCCTTG‑AAGCA‑GGGAACCCCCCAAAAGAtt                                                                                                    >  2:402794/1‑149 (MQ=255)                                             aaaTTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACG‑GCGGTTTTTGCTGTGGAGGGGTGGGGGAGTCAGGGCG‑CCCCTTG‑AAGCA‑GGAAACCCCCGAAAAGAtt                                                                                                    >  1:1121845/1‑147 (MQ=255)                                             aaaTTCAGATGCCAGGCATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCCTTCTCGCCTTTGTGGAGATGGGGGCGTCGGGCCC‑CCCGTTG‑AAGAA‑AGAAACcccccggatg                                                                                                       >  1:587199/1‑138 (MQ=255)                                              aaTTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGc                                                                                               <  2:335787/151‑1 (MQ=255)                                              aaTTCAGATGCCAGGTATGTGAATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTGAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GAAACCCACCGAAATGATTCAGGc                                                                                               <  2:468836/151‑1 (MQ=255)                                               aTTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACG‑CCGGTCTTTCCTGTGGGGGGGAGGGGGAGGCAGGGCG‑CCCGCTG‑GAGCA‑GGAAACCACCGCAAAGATTCAgggg                                                                                              >  2:624818/1‑149 (MQ=255)                                                  cAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCCTTCTTTCCTGGGGGGAGATGGGGGAGTCAGGCCG‑CCCGCTG‑GAGGA‑GGAACCCACCGAAATGAtttaaggg                                                                                              >  2:167649/1‑142 (MQ=255)                                                       gCCAGGTATGTGGTTATACAGCGAGCCCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCTCCGAAATGATTCAGGCGACAGCCt                                                                                       <  2:245630/150‑1 (MQ=255)                                                       gCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTg                                                                                      <  2:629958/151‑1 (MQ=255)                                                       gCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGGTCTTTCCTGTGGGGAGAGGGGGCAGTCAGGCCC‑CCCGGTT‑AAGCA‑GGGACCCACCGAAAAGATTTAGGGGACAGcc                                                                                        >  1:791469/1‑149 (MQ=255)                                                        ccAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCCCG‑CCCTTCTTGCCTGTGGGGAGGTGGGGGCGTCAGGGCC‑CCCGTTG‑GAGCA‑GGGACCCACCCAAAAGATTCAGGCGGCAGCCTg                                                                                      >  1:179093/1‑150 (MQ=255)                                                         cAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACG‑CCGGTCTTGCCTGTGGGGGGATGGGGCAGGCAGGCCG‑CCCCTTG‑AAGCA‑GGAACCCACCCAAATGATTCAGGGGGCAGCCtta                                                                                     >  1:600254/1‑148 (MQ=255)                                                           ggTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACg                                                                                  <  1:307402/151‑1 (MQ=255)                                                            gTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACG‑CCCTTCTTTCCTGTGGGGGGATGGGGCAGTCAGGCCG‑CCCGTTG‑GAGCA‑GGAACCCCCCCCAAAGATTTAAGCGACAGCCTGaa                                                                                    >  2:335792/1‑148 (MQ=255)                                                            gTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACG‑CCCTTCTTGCCTGTGGGGAGGTGGTGCAGGCAGGCCG‑CCCCCTG‑GAGCA‑GGAACCCACCCAAAAGATTCAGGCGACCGCCTgga                                                                                    >  2:712218/1‑146 (MQ=255)                                                            gTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGGTCTTGCCTGTGGGGAGATGGGGGAGTCAGGCCG‑CCCCTTG‑AAGCA‑GGAACCCACCGAAAAGATTCAGGCGACAGCCCGaa                                                                                    >  1:1040020/1‑148 (MQ=255)                                                            gTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACC‑CCCTTCTTTCCTGTGGGGGGGTGGGGGAGTCAGGCCG‑CCCGGTG‑AAGCA‑GGAAACCACCGCAAAGATTTAGGCGGCAGcccgga                                                                                    >  1:976425/1‑144 (MQ=255)                                                                 tgGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGGTCTTTCCTGTGGAGGGATGGGGGAGTCAGGCCC‑CCCCGTG‑GAGCA‑GGGAACCCCCCAAAAGATTTCAGCGACCGCCCGaaagga                                                                                >  1:868846/1‑143 (MQ=255)                                                                  ggATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGCCGCCCGTCTTTTCCTTTGG‑GGGATGGGGGCGTAACGCCG‑CCCGTTG‑AAGGACGGAACCACCCGAAA‑GATTTAGGCGACAGcccggaaacggca                                                                             >  2:1036600/1‑138 (MQ=255)                                                                           gCGAACGCCGATGTAAATGGCGCTGGCAACATGTTAGCGGCGGGGCACG‑GCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGATCCCACCGAAATGATTCAGGCGCCAGGCTGAGCGTAGCAGGGATCCTCGt                                                                  <  2:834561/151‑1 (MQ=255)                                                                           gCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCCCC‑CCGCTCTTTCCCGTGGGGAGAAGGGGGAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCCAAAAGAATCAGGCGACAACCCTAACGTAGCAgggggacac                                                                    >  2:786376/1‑143 (MQ=255)                                                                                cGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTc                                                             <  1:305647/151‑1 (MQ=255)                                                                                   cGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGAAGCCCGTTTTTG‑CTGTGGGGAGGTGGGGGAGGCAGGGCG‑CCCCTTG‑AAACA‑GGAAACCACCGAAATGAATCAAGCGACAACCTGAACGTAACAGGGATGCCCGTCCTTCAgg                                                          >  1:298275/1‑151 (MQ=255)                                                                                    gATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGGTTTTTCCTCTGGGGAGAAGGAGGAGGCAGGGCG‑CCCGCTGAAACCA‑GGAACCACCCCAAAAGGTTCA‑GTGACAGCCGGAAAGGAACAGGGGTTCACCTCCTTCAcgg                                                         >  2:614652/1‑151 (MQ=255)                                                                                         aaaTGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTg                                                     <  2:1087872/150‑1 (MQ=255)                                                                                         aaaTGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTg                                                     <  2:372619/150‑1 (MQ=255)                                                                                          aaTGGCGCTCGCAACATTTTAGCGGCGGGGGACG‑CCGCTCTTTCCCGTGGAGAGAGGGGGGAGTCAGGGCG‑CCCCTTG‑AAGCA‑GGGAACCCCCCAAAAGATTTAGGCGACAGGCCGAAAACAACAGGGGTCCACCCCCTTCAGGGCggg                                                     >  1:1110167/1‑147 (MQ=255)                                                                                             ggCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTggagga                                                <  1:109160/151‑1 (MQ=255)                                                                                                     cAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAat                                         <  1:629510/150‑1 (MQ=255)                                                                                                        cATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCg                                     <  1:277886/151‑1 (MQ=255)                                                                                                            ttAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGtt                                  <  2:1145027/150‑1 (MQ=255)                                                                                                             tAGCAGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCTTGAACGTAGCAGGGATCCACGTGCTTCAGGGCGTGGAGGATGTCAATATCGGtt                                  <  2:1197144/149‑1 (MQ=255)                                                                                                                 ggcggGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGATG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAAAGCAGCAGGGATCCACGCACTTCAGGGCGCGGCGGCTGTCAAAATCGGttt                                 >  2:736624/1‑146 (MQ=255)                                                                                                                  gcggGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCa                           <  1:340576/151‑1 (MQ=255)                                                                                                                  gcggGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCa                           <  1:1189690/151‑1 (MQ=255)                                                                                                                      ggCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGa                       <  1:286962/151‑1 (MQ=255)                                                                                                                           g‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTaa                   >  2:278257/1‑150 (MQ=255)                                                                                                                                 tCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCa              >  1:954919/1‑150 (MQ=255)                                                                                                                                  cTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAac            >  2:400554/1‑151 (MQ=255)                                                                                                                                   ttGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAac            >  2:689599/1‑150 (MQ=255)                                                                                                                                    tGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCaa             <  2:327797/148‑1 (MQ=255)                                                                                                                                     gCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAAcacc         <  1:333432/151‑1 (MQ=255)                                                                                                                                         gtgGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGt     <  1:509722/151‑1 (MQ=255)                                                                                                                                           ggAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTaa   <  2:922097/151‑1 (MQ=255)                                                                                                                                           ggAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTa    >  1:1172134/1‑150 (MQ=255)                                                                                                                                           ggAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTa    <  1:440807/150‑1 (MQ=255)                                                                                                                                             agagaTGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAAc  >  2:848124/1‑150 (MQ=255)                                                                                                                                             |                                                                                                                                                         CGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACG‑CCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTG‑AAGCA‑GGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAAC  >  NC_000913/1504660‑1504948

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