BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,878,894	T→A	7.4%	D409V (GAC→GTC)	gyrB ←	DNA gyrase, subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,878,894	0	T	A	7.4%	423.7 / 18.5	122	D409V (GAC→GTC)	gyrB	DNA gyrase, subunit B
Reads supporting (aligned to +/- strand): ref base T (66/47); new base A (5/4); total (71/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 125 total aligned reads displayed.

GAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTA  >  NC_000913/3878746‑3879042

gAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGt                                                                                                                                                      <  3:1434613/149‑1 (MQ=255)

 aaGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTc                                                                                                                                                     <  1:788173/149‑1 (MQ=255)

  aGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCt                                                                                                                                                    <  2:497095/149‑1 (MQ=255)

  aGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCt                                                                                                                                                    <  1:668767/149‑1 (MQ=255)

        cTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGACTGCCAgg                                                                                                                                              <  2:1005403/149‑2 (MQ=255)

         tttCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGtt                                                                                                                                             >  1:1351667/1‑149 (MQ=255)

         tttCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGtt                                                                                                                                             >  1:242645/1‑149 (MQ=255)

          ttCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGACTGCCAggaa                                                                                                                                            <  4:1118957/149‑4 (MQ=255)

             tcGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGcc                                                                                                                                         >  4:480659/1‑149 (MQ=255)

             tcGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGcc                                                                                                                                         >  3:227523/1‑149 (MQ=255)

               gACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCg                                                                                                                                       >  1:1195927/1‑149 (MQ=255)

                  gTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCa                                                                                                                                    <  1:193548/149‑1 (MQ=255)

                       ggATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGccc                                                                                                                               >  3:1197111/1‑149 (MQ=255)

                          ttttACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCt                                                                                                                            <  4:227523/149‑1 (MQ=255)

                             tACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGttt                                                                                                                                            >  2:828041/1‑130 (MQ=255)

                             tACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGttt                                                                                                                                            <  1:828041/130‑1 (MQ=255)

                             tACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAg                                                                                                                         >  3:687166/1‑149 (MQ=255)

                              aCCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGt                                                                                                                        <  2:733609/149‑1 (MQ=255)

                              aCCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGt                                                                                                                        >  2:1057875/1‑149 (MQ=255)

                                ccTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCg                                                                                                                      >  4:1356078/1‑149 (MQ=255)

                                ccTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCg                                                                                                                      >  2:708299/1‑149 (MQ=255)

                                  ttCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAg                                                                                                                    <  3:921940/149‑1 (MQ=255)

                                  ttCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAg                                                                                                                    >  1:848640/1‑149 (MQ=255)

                                   tCAGCTGCAGAATCGCCTGGTTCTTGCTGTTACTTCTTTGCTTCGCAGAGCCGCCCGCGGAAACTCGTTCCACCAGGTACGCATCGGAAAGCGCCGGATCGCGTCCCCAGCGGTCTGCCAGTTTGCCTGCCACGCCCGCTAAGTCGAgc                                                                                                                   <  4:470504/149‑1 (MQ=255)

                                             aaTCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAc                                                                                                         >  4:473930/1‑149 (MQ=255)

                                             aaTCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAc                                                                                                         <  2:455237/149‑1 (MQ=255)

                                                 gCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAcggcg                                                                                                     >  2:129180/1‑149 (MQ=255)

                                                  ccTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAcggcgg                                                                                                    <  1:818813/149‑1 (MQ=255)

                                                  ccTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTAcggcgg                                                                                                    <  4:1193885/149‑1 (MQ=255)

                                                  acTGGTACTTGCGGTTACGCCCCTGCTTCGCAGAGCCGGCCGCGGAGTCCCCTTCCACCCGGTACAGTTCGGAAAGCGCCGGATCGTGTTTCTGGCAGTCTGGCAGTTTGACCGGCAGGCCCGCTAAGTCGAGCGCACTTTTAccgcgg                                                                                                    >  4:1192092/2‑149 (MQ=255)

                                                   cTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCggg                                                                                                   >  3:786752/1‑149 (MQ=255)

                                                   cTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCggg                                                                                                   >  1:109532/1‑149 (MQ=255)

                                                   cTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCggg                                                                                                   >  1:1084250/1‑149 (MQ=255)

                                                    tGGTTATTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTGCAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGt                                                                                                  >  4:652795/1‑149 (MQ=255)

                                                                ttACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCAcgcgc                                                                                      >  2:965735/1‑149 (MQ=255)

                                                                    gCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCAtt                                                                                              >  3:682328/1‑137 (MQ=255)

                                                                    gCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCAtt                                                                                              <  4:682328/137‑1 (MQ=255)

                                                                      cccTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGt                                                                                                                                              <  1:980689/87‑1 (MQ=255)

                                                                      cccTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGAcg                                                                                >  4:336470/1‑149 (MQ=255)

                                                                      cccTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGACTGCCAggaa                                                                                                                                            >  1:1321989/1‑86 (MQ=255)

                                                                      cccTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGACTGCCAggaa                                                                                                                                            <  2:1321989/89‑4 (MQ=255)

                                                                       ccTGCTTCGCAGAGCCGCCCGCGGAGTCGCCATCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCGGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCCCGCGCGCGAcgc                                                                               >  4:893723/1‑149 (MQ=255)

                                                                        cTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGAcgcg                                                                              >  1:945806/1‑149 (MQ=255)

                                                                        cTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGAcgcg                                                                              >  3:1046499/1‑149 (MQ=255)

                                                                        cTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGAcgcg                                                                              >  2:916482/1‑149 (MQ=255)

                                                                            ttCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGcgccgc                                                                          >  3:190790/1‑149 (MQ=255)

                                                                               gCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTc                                                                       >  2:1439699/1‑149 (MQ=255)

                                                                               gCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTc                                                                       >  2:553770/1‑149 (MQ=255)

                                                                                  gagCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTcacg                                                                    <  4:1046499/149‑1 (MQ=255)

                                                                                   agCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTCTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTcacgg                                                                   <  4:469965/149‑1 (MQ=255)

                                                                                      cgcccgcGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTcacgggca                                                                >  4:1312452/1‑149 (MQ=255)

                                                                                          cgcgGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTcacgggcacggg                                                            <  1:1439699/149‑1 (MQ=255)

                                                                                          cgcgGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTcacgggcacggg                                                            <  3:1001080/149‑1 (MQ=255)

                                                                                           gcgGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCa                                                                                                                <  1:443589/96‑1 (MQ=255)

                                                                                               aGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCATCTTTACGGCGGGTCATTTCACGCGCGCGACACGCCGCTTCACGGGCACGGgcagca                                                       <  2:703607/149‑1 (MQ=255)

                                                                                               aGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGgcagca                                                       <  1:570077/149‑1 (MQ=255)

                                                                                                 tCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATc                                                     >  2:396201/1‑149 (MQ=255)

                                                                                                   cccTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCAc                                                                                          >  4:150127/1‑110 (MQ=255)

                                                                                                   cccTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCAc                                                                                          >  3:343329/1‑110 (MQ=255)

                                                                                                   cccTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCAc                                                                                          <  4:343329/110‑1 (MQ=255)

                                                                                                   cccTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGa                                                   >  1:126853/1‑149 (MQ=255)

                                                                                                    ccTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCCCGGGCACGGGCAGCATCGAt                                                  >  3:873411/1‑149 (MQ=255)

                                                                                                    ccTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGAt                                                  >  1:412447/1‑149 (MQ=255)

                                                                                                     cTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCCTCGATa                                                 >  1:248545/1‑149 (MQ=255)

                                                                                                     cTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATa                                                 >  3:746192/1‑149 (MQ=255)

                                                                                                     cTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATa                                                 >  1:1005324/1‑149 (MQ=255)

                                                                                                     cTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATa                                                 >  1:862461/1‑149 (MQ=255)

                                                                                                      ttCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATaa                                                <  4:361627/149‑1 (MQ=255)

                                                                                                      ttCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATaa                                                >  4:772024/1‑149 (MQ=255)

                                                                                                         caccaGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAAttt                                             >  3:363790/1‑149 (MQ=255)

                                                                                                          accaGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAAtttt                                            <  3:1336331/149‑1 (MQ=255)

                                                                                                           ccaGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTg                                           >  1:38647/1‑149 (MQ=255)

                                                                                                             aGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGcc                                         >  1:463460/1‑149 (MQ=255)

                                                                                                                  cAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCa                                    <  4:86480/149‑1 (MQ=255)

                                                                                                                   aGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCAc                                   >  2:324329/1‑149 (MQ=255)

                                                                                                                    gTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCACTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACg                                  >  2:862618/1‑149 (MQ=255)

                                                                                                                      tCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGAt                                <  4:1197111/149‑1 (MQ=255)

                                                                                                                        ggAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGAttt                              >  4:820573/1‑149 (MQ=255)

                                                                                                                          aaaGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTc                            <  3:1191977/149‑1 (MQ=255)

                                                                                                                           aaGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCCACCACGATTTTcg                           >  1:1409976/1‑149 (MQ=255)

                                                                                                                            aGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTcgc                          <  1:1057875/149‑1 (MQ=255)

                                                                                                                            aGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTcgc                          <  3:1356078/149‑1 (MQ=255)

                                                                                                                            aGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTcgc                          >  2:325019/1‑149 (MQ=255)

                                                                                                                              cgcCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGt                        <  2:412447/149‑1 (MQ=255)

                                                                                                                                 cGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCgg                     <  1:240936/149‑1 (MQ=255)

                                                                                                                                 cGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCgg                     <  3:652795/149‑1 (MQ=255)

                                                                                                                                    aTCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTg                  >  3:619613/1‑149 (MQ=255)

                                                                                                                                     tCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCCCGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTgg                 >  4:1035114/1‑149 (MQ=255)

                                                                                                                                     tCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTgg                 >  4:921350/1‑149 (MQ=255)

                                                                                                                                       gcgTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGt               >  2:458769/1‑149 (MQ=255)

                                                                                                                                          aaacTGGCAGACTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGtttt            >  4:213029/4‑149 (MQ=255)

                                                                                                                                           tCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTc           <  4:1116721/149‑1 (MQ=255)

                                                                                                                                           tCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCCCGTCGGTTGGGGTTTc           >  4:632527/1‑149 (MQ=255)

                                                                                                                                           aacTGGCAGACTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTc           >  1:1267479/3‑149 (MQ=255)

                                                                                                                                             cTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCca         <  2:972656/149‑1 (MQ=255)

                                                                                                                                             cTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCca         >  2:1426583/1‑149 (MQ=255)

                                                                                                                                              tGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTTGCGGCGGTTGGGTTTTCcag        >  1:873074/1‑149 (MQ=255)

                                                                                                                                                gCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcagca      <  3:1268715/149‑1 (MQ=255)

                                                                                                                                                gCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcagca      <  1:1130644/149‑1 (MQ=255)

                                                                                                                                                    tCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTa  >  4:324195/1‑149 (MQ=255)

GAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTA  >  NC_000913/3878746‑3879042

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