Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,133,302	0	T	G	46.8%	99.8 / 12.7	63	H219P (CAC→CCC)	yghS	putative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (5/28);  new base G (29/0);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 7.17e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGG  >  NC_000913/3133152‑3133450                                                                                                                                                       |                                                                                                                                                     gAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgt                                                                                                                                                      <  1:858964/151‑1 (MQ=255)       aGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGc                                                                                                                                                 <  2:855684/150‑1 (MQ=255)                   gcgcGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCaa                                                                                                                                    <  1:264847/151‑1 (MQ=255)                       ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCGCCCATCCAGTTCAAgg                                                                                                                                  >  1:948997/1‑149 (MQ=255)                        gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATc                                                                                                                               <  2:598643/151‑1 (MQ=255)                        gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCaa                                                                                                                                    >  2:56519/1‑146 (MQ=255)                             aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCATCCAGTTCAAAGATtttt                                                                                                                            >  2:379897/1‑149 (MQ=255)                             aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCTCCCATCCCGTTCCAGGGTCttt                                                                                                                            >  1:137158/1‑149 (MQ=255)                                aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTgcg                                                                                                                         <  1:435568/149‑1 (MQ=255)                                   aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGGTCAAGGATCTTTGCGCCAtt                                                                                                                    >  1:864277/1‑151 (MQ=255)                                      agtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTTAAAGATCTTTTCGccc                                                                                                                      >  1:684600/1‑145 (MQ=255)                                         ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATg                                                                                                               <  2:797722/150‑1 (MQ=255)                                          tgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGt                                                                                                              <  2:611486/150‑1 (MQ=255)                                             gtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCACTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGAGCTTTGCGCCATTGAATGTCAg                                                                                                           <  1:934989/150‑1 (MQ=255)                                               ttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGtt                                                                                                         <  1:523754/150‑1 (MQ=255)                                               ttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTGCGCCATTTAATGTCAGtt                                                                                                         >  1:553648/1‑150 (MQ=255)                                               ttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGGTCAAGGAGCTTTGCGCCATTTAATGTAAg                                                                                                           >  1:253822/1‑148 (MQ=255)                                                tgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGGCCCCTCCCGTTCAAGGATCTTTGCGCCATTGGATGTCAGTTg                                                                                                        >  1:462966/1‑150 (MQ=255)                                                            aTTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGATCTTTGCGCCATTTAATGTCAGTTGCGGCGTAACGGCg                                                                                           >  1:586334/1‑151 (MQ=255)                                                            aTTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGAGCTTTTCGCCATTGTAAGTCAGTTGCGGCGTAAAGGCg                                                                                           >  2:491754/1‑151 (MQ=255)                                                             ttCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAACGGCg                                                                                           >  2:2626/1‑150 (MQ=255)                                                              tCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGATCTTTGCGCCCTTGAATGTCAGGTGCGGCGCAACgg                                                                                             >  2:875120/1‑147 (MQ=255)                                                                gTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGGCCATCCAGGTTCAGGGTCTTTGTGGCCTTGTATGTCAGTTGGGGCGTAACGGCGGttt                                                                                       >  1:784898/1‑151 (MQ=255)                                                                   aaaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGAtttttt                                                                                    <  1:191270/151‑1 (MQ=255)                                                                           cccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAAc                                                                            <  1:1012984/151‑1 (MQ=255)                                                                            ccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCCAGGATCTTTTCGCCATTTAATGTCCGTTGCGGCGGAACGGGGGtttttt                                                                                    >  1:378897/1‑142 (MQ=255)                                                                               aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCt                                                                         <  1:240186/150‑1 (MQ=255)                                                                                   aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAg                                                                    <  1:254623/151‑1 (MQ=255)                                                                                    ttttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCGCCCATCCCGTTCAAAGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAACGGCGGTTTTTTTCTGTAACGCTGcc                                                                      >  2:988061/1‑148 (MQ=255)                                                                                     tttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCAGTTCAAGGATCTTTTTGCCCTTGAATGTCAGTTGGGGCGTAAAGGCGATTTTTTTCttttaagccgcccgt                                                                   >  1:395454/1‑136 (MQ=255)                                                                                     tttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGt                                                                   <  2:154884/150‑1 (MQ=255)                                                                                         gTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtggtg                                                              <  1:247231/151‑1 (MQ=255)                                                                                          tGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGATCTTTGCGCCATTGAATGTCAGGTGCGGCGTAACGGCGATTTTTTTCTGTAACGCTGCCCGTtggggg                                                             >  2:34788/1‑151 (MQ=255)                                                                                           gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTTAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGGAACGGCGATTTTTTTCTGTAACGCTGCCAGTtggggg                                                             >  2:687136/1‑150 (MQ=255)                                                                                             tCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTc                                                           <  1:465107/150‑1 (MQ=255)                                                                                               aGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGATATTTTCGCCATTTAATGGCAGGTGCGGCGGAACGGCGATTTTTTTCTGTTAAGCTGCCCGGTGGGGGGCAgg                                                        >  1:634359/1‑151 (MQ=255)                                                                                               aGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGGAACGGCGATTTTTTTCTTTTACGCTGCCCGGtggtggggagg                                                        >  1:367720/1‑146 (MQ=255)                                                                                                  aaaGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGttt                                                     <  1:619631/151‑1 (MQ=255)                                                                                                  aaaGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGTTCAAGGATCTTTTCGCCATTTAATGTCAGTTGCGGCGGAACGGCGGTTTTTTTCTGTTAAGCTGCCCGTTGGTGGTGAGGgtt                                                     >  1:880537/1‑151 (MQ=255)                                                                                                    aGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTTAAGGGACTTTGCGCCATTTAATGTCAGTTGCGGCGGGACGGCGATTTTTTTTTTTAACGCTGCCCGGTGGTGGGCAGGgtttt                                                   >  1:406735/1‑149 (MQ=255)                                                                                                      cgcgGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGGAACGGCGATTTTTTTCTGTAACGCTGCCCGTTGGTGGGCAGGgtttag                                                  >  2:728035/1‑147 (MQ=255)                                                                                                          gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCaa                                             <  2:333752/151‑1 (MQ=255)                                                                                                          gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCaa                                             <  1:624869/151‑1 (MQ=255)                                                                                                          gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCaa                                             <  2:98273/151‑1 (MQ=255)                                                                                                           cgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCaaa                                            <  2:668084/151‑1 (MQ=255)                                                                                                           cgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAAAGCTGCCAGTTGGGGGTCAGGGTTTCGCGCaa                                             >  2:148106/1‑150 (MQ=255)                                                                                                             tGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGGAACGGCGAATTTTTTCTGTAACGCTGCCAGTTGGTGGTCAGGGTTAAGCGCcaa                                            >  2:361083/1‑149 (MQ=255)                                                                                                              gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAAcgc                                         <  2:381965/151‑1 (MQ=255)                                                                                                                       gTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCACCAGGTTAACGGTCTTTGCGCCACTGTATGTGAGCTTCGGCGGGACGGCGATTTTTTCCTGTTACGCTGCCAGGTGGGGGGTAGGGTTTCGCGCAAAAGCGGgttgtt                                 >  2:711707/1‑150 (MQ=255)                                                                                                                         aGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGGATGTTAGTTGCGGCGGAACGGGGATTTTTTTCTGTAACGCTGCCAGTTGGTGGTGAGGTTTACGCGCAAACGCGGGTTGTtcat                              >  1:21468/1‑151 (MQ=255)                                                                                                                         aGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCAACCAGGTCAAGGATCTTTTCGCCCTTGAATGTCAGGTGCGGCGGAACGGCGATTTTTTCCTTTAACGCTGCCAGGTGGTGGGGAGGGTTACGCGCAAAAGCGGgttgtt                                 >  2:978651/1‑148 (MQ=255)                                                                                                                            gACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCaa                           <  2:367793/151‑1 (MQ=255)                                                                                                                            gACGCTTGCAGAATTTCATCGGCGGGGGGTAGCCCCTCCAGTTTTAGGATCTTTTCGCCCTTGAAAGGCAGTTGCGGCGTTACGGCGAATTTTTTCTTTAATGCTGCCAGTTGGGGGTCAGGGTTTCGCGCAAACGCGGgtt                                    >  2:532952/1‑142 (MQ=255)                                                                                                                                 ttGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGGTCTTTTTGCCATTTAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTAACGCTGCCCGTTGGTGGTTAGGGTTACGCGCAAAAGCGGGTTTTTTATCCATGcc                       >  1:302506/1‑150 (MQ=255)                                                                                                                                   gCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCaa                     <  1:158100/150‑1 (MQ=255)                                                                                                                                   gCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAg                    <  1:178362/151‑1 (MQ=255)                                                                                                                                      gAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACg                 <  2:137206/151‑1 (MQ=255)                                                                                                                                        aTTTCATCGGCGGGGGGGCGCCCATCCAGTTCAAGGGTCTTTGCGCCATTGAATGTCAGTTGCGGCGGAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGGTTACGCGCAAACGCGGGTTGTTCATCAATGCCAAGACg                 >  1:162774/1‑149 (MQ=255)                                                                                                                                           tCATCGGCGGGGCGTCGGCCATCCAGTTTAAGGATCTTTGCGCCATTGAATGTTCGTTGCGGCGTTACGGCGGTTTTTTCCTGTAACGCTGCCAGTTGGGGGTCAGGGTTACGCGCAAACGCGGTTTGTTTATCAATGCCCAGACGAATc             >  1:183188/1‑150 (MQ=255)                                                                                                                                               cggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaaca         <  2:365308/150‑1 (MQ=255)                                                                                                                                                 gcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATa      <  1:834501/151‑1 (MQ=255)                                                                                                                                                  cggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATAc     <  2:996792/151‑1 (MQ=255)                                                                                                                                                    gggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACgg   >  1:603987/1‑151 (MQ=255)                                                                                                                                                     ggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACggg  >  2:591104/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                     GAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGG  >  NC_000913/3133152‑3133450

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: ‑