BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,180,506	A→C	6.0%	E29A (GAG→GCG) ‡	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,506	0	A	C	6.0%	1203.7 / 53.5	342	E29A (GAG→GCG) ‡	ygiC	ATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand): ref base A (155/166); new base C (14/7); total (169/173)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 345 total aligned reads displayed.

AACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTC  >  NC_000913/3180357‑3180655

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

      gCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTa                                                                                                                                                  <  2:1599971/150‑1 (MQ=255)

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

         cTATTCAGCGTAGTGCAAGCGGTAACTCTTCTAGTTCTATGGGTGGCTGATAACGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCTCGCCGTAcat                                                                                                                                               <  1:3005530/150‑3 (MQ=255)

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

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

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

            tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGgtg                                                                                                                                           >  2:4380417/1‑150 (MQ=255)

             gCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCATTGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGgtgt                                                                                                                                          <  1:2039812/150‑1 (MQ=255)

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

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

               aGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGa                                                                                                                                        <  2:2049741/150‑1 (MQ=255)

               aGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGa                                                                                                                                        <  2:2824284/150‑1 (MQ=255)

               aGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGa                                                                                                                                        <  2:2465186/150‑1 (MQ=255)

               aGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGa                                                                                                                                        <  2:342997/150‑1 (MQ=255)

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

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

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

                   tAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTTATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGAt                                                                                                                                    <  1:3403709/150‑1 (MQ=255)

                   tAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGAt                                                                                                                                    <  1:3239389/150‑1 (MQ=255)

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

                      tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCt                                                                                                                                 <  2:1526661/150‑1 (MQ=255)

                      tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTAC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCt                                                                                                                                 <  2:1658610/150‑1 (MQ=255)

                        cAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTta                                                                                                                               >  1:2716395/1‑150 (MQ=255)

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

                              gTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCCGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACaa                                                                                                                         <  2:4532626/150‑1 (MQ=255)

                              gTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACaa                                                                                                                         <  2:1044339/150‑1 (MQ=255)

                              gTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACaa                                                                                                                         >  2:2397634/1‑150 (MQ=255)

                                 ccTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGtt                                                                                                                      <  1:2850706/150‑1 (MQ=255)

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

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

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

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

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

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

                                                         gATACCGATGGTATTAGTTAGTATTACCGAGCGCCCGGTCTGGCTTGAGAAAGCCCACGAATACTGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGATGTTGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                              <  1:2517783/150‑1 (MQ=255)

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

                                                              cGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaag                                                                                         >  2:3217339/1‑150 (MQ=255)

                                                               gATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaaga                                                                                        >  2:2279635/1‑150 (MQ=255)

                                                                ctGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                                       <  2:2142246/149‑1 (MQ=255)

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

                                                                   gAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTc                                                                                    >  1:1647482/1‑150 (MQ=255)

                                                                    aaaGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCa                                                                                   <  1:4014520/150‑1 (MQ=255)

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

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

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

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

                                                                            aGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGaa                                                                           >  2:1207909/1‑150 (MQ=255)

                                                                            aGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTCCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGaa                                                                           >  1:447293/1‑150 (MQ=255)

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

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

                                                                                           ccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCCGATGTGc                                                            >  1:3907046/1‑150 (MQ=255)

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

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

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

                                                                                            cGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGcc                                                           <  2:2662396/150‑1 (MQ=255)

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

                                                                                             ggACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTCGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCt                                                          >  1:4808197/1‑150 (MQ=255)

                                                                                             ggACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTCGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCt                                                          >  1:4305833/1‑150 (MQ=255)

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

                                                                                                     gTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGACAgtgg                                                  >  1:2811477/1‑150 (MQ=255)

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

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

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

                                                                                                      tGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                                 >  2:2960304/1‑150 (MQ=255)

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

                                                                                                       gagaAAGCCCACGACTACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggtg                                                <  1:1689133/150‑1 (MQ=255)

                                                                                                              cccACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAg                                         <  1:771439/150‑1 (MQ=255)

                                                                                                               ccACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGt                                        >  1:2597233/1‑150 (MQ=255)

                                                                                                                  cGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAt                                     <  2:2545825/150‑1 (MQ=255)

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

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

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

                                                                                                                        cGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAg                               <  2:4608463/150‑1 (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‑AATTTTCACACCATGTACGGCGAGCCGTAG‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGTAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGATAGTGGTGGAAAAAGTGTTCGCCAGCGa                            >  2:171107/1‑150 (MQ=255)

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

                                                                                                                                 atcttCACACCA‑GTACGGCTCGCCGTACATGGTGTGAAAATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGc                       >  1:3640154/4‑150 (MQ=38)

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

                                                                                                                                 aTTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGc                       <  1:2479224/150‑1 (MQ=255)

                                                                                                                                     tCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgat                   >  2:3858549/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                                                                              gTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAAt           >  1:911015/1‑150 (MQ=255)

                                                                                                                                              gTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAAtt          <  2:192769/150‑1 (MQ=255)

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

                                                                                                                                                aCGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTcc        <  1:4074239/150‑1 (MQ=255)

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

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

                                                                                                                                                      aGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTc  <  1:1430502/150‑1 (MQ=255)

AACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTC‑AATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTC  >  NC_000913/3180357‑3180655

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