BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,928,865	T→G	5.1%	T247P (ACG→CCG)	viaA ←	protein ViaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,928,865	0	T	G	5.1%	474.9 / 12.0	154	T247P (ACG→CCG)	viaA	protein ViaA
Reads supporting (aligned to +/- strand): ref base T (65/81); new base G (5/3); total (70/84)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 156 total aligned reads displayed.

TGCAGGCGATAGGTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACT  >  NC_000913/3928717‑3929014

tGCAGGCGATAGGTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATACCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTc                                                                                                                                                      >  2:754129/1‑150 (MQ=255)

     gCGATAGGTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCgg                                                                                                                                                 <  2:691170/150‑1 (MQ=255)

      cGATAGGTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGt                                                                                                                                                <  2:427599/150‑1 (MQ=255)

        aTAGGTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTc                                                                                                                                              <  1:1090149/150‑1 (MQ=255)

         tAGGTGAGCAACTGTTTTTCCACCAGCCGACGGTAAATCTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTcg                                                                                                                                             <  2:1324029/150‑1 (MQ=255)

         tAGGTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTcg                                                                                                                                             <  1:1201300/150‑1 (MQ=255)

           ggTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTcgcg                                                                                                                                           >  1:1072563/1‑150 (MQ=255)

           ggTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTcgcg                                                                                                                                           >  2:1265593/1‑150 (MQ=255)

           ggTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTcgcg                                                                                                                                           >  1:553830/1‑150 (MQ=255)

           ggTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGGCGACGGTTcctg                                                                                                                                           <  1:996818/150‑4 (MQ=255)

             tGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCa                                                                                                                                         <  1:1172760/150‑1 (MQ=255)

              gAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCa                                                                                                                                         <  1:672087/149‑1 (MQ=255)

              gAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGGCGACGGTTcctgagc                                                                                                                                        <  2:846389/150‑7 (MQ=255)

               aGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCAcc                                                                                                                                       <  1:1212648/150‑1 (MQ=255)

               aGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCAcc                                                                                                                                       <  1:89166/150‑1 (MQ=255)

                        ttttCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCgg                                                                                                                              <  1:223291/150‑1 (MQ=255)

                            ccaccaGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCgcgaac                                                                                                                          <  2:150757/150‑6 (MQ=255)

                            ccaccaGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAgg                                                                                                                          >  1:416714/1‑150 (MQ=255)

                            ccaccaGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAgg                                                                                                                          >  1:1090069/1‑150 (MQ=255)

                            ccaccaGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAgg                                                                                                                          <  1:936085/150‑1 (MQ=255)

                             caccaGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGt                                                                                                                         >  2:1234488/1‑150 (MQ=255)

                              accaGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGtt                                                                                                                        <  1:186402/150‑1 (MQ=255)

                                caGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTc                                                                                                                      <  2:517762/150‑1 (MQ=255)

                                  gCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCa                                                                                                                    <  1:529551/150‑1 (MQ=255)

                                   ccGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCAt                                                                                                                   <  2:1356616/150‑1 (MQ=255)

                                   ccGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCAt                                                                                                                   <  2:1362640/150‑1 (MQ=255)

                                       cGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTgc                                                                                                               >  2:1395474/1‑150 (MQ=255)

                                        ggTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTgcg                                                                                                              <  1:583974/150‑1 (MQ=255)

                                          tAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCa                                                                                                            >  2:76964/1‑150 (MQ=255)

                                                      tCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGt                                                                                                >  1:1083759/1‑150 (MQ=255)

                                                      tCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGt                                                                                                >  1:1176731/1‑150 (MQ=255)

                                                               gTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTg                                                                                       <  2:329236/150‑1 (MQ=255)

                                                                ttATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTgg                                                                                      >  1:145221/1‑150 (MQ=255)

                                                                 tATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGc                                                                                     <  2:854594/150‑1 (MQ=255)

                                                                 aaTCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGc                                                                                     <  2:1315268/149‑1 (MQ=255)

                                                                    cccTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATGAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTc                                                                                  >  1:842841/1‑150 (MQ=255)

                                                                       tAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCg                                                                               >  1:1189156/1‑150 (MQ=255)

                                                                        cctgccGCCAGAACTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCga                                                                              >  2:1238290/6‑150 (MQ=255)

                                                                        aGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTg                                                                                                                                 >  2:539043/1‑99 (MQ=255)

                                                                         gtgtCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGCCCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCgag                                                                             >  2:601533/1‑150 (MQ=255)

                                                                         gtgtCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCgag                                                                             <  2:704018/150‑1 (MQ=255)

                                                                         ctgccGCCAGAACTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCgag                                                                             >  1:368522/5‑150 (MQ=255)

                                                                         ctgccGCCAGAACTGGCGACAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCCGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCgag                                                                             >  1:885329/5‑150 (MQ=255)

                                                                          tgtCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCgaga                                                                            <  2:1189156/150‑1 (MQ=255)

                                                                           gtCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGaa                                                                           >  1:551189/1‑150 (MQ=255)

                                                                             cGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACg                                                                         >  2:708839/1‑150 (MQ=255)

                                                                               ccAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGcc                                                                       <  2:776390/150‑1 (MQ=255)

                                                                               ccAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGcc                                                                       <  1:754129/150‑1 (MQ=255)

                                                                                 aGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGcccc                                                                     <  2:1237205/150‑1 (MQ=255)

                                                                                   ttCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGc                                                                                                                                >  1:155560/1‑89 (MQ=255)

                                                                                   ttCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGc                                                                                                                                <  2:155560/89‑1 (MQ=255)

                                                                                     cTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAgct                                                                 >  2:919805/1‑150 (MQ=255)

                                                                                       ggcggcAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAgctgc                                                               <  1:1395474/150‑1 (MQ=255)

                                                                                        gcggcAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAgctgct                                                              <  1:1265593/150‑1 (MQ=255)

                                                                                         cggcAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGctc                                                             <  1:1050633/150‑1 (MQ=255)

                                                                                           gcAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTg                                                           >  2:965990/1‑150 (MQ=255)

                                                                                             aGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATc                                                                                                                  <  1:67207/93‑1 (MQ=255)

                                                                                             aGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATc                                                                                                                  >  2:67207/1‑93 (MQ=255)

                                                                                             aGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGcc                                                         >  1:617972/1‑150 (MQ=255)

                                                                                              ggAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAg                                                                   >  2:737608/1‑139 (MQ=255)

                                                                                              ggAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCa                                                        <  1:961068/150‑1 (MQ=255)

                                                                                              ggAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCa                                                        <  1:965990/150‑1 (MQ=255)

                                                                                                agaCGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATTTTGCGGAAGGTTTCCATCTGCTCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAgg                                                      <  2:1156982/150‑1 (MQ=255)

                                                                                                   cGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGt                                                   >  2:948846/1‑150 (MQ=255)

                                                                                                     tAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGttt                                                 <  1:1057915/150‑1 (MQ=255)

                                                                                                     tAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGttt                                                 <  2:1071292/150‑1 (MQ=255)

                                                                                                      aaaaTATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTc                                                >  2:293225/1‑150 (MQ=255)

                                                                                                      aaaaTATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTc                                                <  1:728116/150‑1 (MQ=255)

                                                                                                       aaaTATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCa                                               <  1:754896/150‑1 (MQ=255)

                                                                                                        aaTATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAg                                              <  2:1090069/150‑1 (MQ=255)

                                                                                                          tatCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAg                                              >  1:998920/1‑148 (MQ=255)

                                                                                                           atcatcGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTc                                           <  2:416714/150‑1 (MQ=255)

                                                                                                             catcGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGTCGTTTCAGTTCCg                                         >  1:798301/1‑150 (MQ=255)

                                                                                                             catcGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCg                                         >  1:528445/1‑150 (MQ=255)

                                                                                                               tcGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGc                                       <  2:965231/150‑1 (MQ=255)

                                                                                                                  cTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAgctgc                                                               <  2:413997/123‑1 (MQ=255)

                                                                                                                  cTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAgctgc                                                               >  1:413997/1‑123 (MQ=255)

                                                                                                                      gttgCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGt                                >  1:93067/1‑150 (MQ=255)

                                                                                                                        tgCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGtt                                   >  2:338047/1‑145 (MQ=255)

                                                                                                                        tgCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGtt                                   <  1:338047/145‑1 (MQ=255)

                                                                                                                         gCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTaa                             >  1:312409/1‑150 (MQ=255)

                                                                                                                          cAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAg                            <  2:368522/150‑1 (MQ=255)

                                                                                                                          cAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAg                            >  2:871357/1‑150 (MQ=255)

                                                                                                                             aCCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGaaa                         <  2:553830/150‑1 (MQ=255)

                                                                                                                                 tCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCa                     >  1:955862/1‑150 (MQ=255)

                                                                                                                                  cAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCAc                    <  1:130658/150‑1 (MQ=255)

                                                                                                                                   aaCCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCAcc                   <  2:664018/150‑1 (MQ=255)

                                                                                                                                   aaCCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCAcc                   <  2:885329/150‑1 (MQ=255)

                                                                                                                                     ccTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGt                 >  1:443656/1‑150 (MQ=255)

                                                                                                                                       tGCTCAGGAACCGTCGCCTGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGACCGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTAt               <  2:139536/150‑1 (MQ=255)

                                                                                                                                       ggtgcgcgAACCGGCGACGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTATGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTAt               >  2:1049561/8‑150 (MQ=255)

                                                                                                                                        gCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTAtt              >  2:960564/1‑150 (MQ=255)

                                                                                                                                          tCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTc            <  2:1083759/150‑1 (MQ=255)

                                                                                                                                           cgcgAACCGGCGACGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTca           >  1:1253943/4‑150 (MQ=255)

                                                                                                                                           cgcgAACCGGCGACGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTca           >  2:986668/4‑150 (MQ=255)

                                                                                                                                           cAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGTGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTca           >  2:364378/1‑150 (MQ=255)

                                                                                                                                               aaCCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAAt       <  2:506410/150‑1 (MQ=255)

                                                                                                                                                 ccGGCGACGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCa     <  2:1213376/150‑1 (MQ=255)

                                                                                                                                                    tCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACt  <  2:617972/150‑1 (MQ=255)

                                                                                                                                                    tCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACt  <  2:1178139/150‑1 (MQ=255)

TGCAGGCGATAGGTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACT  >  NC_000913/3928717‑3929014

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