BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,180,505	G→T	5.7%	E29* (GAG→TAG) ‡	ygiC →	ATP‑Grasp family ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,180,505	0	G	T	5.7%	1258.2 / 45.6	341	E29* (GAG→TAG) ‡	ygiC	ATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand): ref base G (153/168); new base T (14/6); total (167/174)
Fisher's exact test for biased strand distribution p-value = 6.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 347 total aligned reads displayed.

AAACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATT  >  NC_000913/3180356‑3180654

    aaaGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAg                                                                                                                                                      >  2:1176350/1‑148 (MQ=255)

    aaaGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGcc                                                                                                                                                    >  2:2897492/1‑150 (MQ=255)

        cACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAc                                                                                                                                                <  2:4634638/150‑1 (MQ=255)

          cTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑Tg                                                                                                                                             <  2:2468345/150‑1 (MQ=255)

            aTGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGgt                                                                                                                                           >  2:2740018/1‑150 (MQ=255)

            aTGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGgt                                                                                                                                           >  2:1489516/1‑150 (MQ=255)

             tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCTCGCCGTAcatggt                                                                                                                                           <  1:3660692/150‑6 (MQ=255)

             tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCTCGCCGTAcatggt                                                                                                                                           <  1:2727605/150‑6 (MQ=255)

               cAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGtgtg                                                                                                                                        >  2:4527256/1‑150 (MQ=255)

               cAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGtgtg                                                                                                                                        >  1:874756/1‑150 (MQ=255)

               cAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGtgtg                                                                                                                                        >  2:3511520/1‑150 (MQ=255)

                 gCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGaa                                                                                                                                      <  2:599800/150‑1 (MQ=255)

                 gCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGaa                                                                                                                                      <  1:4322051/150‑1 (MQ=255)

                  cGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAg                                                                                                                                     >  1:913399/1‑150 (MQ=255)

                  cGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAg                                                                                                                                     >  1:1591792/1‑150 (MQ=255)

                   gTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGa                                                                                                                                    >  2:1267228/1‑150 (MQ=255)

                   gTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGa                                                                                                                                    >  2:4283553/1‑150 (MQ=255)

                      gTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGc                                                                                                                                 <  2:2217131/150‑1 (MQ=255)

                      gTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGc                                                                                                                                 <  2:3924973/150‑1 (MQ=255)

                       tGCAACCGGGACCTCTTCTCGTTCAATGGGTGGCTGATACCGAAGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCt                                                                                                                                <  2:1658630/150‑1 (MQ=255)

                        gCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCtt                                                                                                                               >  1:3922610/1‑150 (MQ=255)

                         cAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTta                                                                                                                              <  2:1566935/150‑1 (MQ=255)

                           aCCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTtact                                                                                                                            >  2:3975975/1‑150 (MQ=255)

                           aCCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTtact                                                                                                                            <  1:3267329/150‑1 (MQ=255)

                                  ccTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGtt                                                                                                                     >  2:194334/1‑150 (MQ=255)

                                     cttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGAc                                                                                                                  <  1:4091596/150‑1 (MQ=255)

                                      ttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACGGAGCGCCCGGAGTGGCGTGAGAAAGCCCACGTATACGGTTTG‑AATTTTCACACCATGTACGGCGAGCAGTAC‑TGGTGTGAAGTTGCTTACTACAAGTTGAcc                                                                                                                 <  2:948303/150‑1 (MQ=255)

                                      ttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGAcc                                                                                                                 >  1:1178937/1‑150 (MQ=255)

                                      ttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGAcc                                                                                                                 <  1:3078948/150‑1 (MQ=255)

                                            ttCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGcc                                                                                                           <  2:3266348/150‑1 (MQ=255)

                                            ttCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACAACATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGcc                                                                                                           >  1:2699249/1‑150 (MQ=255)

                                                aTGGGTGGCTGATACCGATGGTAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAgg                                                                                                       <  2:2988191/150‑1 (MQ=255)

                                                aTGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAgg                                                                                                       <  1:2172610/150‑1 (MQ=255)

                                                      ggCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGaaa                                                                                                 <  2:4122485/150‑1 (MQ=255)

                                                        cTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAg                                                                                               >  1:987676/1‑150 (MQ=255)

                                                         tGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGc                                                                                              >  2:2773650/1‑150 (MQ=255)

                                                            tACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTgg                                                                                           <  1:4030480/150‑1 (MQ=255)

                                                            tACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAACATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTgg                                                                                           >  2:1441518/1‑150 (MQ=255)

                                                              ccGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaa                                                                                         <  1:2602522/150‑1 (MQ=255)

                                                                 aTGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                                      <  1:521678/150‑1 (MQ=255)

                                                                  tGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGaagaag                                                                                     <  2:2865612/150‑1 (MQ=255)

                                                                     aaaGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCa                                                                                  >  1:1647576/1‑150 (MQ=255)

                                                                         agTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgc                                                                              <  1:2498279/150‑1 (MQ=255)

                                                                         agTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgc                                                                              <  1:4536278/150‑1 (MQ=255)

                                                                           tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgccg                                                                            >  1:1314316/1‑150 (MQ=255)

                                                                           tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgccg                                                                            >  1:3373348/1‑150 (MQ=255)

                                                                           tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgccg                                                                            <  2:1734023/150‑1 (MQ=255)

                                                                               tATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACt                                                                        >  1:3129858/1‑150 (MQ=255)

                                                                                   aCCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAc                                                                    >  2:1993641/1‑150 (MQ=255)

                                                                                   aCCGAGCGCCCGGACTGGCGCGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAc                                                                    <  1:205329/150‑1 (MQ=255)

                                                                                       aGCGCCCGGACTGGCGTGAGACAGCCCACGAATCCGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGa                                                                <  1:1094733/150‑1 (MQ=255)

                                                                                       aGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGa                                                                >  1:2530816/1‑150 (MQ=255)

                                                                                            ccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                                           >  1:4687809/1‑150 (MQ=255)

                                                                                            ccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                                           <  2:101765/150‑1 (MQ=255)

                                                                                            ccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                                           <  1:2450129/150‑1 (MQ=255)

                                                                                            ccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                                           <  2:4739764/150‑1 (MQ=255)

                                                                                              ggACTGGCGTGAGAAAGCCCACGAATACGGTTTCAAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTCGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGcc                                                          >  1:3515504/1‑150 (MQ=255)

                                                                                                  tGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGaaa                                                     <  2:1046510/150‑1 (MQ=255)

                                                                                                    gCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAACAGCTGGAAGAAGTCACCGCCGACCTGCACCAGATGTGCCTGAAAgt                                                   >  2:3243843/1‑150 (MQ=255)

                                                                                                    gCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgt                                                   >  1:2899623/1‑150 (MQ=255)

                                                                                                      gTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                                 >  2:4452052/1‑150 (MQ=255)

                                                                                                      gTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                                 >  2:2073840/1‑150 (MQ=255)

                                                                                                      gTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                                 >  1:2860959/1‑150 (MQ=255)

                                                                                                       tGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                                <  2:437793/150‑1 (MQ=255)

                                                                                                       tGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                                <  2:256009/150‑1 (MQ=255)

                                                                                                       ggagaAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                                <  2:266629/149‑1 (MQ=255)

                                                                                                          gaAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGa                                             >  2:163233/1‑150 (MQ=255)

                                                                                                              gCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaaa                                         <  2:1829046/150‑1 (MQ=255)

                                                                                                              gCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaaa                                         <  1:4659629/150‑1 (MQ=255)

                                                                                                                ccACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGt                                       <  1:2333902/150‑1 (MQ=255)

                                                                                                                      aTACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGc                                 <  1:3431229/150‑1 (MQ=255)

                                                                                                                       tACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGcc                                >  2:2861550/1‑150 (MQ=255)

                                                                                                                        aCGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACCAGTTGACCCTCGCCCAGGTTGAATAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCa                               <  2:4554468/150‑1 (MQ=255)

                                                                                                                        aCGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGGTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCa                               >  2:3450399/1‑150 (MQ=255)

                                                                                                                         cGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAg                              <  2:4608463/150‑1 (MQ=255)

                                                                                                                         cGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAg                              <  2:2331768/150‑1 (MQ=255)

                                                                                                                         cGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAg                              >  1:794275/1‑150 (MQ=255)

                                                                                                                         cGGTTTC‑AATGTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAg                              <  2:3137695/150‑1 (MQ=255)

                                                                                                                          ggTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                             <  2:4319720/150‑1 (MQ=255)

                                                                                                                           gTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCg                            <  2:3654782/150‑1 (MQ=255)

                                                                                                                            tttC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGa                           >  1:1499979/1‑150 (MQ=255)

                                                                                                                               c‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGa                        >  1:511453/1‑150 (MQ=255)

                                                                                                                                 aaTTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAg                       >  1:640155/1‑150 (MQ=255)

                                                                                                                                  atcttCACACCA‑GTACGGCTCGCCGTACATGGTGTGAAAATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACAGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGc                      >  1:1749639/4‑150 (MQ=38)

                                                                                                                                  aTTTTCCCACCCTGTACGGCGAGCCGTAC‑TGGTGTGAAGAGGCTGACTCCCCGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGc                      <  2:867350/150‑1 (MQ=255)

                                                                                                                                   atttCACACCATGTACGGCGAGCCGTAC‑TGGTGAGAAGATGATAACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt                     <  1:2212759/149‑1 (MQ=255)

                                                                                                                                      tCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgat                  <  2:3845363/150‑1 (MQ=255)

                                                                                                                                        acacCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgatga                >  2:2462387/1‑150 (MQ=255)

                                                                                                                                        acacCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgatga                >  2:768489/1‑150 (MQ=255)

                                                                                                                                        acacCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgatga                <  1:3870656/150‑1 (MQ=255)

                                                                                                                                          acCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAcc              >  2:4086156/1‑150 (MQ=255)

                                                                                                                                           accagTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAcc              >  1:1087245/5‑150 (MQ=255)

                                                                                                                                           accagTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAcc              >  1:67800/5‑150 (MQ=255)

                                                                                                                                           accagTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAcc              >  2:1129676/5‑150 (MQ=255)

                                                                                                                                           accagTACGGCTCGCCGTACATGGTGTGAAAATGCTTACTACAACTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAcc              >  1:4438949/5‑150 (MQ=255)

                                                                                                                                            cATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCaa            >  2:4678476/1‑150 (MQ=255)

                                                                                                                                                tACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTc        >  1:3184775/1‑150 (MQ=255)

                                                                                                                                                  cGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCg      <  2:496843/150‑1 (MQ=255)

                                                                                                                                                     cGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCAt   <  1:4487201/150‑1 (MQ=255)

                                                                                                                                                      gAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCAtt  >  1:2817320/1‑150 (MQ=255)

AAACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATT  >  NC_000913/3180356‑3180654

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`