BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,110,400	A→T	5.5%	noncoding (35/76 nt)	pheV →	tRNA‑Phe

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,110,400	0	A	T	5.5%	374.9 / 8.6	127	noncoding (35/76 nt)	pheV	tRNA‑Phe
Reads supporting (aligned to +/- strand): ref base A (65/55); new base T (4/3); total (69/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

Only 100 of 130 total aligned reads displayed.

TTTGCCTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCAC  >  NC_000913/3110252‑3110548

tttGCCTGACAATGCGTTCAATATCGGCAAAGTGATGATATATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGGTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGaa                                                                                                                                                     <  1:2008313/150‑1 (MQ=255)

 ttGCCTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGaaa                                                                                                                                                    <  1:1892615/150‑1 (MQ=255)

  tGCCTGACAATGCTTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTTACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGTTAGAGCAGGGGATTGaaaa                                                                                                                                                   <  1:1227096/150‑1 (MQ=255)

   gCCTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAAt                                                                                                                                                  >  1:673951/1‑150 (MQ=255)

    ccTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATc                                                                                                                                                 >  1:1187252/1‑150 (MQ=255)

     cTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATcc                                                                                                                                                >  2:828003/1‑150 (MQ=255)

      tGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATccc                                                                                                                                               <  1:2007293/150‑1 (MQ=255)

        aCAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGCTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCg                                                                                                                                             <  1:89437/150‑1 (MQ=255)

         cAATGAGTGAAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTTTCAATcccctg                                                                                                                                            <  1:1515739/150‑3 (MQ=37)

           aTGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTTTCAATcccctgct                                                                                                                                          <  2:1683036/150‑5 (MQ=255)

            tGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTTTCAATcccctgctc                                                                                                                                         <  1:912182/150‑6 (MQ=255)

               gTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCACCGTGTCCtt                                                                                                                                      >  2:197134/1‑150 (MQ=255)

               gTGCAATATCGGCAAACTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGGCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCtt                                                                                                                                      <  2:1730591/150‑1 (MQ=255)

                    atatCTGCAAAGTGATTATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTc                                                                                                                                 <  2:1016219/150‑1 (MQ=255)

                    atatCGGTAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGTTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTTAAAATCCCCGTGTCCTTTGTTc                                                                                                                                 <  2:1560442/150‑1 (MQ=255)

                    atatCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTc                                                                                                                                 >  2:1395169/1‑150 (MQ=255)

                    atatCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTc                                                                                                                                 <  1:1486458/150‑1 (MQ=255)

                        cGGCAAAGTGATGAAAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGAAGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGGAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGAtt                                                                                                                             <  2:1585293/150‑1 (MQ=255)

                         ggCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTc                                                                                                                            >  1:1650755/1‑150 (MQ=255)

                         ggCAAAGTAATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTc                                                                                                                            <  2:418663/150‑1 (MQ=255)

                           cAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCGCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCg                                                                                                                          <  2:1133966/150‑1 (MQ=255)

                            aaaGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGGTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGa                                                                                                                         <  1:829923/150‑1 (MQ=255)

                               gTGATGATCGTTTGTGCAGGCTGCAGTAAATTGAAGAAATATGATTGACGAGACGCGGCGAATCTGGTTTAATGCGCCCCGTTGCCCTGCTAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTc                                                                                                                      <  2:275038/150‑1 (MQ=21)

                               gTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTc                                                                                                                      <  1:1252425/150‑1 (MQ=255)

                               gTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATCGACGAGACGAGGCGAATCAGGTTTACTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAggc                                                                                                                      <  1:301913/150‑3 (MQ=255)

                                  atgatAGATTGTGCAGTCTCCAGTAAATTGAAGAAATTTTATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCgg                                                                                                                   <  2:814343/150‑1 (MQ=255)

                                    gatAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGc                                                                                                                 >  1:629882/1‑150 (MQ=255)

                                    gatAGATTGTGCAGGCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGc                                                                                                                 <  2:1463184/150‑1 (MQ=255)

                                      tAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGcac                                                                                                               >  1:1673822/1‑150 (MQ=255)

                                       aGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGcacc                                                                                                              <  1:868936/150‑1 (MQ=255)

                                       aGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGcacc                                                                                                              <  1:1653227/150‑1 (MQ=255)

                                       aGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGGTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGcacc                                                                                                              <  2:1625490/150‑1 (MQ=255)

                                          ttGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGCCCGGGCACCACt                                                                                                           <  1:1192440/150‑1 (MQ=255)

                                           tgtgCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTa                                                                                                          >  1:1882913/1‑150 (MQ=255)

                                              gCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAAtt                                                                                                       >  1:1596651/1‑150 (MQ=255)

                                              gCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTATAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAAtt                                                                                                       >  2:183163/1‑150 (MQ=255)

                                               ccgTCTGCAGTAAAGTGAAGAAAATTGATTGACCCGACGACGCGACTCAGGTTAAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAACTCCCCGTGTCCTTGGTTCGCTTCCGAGTCCGGGCACCACTAATTc                                                                                                      <  2:1276747/148‑1 (MQ=18)

                                                  tCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTa                                                                                                   <  2:863344/150‑1 (MQ=255)

                                                   cTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTaa                                                                                                  <  1:1280318/150‑1 (MQ=255)

                                                    tGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAg                                                                                                 <  2:1687072/150‑1 (MQ=255)

                                                       aGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAc                                                                                              >  2:1813201/1‑150 (MQ=255)

                                                        gTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAACTCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAcc                                                                                             >  2:1728548/1‑150 (MQ=255)

                                                            aTTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGcc                                                                                         <  1:204420/150‑1 (MQ=255)

                                                              tGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCca                                                                                       <  1:854397/150‑1 (MQ=255)

                                                               gaagaaATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcac                                                                                      >  1:1135881/1‑150 (MQ=255)

                                                                aagaaATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcaca                                                                                     <  2:723004/150‑1 (MQ=255)

                                                                 agaaATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                                                                                    >  2:1904302/1‑150 (MQ=255)

                                                                 agaaATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACACGCCCACaa                                                                                    >  1:61221/1‑150 (MQ=255)

                                                                  gaaATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGACCCGTTGACCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAg                                                                                   >  2:1182858/1‑150 (MQ=255)

                                                                     aTTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCg                                                                                >  1:1141043/1‑150 (MQ=255)

                                                                      tttGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCgg                                                                               >  1:699380/1‑150 (MQ=255)

                                                                       ttgattgaCGAGTCGAGGGGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTGCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCggg                                                                              <  1:1523562/150‑1 (MQ=255)

                                                                        tgattgaCGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                                                                             <  1:1107493/150‑1 (MQ=255)

                                                                        tgattgaCGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                                                                             >  2:726843/1‑150 (MQ=255)

                                                                             gacgagacgcGGCGAATCAGGTTTAATCCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGGCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTg                                                                        <  1:140360/150‑1 (MQ=255)

                                                                             gacgagacgaGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGttttag                                                                        >  1:621313/1‑148 (MQ=255)

                                                                             gacgagacgaGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTg                                                                        >  1:1599178/1‑150 (MQ=255)

                                                                             gacgagacgaGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTg                                                                        >  2:1246785/1‑150 (MQ=255)

                                                                             gacgagacgaGGCGAATCAGGTTTAATGCGCCCCGATGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTg                                                                        >  2:632425/1‑150 (MQ=255)

                                                                             cacgagacgagGCGAATCAGGTGTACTTCGCCCCCTAGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTg                                                                        <  1:308404/149‑1 (MQ=21)

                                                                                 agacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCttt                                                                    <  1:406018/150‑1 (MQ=255)

                                                                                  gacgagGCGAATCAGGGTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtttt                                                                   <  1:1167017/150‑1 (MQ=255)

                                                                                        gCGAATCAGGATTAATGCGCCCCGTTGCCCGGATACCTCATTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt                                                             <  1:267550/150‑1 (MQ=255)

                                                                                                ggTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGGGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGATTTTGCTTTTGGATCTGCCAAGaa                                                     >  1:82023/1‑150 (MQ=255)

                                                                                                  tttAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAAcc                                                   <  2:1756569/150‑1 (MQ=255)

                                                                                                   ttAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCt                                                  <  1:1773055/150‑1 (MQ=255)

                                                                                                     aaTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTc                                                <  1:1770253/150‑1 (MQ=255)

                                                                                                      aTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCa                                               >  1:459603/1‑150 (MQ=255)

                                                                                                      aTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTGTTTGCTTTTGGATCTGACAATAACCTTCa                                               >  1:975125/1‑150 (MQ=255)

                                                                                                      aTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTACAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCa                                               <  2:1852980/150‑1 (MQ=255)

                                                                                                          gcCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGCGTTTTTGCATTTGGATCTGACAATAACCTTCACGaa                                           >  1:246715/1‑150 (MQ=255)

                                                                                                           ccccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCAGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaa                                          >  1:1534900/1‑150 (MQ=255)

                                                                                                            cccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaa                                         >  1:1404404/1‑150 (MQ=255)

                                                                                                                   cccGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAg                                  <  2:1596651/150‑1 (MQ=255)

                                                                                                                    ccGGATAGCTAAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGATCGATTCCGAGTCCGGGCCCCACTAATACTTAAGAACCCGCACACAAGGCGGGTTGTTGCTTTTGGATCTGACAATAACATTCACGAAAAAAATTAgg                                 >  1:1688366/1‑149 (MQ=25)

                                                                                                                      ggATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCtt                               <  1:1607139/150‑1 (MQ=255)

                                                                                                                         tAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTtata                            >  1:409574/1‑150 (MQ=255)

                                                                                                                           gCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATaaa                          >  2:1572850/1‑150 (MQ=255)

                                                                                                                            cTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAg                         >  2:1043381/1‑150 (MQ=255)

                                                                                                                             tCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTTGCTTATAAAGt                        >  1:693186/1‑150 (MQ=255)

                                                                                                                             tCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGt                        >  1:125473/1‑150 (MQ=255)

                                                                                                                              cAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTc                       >  2:1244036/1‑150 (MQ=255)

                                                                                                                                 tCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACGAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTgg                    <  2:1650755/150‑1 (MQ=255)

                                                                                                                                  cGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTggg                   >  1:875354/1‑150 (MQ=255)

                                                                                                                                  cGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTggg                   >  1:882930/1‑150 (MQ=255)

                                                                                                                                  cGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTggg                   >  2:415235/1‑150 (MQ=255)

                                                                                                                                    gTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCGCCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGGTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATACAGTCTggggg                 >  1:1939763/1‑150 (MQ=255)

                                                                                                                                        acacggggATTTTCAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAAtt             >  2:623023/5‑150 (MQ=255)

                                                                                                                                         cacggggATTTTCAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGACCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTa            >  2:570653/4‑150 (MQ=39)

                                                                                                                                          gaGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTAc           <  1:797804/149‑1 (MQ=255)

                                                                                                                                          cAGGGGATTGAAAATCCCCGTGGCCTTTGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTAc           <  1:1644534/150‑1 (MQ=255)

                                                                                                                                           cggggATTTTCAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGCCAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTAct          >  2:1975014/2‑150 (MQ=255)

                                                                                                                                            ggggATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTActc         >  2:1028185/1‑150 (MQ=255)

                                                                                                                                            ggggATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTActa         >  2:1632469/1‑149 (MQ=255)

                                                                                                                                              ggATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAATTCTGGGGGAATTACtctc       >  1:1439876/1‑150 (MQ=255)

                                                                                                                                               gATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCg      <  2:1673822/150‑1 (MQ=255)

                                                                                                                                               gATTGAAAATCACCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCg      <  2:673951/150‑1 (MQ=255)

                                                                                                                                                  tGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa   >  2:552980/1‑150 (MQ=255)

                                                                                                                                                  tGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa   >  2:421706/1‑150 (MQ=255)

                                                                                                                                                   gAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCTAc  >  1:638788/1‑150 (MQ=255)

TTTGCCTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCAC  >  NC_000913/3110252‑3110548

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