BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,180,506	A→C	10.1%	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	10.1%	414.0 / 32.8	127	E29A (GAG→GCG) ‡	ygiC	ATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand): ref base A (51/63); new base C (5/8); total (56/71)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 127 total aligned reads displayed.

AAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGC  >  NC_000913/3180360‑3180651

 aaGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGcc                                                                                                                                                 <  2:618828/149‑1 (MQ=255)

 aaGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCg                                                                                                                                                >  4:511009/1‑149 (MQ=255)

    cACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATTCGGTTTCAATTTTCACACCATGTACGGCGAGCCGTa                                                                                                                                              <  2:606572/149‑1 (MQ=255)

     aCTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAc                                                                                                                                             <  4:577422/149‑1 (MQ=255)

         tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGgt                                                                                                                                        <  1:1105259/149‑1 (MQ=255)

         tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGgt                                                                                                                                        >  4:9274/1‑149 (MQ=255)

         tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGgt                                                                                                                                        <  3:43742/149‑1 (MQ=255)

         tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGgt                                                                                                                                        <  1:652219/149‑1 (MQ=255)

          gCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTACATGgtg                                                                                                                                       <  1:14602/149‑1 (MQ=255)

           cAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTACATGgtgt                                                                                                                                      <  1:1137693/149‑1 (MQ=255)

           cAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTACATGgtgt                                                                                                                                      <  3:511009/149‑1 (MQ=255)

            aGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGtgtg                                                                                                                                     >  4:208533/1‑149 (MQ=255)

               gTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGgt                                                                                                                                        >  4:613141/1‑143 (MQ=255)

                  gTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATg                                                                                                                               >  1:527745/1‑149 (MQ=255)

                  gTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTACATGGTGTGAAAAtt                                                                                                                               <  1:47695/149‑2 (MQ=38)

                   tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGc                                                                                                                              <  3:9274/149‑1 (MQ=255)

                   tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGc                                                                                                                              <  4:52777/149‑1 (MQ=255)

                   tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGc                                                                                                                              <  3:444234/149‑1 (MQ=255)

                   tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGc                                                                                                                              >  2:327786/1‑149 (MQ=255)

                   tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGc                                                                                                                              <  3:186860/149‑1 (MQ=255)

                      aaCCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTta                                                                                                                           >  2:488422/1‑149 (MQ=255)

                            tACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACaa                                                                                                                     <  4:718356/149‑1 (MQ=255)

                            tACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTAcatggt                                                                                                                                        <  4:154008/130‑6 (MQ=255)

                            tACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTAcatggt                                                                                                                                        >  3:154008/1‑125 (MQ=255)

                              ccTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGt                                                                                                                   >  2:141290/1‑149 (MQ=255)

                               ctctTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGtt                                                                                                                  <  4:60884/149‑1 (MQ=255)

                                tcttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTg                                                                                                                 >  2:535956/1‑149 (MQ=255)

                                    ctcGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGAccc                                                                                                             >  2:525638/1‑149 (MQ=255)

                                     tcGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGtt                                                                                                                  >  3:968757/1‑143 (MQ=255)

                                     tcGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCt                                                                                                            >  2:24896/1‑149 (MQ=255)

                                          cAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGccc                                                                                                       <  1:652661/149‑1 (MQ=255)

                                             tGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAgg                                                                                                    >  3:373768/1‑149 (MQ=255)

                                                gTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTg                                                                                                 >  4:549185/1‑149 (MQ=255)

                                                      gATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGc                                                                                           <  3:573518/149‑1 (MQ=255)

                                                       aTACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                          <  2:795848/149‑1 (MQ=255)

                                                        tACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTg                                                                                         <  4:569586/149‑1 (MQ=255)

                                                           cGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaa                                                                                      <  3:649854/149‑1 (MQ=255)

                                                             aTGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAg                                                                                                     >  4:90014/1‑132 (MQ=255)

                                                              tGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCTCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGCTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                                   <  3:801147/149‑1 (MQ=255)

                                                              tGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                                   <  3:208533/149‑1 (MQ=255)

                                                                gAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGt                                                                                 >  1:240849/1‑149 (MQ=255)

                                                                   agagTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                          <  1:782866/137‑1 (MQ=255)

                                                                       tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgcc                                                                          >  1:559544/1‑149 (MQ=255)

                                                                         aGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                          <  3:578138/131‑1 (MQ=255)

                                                                           tATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTta                                                                                                                           >  4:735978/1‑96 (MQ=255)

                                                                           tATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTta                                                                                                                           <  3:735978/96‑1 (MQ=255)

                                                                           tATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAAc                                                                      >  3:964356/1‑149 (MQ=255)

                                                                            aTTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACt                                                                     >  4:427057/1‑149 (MQ=255)

                                                                             ttACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTg                                                                    >  4:1126427/1‑149 (MQ=255)

                                                                             ttACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTg                                                                    >  2:909451/1‑149 (MQ=255)

                                                                                ccGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAc                                                                 <  1:842866/149‑1 (MQ=255)

                                                                                ccGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAc                                                                 <  1:535956/149‑1 (MQ=255)

                                                                                        ccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACt                                                                     <  2:157319/137‑1 (MQ=255)

                                                                                        ccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACt                                                                     >  1:157319/1‑137 (MQ=255)

                                                                                                gCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAg                                                 >  4:904348/1‑149 (MQ=255)

                                                                                                   tGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                              <  4:373768/149‑1 (MQ=255)

                                                                                                    gagaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                             >  2:395758/1‑149 (MQ=255)

                                                                                                     agaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggtg                                            <  3:983733/149‑1 (MQ=255)

                                                                                                     agaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggtg                                            >  2:516896/1‑149 (MQ=255)

                                                                                                     agaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggtg                                            <  4:867101/149‑1 (MQ=255)

                                                                                                       aaaGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGa                                          >  1:52904/1‑149 (MQ=255)

                                                                                                       aaaGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGa                                          >  1:907238/1‑149 (MQ=255)

                                                                                                           cccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaaa                                      <  1:24896/149‑1 (MQ=255)

                                                                                                            ccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGc                                                                                                                              >  4:429223/1‑60 (MQ=255)

                                                                                                            ccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGAccc                                                                                                             <  1:194841/77‑1 (MQ=255)

                                                                                                            ccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGAccc                                                                                                             >  2:194841/1‑77 (MQ=255)

                                                                                                            ccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAg                                     >  4:330641/1‑149 (MQ=255)

                                                                                                            ccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAg                                     >  4:1050847/1‑149 (MQ=255)

                                                                                                              aCGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                                   >  2:520388/1‑149 (MQ=255)

                                                                                                                gAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAt                                 <  1:141290/149‑1 (MQ=255)

                                                                                                                 aaTACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATc                                <  1:520388/149‑1 (MQ=255)

                                                                                                                       gTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAGGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAt                                 <  1:247025/142‑1 (MQ=255)

                                                                                                                       gTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAGGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAt                                 >  2:247025/1‑142 (MQ=255)

                                                                                                                       gTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                          <  2:907238/149‑1 (MQ=255)

                                                                                                                         ttCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAcc                                                                <  3:1067002/109‑1 (MQ=255)

                                                                                                                         ttCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAcc                                                                >  4:1067002/1‑109 (MQ=255)

                                                                                                                          tCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGAt                       <  3:427057/149‑1 (MQ=255)

                                                                                                                           cAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATg                      <  3:330641/149‑1 (MQ=255)

                                                                                                                             aTTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtg                                               <  3:531593/122‑1 (MQ=255)

                                                                                                                              tcttCACACCA‑GTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGc                   >  1:1052358/3‑149 (MQ=38)

                                                                                                                               tttCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt                  <  2:559544/149‑1 (MQ=255)

                                                                                                                                   acacCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATg                      <  1:466284/141‑1 (MQ=255)

                                                                                                                                    cacCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt                  >  1:638410/1‑144 (MQ=255)

                                                                                                                                    cacCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt                  <  2:638410/144‑1 (MQ=255)

                                                                                                                                     acCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGa                        <  1:605811/137‑1 (MQ=255)

                                                                                                                                     acCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgat               >  2:200147/1‑146 (MQ=255)

                                                                                                                                     acCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgat               <  1:200147/146‑1 (MQ=255)

                                                                                                                                     acCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAc            >  1:378356/1‑149 (MQ=255)

                                                                                                                                     acCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAc            >  2:532674/1‑149 (MQ=255)

                                                                                                                                      ccATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAt                                                            <  3:816698/100‑1 (MQ=255)

                                                                                                                                      ccATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAt                                                            >  4:816698/1‑100 (MQ=255)

                                                                                                                                       cATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtg                                                           >  3:418443/1‑100 (MQ=255)

                                                                                                                                       cATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATg                      <  3:839466/137‑1 (MQ=255)

                                                                                                                                       cATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATg                      >  4:839466/1‑137 (MQ=255)

                                                                                                                                        cagTACGGCTCGCCGTACATGGTGTGAAAATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaa                                                                                      >  3:728881/3‑73 (MQ=21)

                                                                                                                                        cagTACGGCTCGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAt                                                            <  4:610013/96‑1 (MQ=255)

                                                                                                                                        cagTACGGCTCGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAt                                                            >  3:610013/3‑98 (MQ=255)

                                                                                                                                        aTGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCaa         >  2:651413/1‑149 (MQ=255)

                                                                                                                                             cGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTcc    <  3:898470/149‑1 (MQ=255)

                                                                                                                                              ggCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCg   >  3:759237/1‑149 (MQ=255)

AAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGC  >  NC_000913/3180360‑3180651

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