BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,058,621	T→A	5.9%	I72N (ATC→AAC) ‡	typA →	GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,058,621	0	T	A	5.9%	395.6 / 10.1	136	I72N (ATC→AAC) ‡	typA	GTP‑binding protein
Reads supporting (aligned to +/- strand): ref base T (65/63); new base A (3/5); total (68/68)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

Only 100 of 137 total aligned reads displayed.

ACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGA  >  NC_000913/4058474‑4058770

aCAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCg                                                                                                                                                     >  1:1181742/1‑150 (MQ=255)

  aaGCTGCTCCAACAATCCGGTACGTTCGCCTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGtt                                                                                                                                                   >  1:611942/1‑150 (MQ=255)

   aGCTGCTCCACCATTCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTg                                                                                                                                                  <  2:1326153/150‑1 (MQ=255)

   aGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTg                                                                                                                                                  >  2:736639/1‑150 (MQ=255)

   aGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTg                                                                                                                                                  >  1:213242/1‑150 (MQ=255)

   aGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTg                                                                                                                                                  <  2:1075123/150‑1 (MQ=255)

   aGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTg                                                                                                                                                  <  2:1653442/150‑1 (MQ=255)

   aGCTGCTCCAAAATTCCGTTACGTGCGACTCTCGTGCCGAAACCCAAGAGAGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTg                                                                                                                                                  <  2:348660/150‑1 (MQ=255)

     ctgctCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGAt                                                                                                                                                >  1:664031/1‑150 (MQ=255)

       gctCCAACAAACCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACGATGTTGATAc                                                                                                                                              <  2:210136/150‑1 (MQ=255)

        ctCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACGATGTTGATAcg                                                                                                                                             <  2:909136/150‑2 (MQ=255)

         tCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACTCCGCTATCAAATGGAATGATTACCGTATCAACGATGTTGATAcgg                                                                                                                                            <  1:619963/150‑3 (MQ=255)

         tCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACGATGTTGATAcgg                                                                                                                                            <  1:1566557/150‑3 (MQ=255)

         tCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATAccc                                                                                                                                            >  2:338966/1‑150 (MQ=255)

          ccAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACGATGTTGATAcggt                                                                                                                                           <  2:262360/150‑4 (MQ=255)

          ccAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATAcccc                                                                                                                                           <  1:736639/150‑1 (MQ=255)

            aataaTCCGGTAGGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTTGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCGATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCgg                                                                                                                                         <  1:675442/147‑1 (MQ=255)

             acaaTCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCggg                                                                                                                                        <  1:123488/150‑1 (MQ=255)

              caaTCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCgggg                                                                                                                                       <  2:56561/150‑1 (MQ=255)

               aaTCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGc                                                                                                                                      <  1:1489449/150‑1 (MQ=255)

                aTCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCa                                                                                                                                     >  2:1431883/1‑150 (MQ=255)

                aTCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTAGATACCCCGGCGCa                                                                                                                                     <  2:1606356/150‑1 (MQ=255)

                   cGGTACGTTCGACTCTCGTGCCGAAAGCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGc                                                                                                                                  <  2:89752/150‑1 (MQ=255)

                    ggTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTCCCATCCTCGCGAAAAACACCGCTATCCAATGGAATGATTACCGTATCAACATCGTTGATCCCCCGGGGCCCggc                                                                                                                                 >  1:1080020/1‑148 (MQ=255)

                        cGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATTCCCCGGGGCACGCCGACt                                                                                                                             >  1:1278407/1‑150 (MQ=255)

                         gTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACtt                                                                                                                            >  2:94993/1‑150 (MQ=255)

                           tCGACGCTCGTGCCGAAACCCAAGAGCGCGTGATGGACGCCACCGATTTGGAGAAAGAGCGTGGGAGTACCATCCTCGCGAAATACACCGCTATCAAATGGAATGATTACAGTAGCAACATCGTTGATACCCCGGGGCACGCCGACTTCg                                                                                                                          <  1:257027/150‑1 (MQ=255)

                             gACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATAc                                                                                                                                              >  2:248190/1‑128 (MQ=255)

                              aCTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGCCTTCggtg                                                                                                                       >  1:1086386/1‑150 (MQ=255)

                               ctctCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCggtgg                                                                                                                      <  2:1462906/150‑1 (MQ=255)

                                 ctcGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGgtggtg                                                                                                                    >  1:72410/1‑150 (MQ=255)

                                  tcGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGa                                                                                                                   <  1:742074/150‑1 (MQ=255)

                                  tcGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGGGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGa                                                                                                                   <  2:1532286/150‑1 (MQ=255)

                                   cGTGCCGAAACCCAAGCGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGCTGATACCCCGGGGCACGCCGACTTCGGTGGTGaa                                                                                                                  <  2:164563/150‑1 (MQ=255)

                                   cGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGaa                                                                                                                  <  2:1187124/150‑1 (MQ=255)

                                    gTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAg                                                                                                                 >  1:1222609/1‑150 (MQ=255)

                                     tGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGt                                                                                                                >  2:496803/1‑150 (MQ=255)

                                        cgcaaCCCAAGAGCGCGGGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGACTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGCGCACGCCGACTTCGGTGGTGAAGTTGa                                                                                                             <  2:1339722/147‑1 (MQ=255)

                                           aaCCCAAGAGCGCGTGATGGACTCCACCGATTTGGAGAAAGAGCGTGGGATTACCCTCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCCCGCCGACTTCGGTGGGGAAGTTGCACg                                                                                                          >  2:1080325/1‑150 (MQ=255)

                                             cccAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACgtg                                                                                                        >  2:329884/1‑150 (MQ=255)

                                              ccAAGAGCGCGTGATGGACTCCAACGATTTGTAGAAAGAGCGTGGGATTACCAACCTCGCGAAAAACACCGCTATCAAATGGCAGGATTACCGTCTCACCATCGTTGATACCCCGGGGCACGACGACTTCGGTGGTGAAGTTGAACttgt                                                                                                       >  1:306552/1‑150 (MQ=255)

                                              ccAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACgtgt                                                                                                       >  1:399064/1‑150 (MQ=255)

                                                 agagCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAAt                                                                                                    >  1:1429544/1‑150 (MQ=255)

                                                     cgcgTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTcc                                                                                                <  2:664031/150‑1 (MQ=255)

                                                     cgcgTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTcc                                                                                                >  1:558158/1‑150 (MQ=255)

                                                     cgcgTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTcc                                                                                                >  1:541111/1‑150 (MQ=255)

                                                       cgTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCAt                                                                                              >  1:248349/1‑150 (MQ=255)

                                                       cgTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCAt                                                                                              >  1:939687/1‑150 (MQ=255)

                                                         tGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATgg                                                                                            <  1:329884/150‑1 (MQ=255)

                                                                 tCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCa                                                                                    <  2:939687/150‑1 (MQ=255)

                                                                  ccAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCATCATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAg                                                                                   <  2:611942/150‑1 (MQ=255)

                                                                    aaCGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGtc                                                                                 >  2:177638/1‑149 (MQ=255)

                                                                    aaCGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTAAGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGtg                                                                                 >  2:319457/1‑150 (MQ=255)

                                                                     aCGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGtgc                                                                                <  2:544404/150‑1 (MQ=255)

                                                                     aCGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGtgc                                                                                <  2:964685/150‑1 (MQ=255)

                                                                       gATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGtgctg                                                                              <  1:1439368/150‑1 (MQ=255)

                                                                         atttggtaaaGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTgctgct                                                                            <  1:553204/143‑1 (MQ=255)

                                                                           tGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCtgg                                                                          >  1:1597844/1‑150 (MQ=255)

                                                                                  agagCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGa                                                                   <  2:1608870/150‑1 (MQ=255)

                                                                                   gagCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGAc                                                                  <  1:724578/150‑1 (MQ=255)

                                                                                    agCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACg                                                                 <  1:500778/150‑1 (MQ=255)

                                                                                       gTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCAt                                                              >  2:1231594/1‑150 (MQ=255)

                                                                                       gTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCAt                                                              >  1:1073830/1‑150 (MQ=255)

                                                                                        ggggATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGATGATACCCCGGGGCACGCCGACTTCGGGGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGGGGTTGACGCAtt                                                             <  2:1399305/149‑1 (MQ=255)

                                                                                            aTTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGAc                                                         <  1:95591/150‑1 (MQ=255)

                                                                                             ttACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACg                                                        >  1:183293/1‑150 (MQ=255)

                                                                                              tACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGTACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACgg                                                       >  2:784344/1‑150 (MQ=255)

                                                                                              tACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACgg                                                       >  2:110142/1‑150 (MQ=255)

                                                                                                ccATCCTCTCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGcc                                                     <  1:230237/150‑1 (MQ=255)

                                                                                                 cATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTCTCAACATCGTTGATACCCCGGGGCACGCCGACTGCGGTGGTGAAGTTGAACGGGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTGGACGCATTTGACGGccc                                                    <  1:1655722/150‑1 (MQ=255)

                                                                                                    ccTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGAt                                                 >  1:905656/1‑150 (MQ=255)

                                                                                                     cTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTTACGCATTTGACGGCCCGATg                                                >  1:1366652/1‑150 (MQ=255)

                                                                                                     cTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGagg                                                >  2:1036525/1‑148 (MQ=255)

                                                                                                       cgcgAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGcc                                              >  2:384084/1‑150 (MQ=255)

                                                                                                       cgcgAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGcc                                              >  2:1398213/1‑150 (MQ=255)

                                                                                                       cgcgAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGcc                                              <  1:1686623/150‑1 (MQ=255)

                                                                                                          gAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCa                                           >  2:1130560/1‑150 (MQ=255)

                                                                                                             aaaTACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGACGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAAc                                        <  1:319457/150‑1 (MQ=255)

                                                                                                                cacCGCTATCAAATGGAATGATTACCGTATCAACATCGATGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGATGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGCCCCGATGCCGCAAAcgcg                                     >  2:1683322/1‑150 (MQ=255)

                                                                                                                  ccGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGCCGCATTTGACGGCCCGATGCCGCAAACgcgcg                                   >  1:258198/1‑149 (MQ=255)

                                                                                                                  ccGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTTACGGCCCGATGCCGCAACCGCGCt                                   >  2:280659/1‑150 (MQ=255)

                                                                                                                  ccGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCt                                   >  1:638106/1‑150 (MQ=255)

                                                                                                                         cAAATTGAATTATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAAc                            <  1:995366/150‑1 (MQ=255)

                                                                                                                          aaaTGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAAcc                           <  2:248349/150‑1 (MQ=255)

                                                                                                                          aaaTGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAAcc                           <  2:541111/150‑1 (MQ=255)

                                                                                                                              ggAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCaaaa                       >  2:787224/1‑150 (MQ=255)

                                                                                                                                 aTGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAg                    <  2:1322710/150‑1 (MQ=255)

                                                                                                                                   gATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCg                  >  1:158423/1‑150 (MQ=255)

                                                                                                                                    aatACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGGTGAACGTGTAATTTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCTCTTCGTAACCAAAAAAGCGt                 <  2:56071/148‑1 (MQ=255)

                                                                                                                                     ttACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCAATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGAATCCTAACCAAAAAAGCGtt                >  2:286624/1‑150 (MQ=255)

                                                                                                                                       ggggTATCAACGATGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTg              >  1:914083/4‑150 (MQ=255)

                                                                                                                                       aCCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTg              >  1:1418364/1‑150 (MQ=255)

                                                                                                                                        ccGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGc             >  1:708396/1‑150 (MQ=255)

                                                                                                                                         cGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCt            >  2:521031/1‑150 (MQ=255)

                                                                                                                                           tATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTa          <  2:558158/150‑1 (MQ=255)

                                                                                                                                               aaCATCGTTGATACCCCGGGGCAGGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGGCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGc      <  2:1060001/150‑1 (MQ=255)

                                                                                                                                               aaCATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGc      <  2:1237541/150‑1 (MQ=255)

                                                                                                                                                ccATCGTTGATACCCCGGGGCACGCCGACTTCGTGGGGGAAGTTGACCGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGcc     <  1:1148597/149‑1 (MQ=255)

                                                                                                                                                  aTCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTg   <  2:1311696/150‑1 (MQ=255)

                                                                                                                                                   tCGTTGATACCCCGGGGCACGCCGACTTCGGGGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGa  <  1:1242394/150‑1 (MQ=255)

ACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGA  >  NC_000913/4058474‑4058770

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