Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,157,271	0	T	G	38.6%	148.7 / 13.3	88	G165G (GGT→GGG)	ycfH	putative DNase
Reads supporting (aligned to +/- strand):  ref base T (10/44);  new base G (34/0);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.82e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCAT  >  NC_000913/1157123‑1157418                                                                                                                                                      |                                                                                                                                                    atatCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑Taa                                                                                                                                                   <  2:431846/151‑1 (MQ=255)  tatCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑Taa                                                                                                                                                   >  1:236655/1‑150 (MQ=255)   atCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑Taaa                                                                                                                                                  <  2:64415/150‑1 (MQ=255)     ccAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAAtt                                                                                                                                                <  2:591682/150‑1 (MQ=255)        gATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTg                                                                                                                                            <  1:995625/151‑1 (MQ=255)        gATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTg                                                                                                                                            <  2:563651/151‑1 (MQ=255)        gATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑Gaaa                                                                                                                                                  >  1:1212918/1‑145 (MQ=255)         aTTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTgg                                                                                                                                           <  2:666147/151‑1 (MQ=255)          ttGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAAATACTgg                                                                                                                                           >  1:1219390/1‑150 (MQ=255)            ggTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGAt                                                                                                                                         <  1:172370/150‑1 (MQ=255)            ggTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAATTAATGGAtt                                                                                                                                        >  1:990328/1‑150 (MQ=255)                  gaactgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TGAATTATTTGATATCggt                                                                                                                                   >  2:1132524/1‑149 (MQ=255)                   aactgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGAtt                                                                                                                                 <  2:515524/151‑1 (MQ=255)                   aactgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGAtt                                                                                                                                 <  2:339753/151‑1 (MQ=255)                    actgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGGTCTCGGAttg                                                                                                                                <  2:394131/151‑2 (MQ=255)                      tgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAAATAATTGATCTCGGGtttt                                                                                                                               >  2:936763/1‑150 (MQ=255)                      tgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAATTACTGGGTCTCGGGtttt                                                                                                                               >  1:903414/1‑150 (MQ=255)                        aacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAATTACTGGATCTCGGAtttt                                                                                                                               >  1:212062/1‑148 (MQ=255)                          cAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAATTAATGGGTCTCGGGTTTTata                                                                                                                            >  1:591670/1‑147 (MQ=255)                           aaGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACAtct                                                                                                                         <  1:767219/151‑1 (MQ=255)                            aGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAAATACTGGGTCTCGGGtttt                                                                                                                               >  1:622347/1‑144 (MQ=255)                            aGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAAATACTGGGTCTCGGATTTTACAtctc                                                                                                                        >  2:539454/1‑151 (MQ=255)                            aGCCGGTTATCGTCCATACCCGTGACGCCCCTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TCAATTACTGGGATTCGGGtttt                                                                                                                               >  2:907752/1‑144 (MQ=255)                              ccGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCt                                                                                                                      <  2:1216705/151‑1 (MQ=255)                                ggTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCttt                                                                                                                    <  2:518848/151‑1 (MQ=255)                                ggTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCttg                                                                                                                    <  2:662253/151‑2 (MQ=255)                                  ttATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCtttt                                                                                                                   <  2:1197241/150‑1 (MQ=255)                                   tATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTc                                                                                                                  <  2:658482/150‑1 (MQ=255)                                   tATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTc                                                                                                                  <  1:492037/150‑1 (MQ=255)                                    aTCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTTGATCTCGGGTTTTACATCTCCtttt                                                                                                                   >  1:699007/1‑148 (MQ=255)                                      cGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCgg                                                                                                               <  1:126885/150‑1 (MQ=255)                                         ccATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGAGTTTATATCTCCTTTTTCCGCAtt                                                                                                           >  1:1323211/1‑151 (MQ=255)                                          cATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATtg                                                                                                          <  1:753917/151‑1 (MQ=255)                                          cATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGGTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATtg                                                                                                          <  2:1208045/151‑1 (MQ=255)                                           aTACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAAATACTGGATCCCCGATTTTATATCCCCTTTTCCCGCCtt                                                                                                           >  2:561468/1‑149 (MQ=255)                                             aCCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATtgtg                                                                                                        <  1:190267/150‑1 (MQ=255)                                               ccGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GGAATTACTTGATCTCGGGTTTTTCATCTCCCTTTTCGGCCTTGTGAcc                                                                                                     >  2:1167761/1‑151 (MQ=255)                                               ccGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAGATTTTTGGGGTTCGGGTTTTACAACTCCTTTTTCggtcat                                                                                                           >  2:113080/1‑141 (MQ=255)                                               ccGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAATTACTGGGTCTCGGGTTTTACATCTCCTTTTTCGGCATTGTGAcc                                                                                                     >  1:580303/1‑151 (MQ=255)                                                 gTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAATTTCTGGGTCTCGGATTTTAAATCTCCTTTTTCGGCATTGTGACCtt                                                                                                   >  1:525865/1‑151 (MQ=255)                                                    aCGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCg                                                                                                <  1:105836/151‑1 (MQ=255)                                                      gCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTa                                                                                              <  2:652798/151‑1 (MQ=255)                                                       cccGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTaa                                                                                             <  1:645435/151‑1 (MQ=255)                                                       cccGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTaa                                                                                             <  1:127065/151‑1 (MQ=255)                                                       cccGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAAATACTGGGTCTCGGGTTTTTAATCTCCTTTTTCCGCATTTTGACCTTCCCTa                                                                                              >  1:146815/1‑150 (MQ=255)                                                        cagTGCCGATACGCTGGCAATTCTGCTCGAAGAAAGAGTGACGGATTGCGGTGGCGTACTACGCTGTAGCACAGAGGACAGAGAAACGGCTTG‑TAAATTGCTGCATCTCTGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCTTAAt                                                                                            <  2:26253/149‑1 (MQ=255)                                                         cGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATg                                                                                           <  2:1313926/151‑1 (MQ=255)                                                              cGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAAATACTGGATCTCGGGTTTTAAAACTCCTTTTTCGGCATTGTGAACTTCCGTaa                                                                                             >  1:1092760/1‑144 (MQ=255)                                                               gATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTTATGGATCTCCGGTTTTTCATTTCCTTTTTCGGCATTGTGACCTTTCGCAATgaggag                                                                                      >  1:886604/1‑150 (MQ=255)                                                               gATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGc                                                                                     <  1:139304/151‑1 (MQ=255)                                                                aTACGCTGGCAATTCTGCGCGAAGAAAAAGCGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GGAAATACTTGCCCCCCCGCTTTACATCCCCCTTTTCCCCCTTGTGACCCCCCCCAATGCGGAGc                                                                                     >  2:638882/1‑150 (MQ=255)                                                                      tGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAAATAATTGGTCTCGGGTTTTACATCTCCTTTTTCCGCCTTGTGGCCTTTCCTTATGGGGAGGaa                                                                                   >  2:648523/1‑146 (MQ=255)                                                                       ggCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTgc                                                                               <  2:1236822/149‑1 (MQ=255)                                                                         cAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAAATACTGGGTCTCGGGTTTTTCATCTCCTTTTTCGGCATTGTGACCTTTCCTAAATCGGGGGAAATgcg                                                                              >  1:1259397/1‑148 (MQ=255)                                                                              cTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAATTACTGGGTCTCGGGGTTTACATCTCCTTTTTCGGCATTGTGACCTTCCGTAATGAGGAGCCACTGCGCGGtg                                                                         >  1:1026022/1‑148 (MQ=255)                                                                               tGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAATTAATGGGTCTCGGATTTTAACTCTCCTTTTTCCGCCCTGTGACCTTCCGTAATGCGGGGCAACTGCGCGAt                                                                          >  2:679939/1‑146 (MQ=255)                                                                                  gcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGcgcg                                                                  <  1:8622/151‑1 (MQ=255)                                                                                    gaagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTAACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCACGt                                                                 <  2:107745/150‑1 (MQ=255)                                                                                     aggaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGAGCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAAGGCGGAGCAACTGCGCGATGCTGCGCGttt                                                               <  2:347033/149‑2 (MQ=255)                                                                                      agaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAATTACTTGGTCTCGGGTTTTTTATCTCCTTTTTCCGCATTGTGACCTTCCCTAATGCGGGGGAACTGCGCGGTGCTGCGCGGTAt                                                              >  1:106416/1‑151 (MQ=255)                                                                                          aaaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAAATACTGTTTCTTGGATTTTTTATCTCCCTTTTCGCGGTTGTGACCCTTCTTAATTCgggggaa                                                                                   >  2:1198837/1‑121 (MQ=255)                                                                                           aaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAAATAATGGATCTCGGGTTTTACATCTCCTTTTCCCGCATTGTGACCTTCCCTAATTCGGAGCAACAGCGCGATGCTGCGCGTTATGTccc                                                         >  1:1353314/1‑151 (MQ=255)                                                                                            aGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GGGATTTCTGGAGATCGTCTTTTATACTTCCTTTTTTGGCATTGTGAAATTCCCTAATGTGGGGCAAATGCGAGCTGTGGCGCGTTTTGTcccc                                                        >  2:1178357/1‑151 (MQ=255)                                                                                              tGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAATTACTGGATGCCGGATTTTAAAACTCCTTTTTCGGCATTGTGGACTTCCGTAATGAGGAGGAACAGCGCGATGCTGCGCGGTATGTccccc                                                       >  2:1120677/1‑150 (MQ=255)                                                                                              tGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAATTAATGGATCTCGGGTTTTATATCTCCTTTTCCGGCATTTTGGCCTTCCGTAATGCGGAGCAAATGCGCGATGCTTCGCGGTTTGTccccc                                                       >  2:561804/1‑150 (MQ=255)                                                                                                aCGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTgg                                                    <  1:305872/151‑1 (MQ=255)                                                                                                aCGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GTAATTACTGGATCTCGGGTTTTAAATTTCCTTTTTCGGCGTTGTGACATTTCGTAATGCGGAGCAACTGAGCGCTGCTGCGCGTTATGTccccc                                                       >  2:994381/1‑148 (MQ=255)                                                                                                  ggATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGAt                                                  <  2:208972/151‑1 (MQ=255)                                                                                                      tGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCgg                                               <  2:1258693/150‑1 (MQ=255)                                                                                                       gCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGtt                                             <  1:438309/151‑1 (MQ=255)                                                                                                       gCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGt                                              <  2:424333/150‑1 (MQ=255)                                                                                                       gCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAATTACTTGGTCTTCGATTTTAAATCTCCCTTTTCCGCATTTTGAACTTCCCTAATTCGGAGCCAATTAGCGATGCTGCGCGGTATTTCCCCCCTGATCgg                                               >  1:751174/1‑149 (MQ=255)                                                                                                              cGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑GAAATTACTGGATCTCGGGTTTTAAAACTCCTTTTCCGGGATTGTGACCTTTCCTAATGCGGGGGAAATACGCGATTTTGCGCGTTTTGTCCCCCCGGGTCGGGTACt                                          >  2:885094/1‑147 (MQ=255)                                                                                                                  ctacACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAAATACTTGATCTCCGACTTTAAATCTCCTTTTTCGGCCTTTTGGCCTTCCTTAAAGCGGGGCAAATGCGCGCTGCTGCGCGTTTTGTCCCCCCCGATCGGTTACtggtggga                                   >  2:1240958/1‑148 (MQ=255)                                                                                                                        tGTTTTACAGAGGACAGAGAAACGGCGGGTTAAATTAC‑GGGTTTCCGATTTTATATTTCCTTTTTCGGGCTTTTTATCTTCCGTAATGCGGGGCAAACGCTCGGTGGTGCTCGGTATGTTCCCCTCGATCGGttgc                                           >  1:181617/1‑134 (MQ=255)                                                                                                                           tttACAGAGGACAGAGAAACGGCGGG‑GAAAATACTTGGTTTCTGTTTTTTAATCACCTTTTCCGGCATTATGAACTTCCGTGATGACGGGGACCCGCGCGCAGCTGCGCGGTTTGTGCCCCCGGGTCGGGTCCTGGTGGGAACAGaaccac                         >  2:687865/1‑146 (MQ=255)                                                                                                                               cAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCacctta                     <  1:955258/151‑1 (MQ=255)                                                                                                                               cAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCacctta                     <  2:1212203/151‑1 (MQ=255)                                                                                                                                agagGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCacctta                     <  1:700140/150‑1 (MQ=255)                                                                                                                                agagGACAGAGAAACGGCGGG‑GAAATTACTTGGTCTCGGGTTTTAAAACTCCTTTTTCGGCATTGTGGCCTTCCGTAATGCGGAGCAACCGCGCGGTGCTGCGCGTTATGTCCCCCCGGGTCGGGTACTGGTGGGAACTGACTCacctt                      >  1:473222/1‑149 (MQ=255)                                                                                                                                 gagGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCacctt                      <  1:813603/148‑1 (MQ=255)                                                                                                                                 gagGACAGAGAAACGGCGGG‑GAAAATAATGGGGCTCGGGTTTTAAAACTCCCTTTTCGGCATTGTGACACTCCGTAATTCGGAGGAAACGCGCGACGCTGCGCGGTTTGTCCCCCCGGGTCGGGTACTGGGGGGAAATGAcatacc                        >  2:813032/1‑141 (MQ=255)                                                                                                                                     aCAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTgc               <  1:926911/151‑1 (MQ=255)                                                                                                                                       agagAAACGGCGGT‑GAAATTAATGGGTCTCGGGTTTTAAAATTCCTTTTCCGGGCTTTTGACCCTCCGTAATGCGGAGCAAACGCGCGATGCTGCGGGTTTTTTGCCCCCCGGTCGGTTACTGGTGGAAACTGACTCcccttacccttc               >  2:431108/1‑145 (MQ=255)                                                                                                                                           aaaCGGCGGG‑GAAATTACTGGATCTCGGGTTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTTTGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCCCCTTACCTTGCGcccg          >  2:976672/1‑148 (MQ=255)                                                                                                                                                 cggG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCa   <  1:1011292/151‑1 (MQ=255)                                                                                                                                                  ggg‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCAt  <  2:580173/151‑1 (MQ=255)                                                                                                                                                  ggg‑TAAATTACTGGATATCGGATTTTACATCTCCTTTTCCCGCATTGTGACCTTCCGTAATGCGGATCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCa   >  2:998747/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                    ATATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGG‑TAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCAT  >  NC_000913/1157123‑1157418

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