BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,058,620	A→G	5.9%	I72V (ATC→GTC) ‡	typA →	GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,058,620	0	A	G	5.9%	394.1 / 11.7	135	I72V (ATC→GTC) ‡	typA	GTP‑binding protein
Reads supporting (aligned to +/- strand): ref base A (65/62); new base G (3/5); total (68/67)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.79e-01

Only 100 of 138 total aligned reads displayed.

TAGACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTG  >  NC_000913/4058471‑4058769

tAGACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACa                                                                                                                                                       >  1:544404/1‑150 (MQ=255)

tAGACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACa                                                                                                                                                       >  1:1608870/1‑150 (MQ=255)

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

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

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

      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)

            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)

                       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)

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

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

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

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

                                                 ccAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACgtgt                                                                                                      <  2:72410/150‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                  gaAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGTATGATTACCGTATCAACATCGTTGATACCCCGGGTATCGCCGACTTCGGTGGTGAAGTTGACCGTGTAATGTCCATGGTAGACTCAGTGCTGCTggtggt                                                                     >  2:1039291/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: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)

                                                                                                                     ccGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCt                                  >  1:638106/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)

TAGACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCCATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTG  >  NC_000913/4058471‑4058769

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