Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	659,934	0	T	G	44.6%	122.1 / 11.9	74	T95P (ACC→CCC)	lipA	lipoate synthase
Reads supporting (aligned to +/- strand):  ref base T (5/36);  new base G (33/0);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 4.47e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCGCGCAGGTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGA  >  NC_000913/659786‑660082                                                                                                                                                     |                                                                                                                                                       ccATCGCGCAGGTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTAc                                                                                                                                                      <  2:944849/151‑1 (MQ=255)   aTCGCGCAGGTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACa                                                                                                                                                     <  1:312993/150‑1 (MQ=255)           ggTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACaaaaaaagcc                                                                                                                                            >  1:567583/1‑144 (MQ=255)              catcaCGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACaaaaa                                                                                                                                                 >  2:213267/1‑141 (MQ=255)               atcaCGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACaaaaaaagccccgg                                                                                                                                        >  1:1098498/1‑140 (MQ=255)                 caCGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGa                                                                                                                                       <  1:1024531/150‑1 (MQ=255)                 caCGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGAt                                                                                                                                      <  2:1072835/151‑1 (MQ=255)                  aCGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATc                                                                                                                                     <  1:504247/151‑1 (MQ=255)                   cGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAGCGCCGAGGGTc                                                                                                                                     >  1:395945/1‑150 (MQ=255)                    ggTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCAt                                                                                                                                   <  1:522439/151‑1 (MQ=255)                    ggTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGGCAAAAAGCGCCGCgg                                                                                                                                        >  1:304107/1‑146 (MQ=255)                    ggTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGCCaaaaaaagccgcggggc                                                                                                                                     >  2:256133/1‑135 (MQ=255)                        aaCGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAAc                                                                                                                               <  1:348520/151‑1 (MQ=255)                        aaCGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGGCAAAAAGAGCCGGGGGTTATaaa                                                                                                                                >  2:61710/1‑150 (MQ=255)                        aaCGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGGCAAAAAAAGCCGAGGGACATaaa                                                                                                                                >  2:941424/1‑150 (MQ=255)                          cGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGt                                                                                                                             <  2:812912/151‑1 (MQ=255)                          cGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAGCGCCGAGGGGCATaaa                                                                                                                                >  1:847695/1‑148 (MQ=255)                           ggAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGtt                                                                                                                            <  2:579209/151‑1 (MQ=255)                               gTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCtg                                                                                                                        <  1:215848/151‑1 (MQ=255)                               gTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCGAGGATCATAAACAATGCt                                                                                                                         >  1:952202/1‑150 (MQ=255)                                    aaCCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAGAGCCGAGGGGCATAAAAGTTGCtgt                                                                                                                       >  1:978512/1‑147 (MQ=255)                                       cacaTAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTgg                                                                                                                <  2:363829/151‑1 (MQ=255)                                        acaTAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGAAAAAAAGCGCCCGGGGAAATAAAAATTGCTGTGCCgggg                                                                                                                >  1:264616/1‑150 (MQ=255)                                          aTAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAAAGCCCCGGGGCATAAACATTGCTGTGGCGTgg                                                                                                                >  1:513299/1‑148 (MQ=255)                                            aaCGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGaa                                                                                                           <  1:895692/151‑1 (MQ=255)                                             aCGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAGAGCCGCGGAAAATAAAAGTTGCTGTTCCgggggt                                                                                                              >  2:319218/1‑145 (MQ=255)                                              cGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTATCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑c                                                                                                        <  2:190266/151‑1 (MQ=255)                                              cGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAAAAACGCCGAGGATCATTAAACTTGCTGTGCCgggggt                                                                                                              >  1:21567/1‑144 (MQ=255)                                                 aGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CAtt                                                                                                     <  2:27981/151‑1 (MQ=255)                                                   cgcCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCGAGGGGCATAAAAGTTGCTGTGGCGTGGGTGGAG‑CATTcc                                                                                                   >  1:428277/1‑151 (MQ=255)                                                     ccATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAATAACGACCCGGATACAAAAAATTGTTGTGCCGGGGGTGGAA‑CAtttc                                                                                                   >  1:668803/1‑147 (MQ=255)                                                     ccATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCGAGGGGCATAAAAATTTGTGTGGCGGGGGTGGAG‑AATTcccc                                                                                                 >  1:769776/1‑149 (MQ=255)                                                     ccATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGAAAAATAATATCGCGGGAAATAAACGTTTGTGTTCCGTCGTTGaggagt                                                                                                       >  2:616707/1‑141 (MQ=255)                                                              cAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAgg                                                                                       <  1:58624/151‑1 (MQ=255)                                                                aTGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGc                                                                                     <  2:1007138/151‑1 (MQ=255)                                                                  ggTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCa                                                                                    <  2:239362/150‑1 (MQ=255)                                                                       gcgcCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCaggag                                                                                <  2:592244/149‑1 (MQ=255)                                                                       gcgcCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCagga                                                                                 <  2:19475/148‑1 (MQ=255)                                                                        cgcCAGTTTCACTGGTTCATTGGCATC‑GGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGTCAAAAAGCGCCGGGGTGCAAAAAAATTGGTGTGCCGTGGGTGAAG‑GATTTCG‑CCCGGTTTGGGCgggagg                                                                               >  1:193775/1‑148 (MQ=255)                                                                         gcCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGGCAAAAAGCGCCCAGGGTCATAAAAAATGCTGTGGCGGGGGGGGAG‑CATTCCG‑CCCGGGTTGGGGagggggggt                                                                            >  2:132489/1‑148 (MQ=255)                                                                               ttCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCg                                                                       <  1:396902/150‑1 (MQ=255)                                                                                tCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGGCAAAAAACGGCCAGGGACATAAAAGTTGCTGTGGCCCGGGTGGAG‑CATTTCGCCCGGGTTAGGG‑AGGGGGGTTTCTCGc                                                                      >  1:1050457/1‑150 (MQ=255)                                                                                 cACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGca                                                                     <  1:616733/150‑1 (MQ=255)                                                                                 cACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGca                                                                     <  2:395929/150‑1 (MQ=255)                                                                                 cACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCCAGGAGCATAAACATTGCTGTGCCGTGGGTTAAG‑CAATCCG‑CCCGGGTAGGGGGGGAGGcgttc                                                                          >  1:74730/1‑141 (MQ=255)                                                                                     ggTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGcagac                                                                  <  1:91110/149‑1 (MQ=255)                                                                                       tacATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACGTCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGcagacaga                                                               <  1:961204/148‑1 (MQ=255)                                                                                          aTTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAGCGCCCACGATCATAAACCATGCTGTGCCGTGGGTGGAG‑GATTTCG‑CCCGGGTAGGGCCGGGGGGTTCCTCCCCGCCAGaa                                                              >  1:464849/1‑148 (MQ=255)                                                                                                  cAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCAtt                                                   <  1:400792/151‑1 (MQ=255)                                                                                                  cAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCAtt                                                   <  2:288383/151‑1 (MQ=255)                                                                                                  cAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAAAAACGCCCAGGAGCCTAAAAATTCCGTTGCCGTGGGTGGAG‑GATTCCG‑CCCGGGTAGGGCAGGAGGGTTCCTCCCCGACAGaaaaaagggccct                                                   >  2:289088/1‑140 (MQ=255)                                                                                                  cAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAAAAAAGCCGAGGGACATAAAAATTGCTGTGGCGGGGGTGAAG‑CATTCCC‑CCCGGGTAGGGGAGGAGGCTTCCTCGCAGACAGAAAGCAgggcctt                                                   >  1:45140/1‑146 (MQ=255)                                                                                                   aGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCGAGGGTCATAAAAGTTGCTGTGGCGGGGGTGAAG‑GAATCCG‑CCCGGGTTGGGGGGGGGGGGTCCTCCCCGACAGAAAGCAGGCCCttt                                                  >  1:545012/1‑151 (MQ=255)                                                                                                      aGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTg                                               <  1:978965/151‑1 (MQ=255)                                                                                                      aGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTg                                               <  1:661299/151‑1 (MQ=255)                                                                                                         tACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCGAGGATCATAAAAGTTGCTGTGCCGTGGGTTAAG‑CATTTCG‑CCCGGGTTGGGGAGGGGGGGTCCTCGCCGACAGAAAGCAGGCCCttttt                                                >  1:358865/1‑147 (MQ=255)                                                                                                          aCCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAGAGCCGGGGATCATAAAAGTTTGTGTGCCGGGGGTGAAG‑CAATCCG‑CCCGGGTAGGGGAGGGGGGGTCCTCGCCGACAGAATGGAGGGCAttttt                                                >  1:763182/1‑146 (MQ=255)                                                                                                             gggCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAGCGCCGCGGGGGATAAAAGCTGCTGTGCCCCGGGTTTAG‑CACTTCC‑CCCGGCTAGGGGGGGGGGGTTCCCCCCCGGCAGAAAGCAGGGCCtttttttgca                                           >  1:767339/1‑143 (MQ=255)                                                                                                                 gACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACTGCGAGGAGCATAAAAGCTGCTGTGCCGTGGGTGAAG‑AATTCCG‑CCACGTTGGGGGCGGAGGCGTTCCCGCCGAAAAAAAACAGGGCACTTTTTTGCAt                                          >  2:374626/1‑145 (MQ=255)                                                                                                                    cGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCttt                                 <  2:1018126/151‑1 (MQ=255)                                                                                                                    cGTGGGCAACATCACAGAACGGACAACGGCGGGGAGAAAAAAAGCCCAGGAACATAAAAATTGCTGTGCCGGGGGTGGAGACATT‑CG‑CCCGGGTAGGGCAGGGGGGTTCCTCCCCGACAGAAAGGAGGGCCtttttttggattttgggttt                                 >  1:879513/1‑136 (MQ=255)                                                                                                                     gTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCCCGGGGCACAAAAAGTGCTGTGGCGGGGGTGGAA‑CAATCCG‑CCCCGGTTGGGGAGGGGGGGTCCCCCCCGGCAGAATGCAGGCCCCTTTTTCGCATTTCggggtt                                 >  2:1018529/1‑146 (MQ=255)                                                                                                                     gTGGGCAACATCACAGAACGGACAACGGCGGGGAAAAATAAGGCCGCGGATCATAAAAATCGCTGTGCCCCGGGTGTAG‑CATTCCC‑CCCCGCTTGGGCAGGGGGGTTCCTCCCCGACCGAAAAGCGGCCAATTTTTTGCATTTCGGCctt                                 >  2:532743/1‑150 (MQ=255)                                                                                                                         gCAACATCACAGAACGGACAACGGCGGGGACAAAAAAAGCCGAGGGTCATAAAAATTGGTGTGCCGGGGGTGGAG‑GATTCCC‑CCCGGGGAGGGGGGGGGGGTTCCTCCCCGACAGAAAGAAGGGCATTTTTTTGCATTTCggggttttt                              >  2:959306/1‑142 (MQ=255)                                                                                                                           aaCATCACAGAACGGACAACGGCGGGGACAAAAAGCGCCGAGGGTCATAAAAGGTTTTGTGGCGTGGGTTGAG‑CATTCCG‑CCCGGGTAGGGGAGGGGGGGTCCTCGCAGACAGAATGCAGGGCATTTTTTTTCATTGCGGCGTTTATGccc                          >  1:570369/1‑151 (MQ=255)                                                                                                                            aCATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCt                         <  2:548059/151‑1 (MQ=255)                                                                                                                            aCATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCt                         <  2:1033160/151‑1 (MQ=255)                                                                                                                                   aGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACg                  <  1:729398/151‑1 (MQ=255)                                                                                                                                     aaCGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACgtg                <  2:978475/151‑1 (MQ=255)                                                                                                                                        ggACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTag             <  1:1063143/151‑1 (MQ=255)                                                                                                                                          aCAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTagag           <  1:1024355/151‑1 (MQ=255)                                                                                                                                              cggcggGTACAAAAAGCGCCGAGGATCATAAAAGTTTCTGTGCCGGGGGTGAAG‑CATTCCG‑CCCGGTTAGGGCAGGAGGCTTCCTCGCCGACAGAAAGCAGGCCATTTTTTCGCATTGCGGCGTTGTTGCCCTGGATACGGGTAGAGTcc        >  2:586256/1‑150 (MQ=255)                                                                                                                                                gcggGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCt     >  2:746772/1‑151 (MQ=255)                                                                                                                                                   ggTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGa  <  2:450529/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       CCATCGCGCAGGTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGTGCCGTGGTTGAAG‑CATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCAGGCCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGA  >  NC_000913/659786‑660082

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