BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,058,622	C→T	5.8%	I72I (ATC→ATT) ‡	typA →	GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,058,622	0	C	T	5.8%	441.4 / 9.3	137	I72I (ATC→ATT) ‡	typA	GTP‑binding protein
Reads supporting (aligned to +/- strand): ref base C (65/64); new base T (3/5); total (68/69)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-01

Only 100 of 138 total aligned reads displayed.

ACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAA  >  NC_000913/4058474‑4058771

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

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

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

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

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

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

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

   aGCTCCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTTGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTg                                                                                                                                                   <  2:271802/150‑1 (MQ=255)

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

        ctCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACCCCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATAcc                                                                                                                                              >  2:1083685/1‑150 (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)

            aacaaTCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCgg                                                                                                                                          >  1:730867/1‑150 (MQ=255)

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

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

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

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

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

                   cGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGc                                                                                                                                   <  1:1328036/150‑1 (MQ=255)

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

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

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

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

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

                               ctctCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCggtgg                                                                                                                       <  2:1214537/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)

                                       ccGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGCGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTg                                                                                                               >  2:1242394/1‑150 (MQ=255)

                                         gAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGaa                                                                                                             <  1:912313/150‑1 (MQ=255)

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

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

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

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

                                              ccAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACgtgt                                                                                                        <  2:72410/150‑1 (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:541111/1‑150 (MQ=255)

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

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

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

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

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

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

                                                                    aaatcgttGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGtg                                                                                  >  2:1327620/7‑150 (MQ=255)

                                                                    aaCGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGGTGAACGTGTAATGTCCATGGTAGACTCAGtg                                                                                  <  2:1685120/150‑1 (MQ=255)

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

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

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

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

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

                                                                               gaAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGTATGATTACCGTATCAACATCGTTGATACCCCGGGTATCGCCGACTTCGGTGGTGAAGTTGACCGTGTAATGTCCATGGTAGACTCAGTGCTGCTggtggt                                                                       >  2:1039291/1‑150 (MQ=255)

                                                                               gaAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTggtggt                                                                       >  1:1598902/1‑150 (MQ=255)

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

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

                                                                                     gCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGc                                                                 >  1:577456/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)

                                                                                                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:1398213/1‑150 (MQ=255)

                                                                                                       cgcgAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGcc                                               >  2:384084/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)

                                                                                                                         cAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAAc                             <  2:1411204/150‑1 (MQ=255)

                                                                                                                         cAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAAc                             >  1:375305/1‑150 (MQ=255)

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

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

                                                                                                                           aaTGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCa                           >  2:514585/1‑150 (MQ=255)

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

                                                                                                                                aaTGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCaaaaaa                      <  1:94993/150‑1 (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)

                                                                                                                                       ggggTATCAACGATGTTGATACCCTGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTg               >  2:1645229/4‑150 (MQ=255)

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

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

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

                                                                                                                                              cAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACgg        <  1:655988/150‑1 (MQ=255)

                                                                                                                                              cAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACgg        >  1:269210/1‑150 (MQ=255)

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

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

                                                                                                                                                aCATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGcc      <  1:514585/150‑1 (MQ=255)

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

                                                                                                                                                    cGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGaa  >  1:642412/1‑150 (MQ=255)

ACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAA  >  NC_000913/4058474‑4058771

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