BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,110,401	A→C	6.0%	noncoding (36/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,401	0	A	C	6.0%	497.9 / 15.2	161	noncoding (36/76 nt)	pheV	tRNA‑Phe
Reads supporting (aligned to +/- strand): ref base A (79/71); new base C (3/7); total (82/79)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

TGCCTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGT  >  NC_000913/3110254‑3110550

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

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

   cTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATcc                                                                                                                                                  <  2:95523/150‑1 (MQ=255)

   cTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATcc                                                                                                                                                  <  2:802827/150‑1 (MQ=255)

    tGACAATGCGTGCAATATCGGCAAAGTGTTGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCTCCGCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATccc                                                                                                                                                 <  1:1616385/150‑1 (MQ=255)

    tGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATccc                                                                                                                                                 >  2:1068342/1‑150 (MQ=255)

     taCAATTCGTGCAATATCGGCAAAGTGATGATAGATTGTGAAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGACGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATcccc                                                                                                                                                <  1:942300/149‑1 (MQ=255)

     gACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATcccc                                                                                                                                                >  1:2112873/1‑150 (MQ=255)

      aCAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTTTCAATcccct                                                                                                                                               <  1:950290/150‑2 (MQ=255)

      aCAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTTTCAATcccct                                                                                                                                               <  1:2025572/150‑2 (MQ=255)

      aCAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCg                                                                                                                                               <  1:1920400/150‑1 (MQ=255)

       cAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTTTCAATcccctg                                                                                                                                              <  1:1722017/150‑3 (MQ=255)

         aTGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTTTCAATcccctgct                                                                                                                                            <  1:1214889/150‑5 (MQ=255)

         aTGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCgtgt                                                                                                                                            <  1:479637/150‑1 (MQ=255)

         aTGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCgtgt                                                                                                                                            >  2:1480744/1‑150 (MQ=255)

         aTGCGTGCAATATCGGCAAAGGGTTGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTTTCAATcccctgct                                                                                                                                            <  2:2065419/150‑5 (MQ=37)

               gCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTgg                                                                                                                                      <  1:1510644/150‑1 (MQ=255)

               gCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTgg                                                                                                                                      <  2:1090836/150‑1 (MQ=255)

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

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

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

                       ggCAAAGTGATGATAGATGGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTc                                                                                                                              <  1:1693402/150‑1 (MQ=255)

                         cAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCg                                                                                                                            <  2:790285/150‑1 (MQ=255)

                           aaGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAg                                                                                                                          <  1:1183758/150‑1 (MQ=255)

                                 tgatAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGAAAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCggg                                                                                                                    >  2:95063/1‑150 (MQ=255)

                                  tatAGATTGTGCAGTCTGCTGTAAATTGAAGAATTGTGATTGACGAGCCGAGGCGACTCAGGTTTAATGCGCCGCGTTGCCTGGATAGCTCAGTCGGTAGAGCAGGGGATTGAATATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGc                                                                                                                   <  1:2148292/149‑1 (MQ=18)

                                   atAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGAAGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGTATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCACGTGTCCTTGGTTCGATTCCGAGTCAGGGca                                                                                                                  >  1:46253/1‑150 (MQ=255)

                                          gtgCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTaa                                                                                                           >  2:439149/1‑150 (MQ=255)

                                           tgCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATTCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGTTTCGATTCCGAGTCCGGGCACCACTAAt                                                                                                          <  2:37252/150‑1 (MQ=255)

                                              aGTCTGCAGTATATTGAAGAAATTTGTTTGACGAGACGAGGGGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTG‑AAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCtt                                                                                                      <  2:2023602/150‑1 (MQ=255)

                                              aGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCt                                                                                                       <  1:973541/150‑1 (MQ=255)

                                              aGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTATAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCt                                                                                                       >  2:113160/1‑150 (MQ=255)

                                              aGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAACCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCt                                                                                                       <  1:847243/150‑1 (MQ=255)

                                               gTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCtt                                                                                                      <  2:2112873/150‑1 (MQ=255)

                                               gTCTGCAGTAAATTGAAGAAATATGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCAGGGCACCACTAATTCtt                                                                                                      >  1:1372139/1‑150 (MQ=255)

                                                 cTGCAGTAAATTGACGAACTTTGATTGACGACACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCACAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTGCGATTCCGAGTCCGGGCACCACTAATTCTTaa                                                                                                    <  1:1144670/150‑1 (MQ=21)

                                                 cTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGCTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCATaa                                                                                                    <  2:1126578/150‑1 (MQ=255)

                                                    cAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGaa                                                                                                 >  2:1270552/1‑150 (MQ=255)

                                                      gTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAcc                                                                                               >  1:641382/1‑150 (MQ=255)

                                                        aaaTTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAcccc                                                                                             >  2:2040555/1‑149 (MQ=255)

                                                            tGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCca                                                                                         >  2:625670/1‑150 (MQ=255)

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

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

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

                                                                      tgattgaCGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCAACACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                                                                               >  1:2090646/1‑150 (MQ=255)

                                                                       gattgaCGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtt                                                                              >  2:1989918/1‑150 (MQ=255)

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

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

                                                                            acgagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGc                                                                         >  2:973498/1‑150 (MQ=255)

                                                                             cgagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCt                                                                        <  1:1448239/150‑1 (MQ=255)

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

                                                                                  cgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCAACAAGTCGGGAGTTTGGTTTTgg                                                                   >  2:320020/1‑150 (MQ=18)

                                                                                   gagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGa                                                                  >  1:1893904/1‑150 (MQ=255)

                                                                                   gagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGa                                                                  <  2:1504900/150‑1 (MQ=255)

                                                                                     ggCGACTCAGGTTTAATGCGGCCCGGTGCCCGGATACCTCAGTCGGTAGAGCAGGGGATTGACAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                                                                <  1:1677209/150‑1 (MQ=21)

                                                                                     ggCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                                                                >  1:1636670/1‑150 (MQ=255)

                                                                                     ggCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                                                                >  1:577278/1‑150 (MQ=255)

                                                                                     ggCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                                                                <  1:1703065/150‑1 (MQ=255)

                                                                                     ggCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                                                                >  2:2072491/1‑150 (MQ=255)

                                                                                      gcttaTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCAGAGTCCGGTCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt                                                               <  2:377680/146‑1 (MQ=255)

                                                                                      gCGATTCATGTTTAATGCGCCCCGTTGCCAGGATAGATCATTCGGTAGAGCAGGGGATTGAAAATCCACGTGTCCTTTGTTCGATTCCTAGTCCTGGCACCACTAATTCTTAAGAACCCGCCCACAAGTCGGGTTTTTGCTTTTGGATCt                                                               <  2:1942371/150‑1 (MQ=21)

                                                                                      gCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt                                                               <  2:864247/150‑1 (MQ=255)

                                                                                       ggAATCAGGTTTAATGCTCCCCGTTGCGCGTATAGCTCAGTCGGTAGAGCAGGGGTTTGAAAATCCCCGTGTCGTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTg                                                              <  2:649725/149‑1 (MQ=21)

                                                                                       cGAATCAGGTTTACTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCATGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTg                                                              <  1:2124700/150‑1 (MQ=255)

                                                                                         acgcAGGTCTAATGCGCCCCGTTGCCCGGATCGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGATCGATGCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGAc                                                            <  2:1658668/147‑1 (MQ=21)

                                                                                          aTCAGGTTTAATGTGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACa                                                           <  2:6112/150‑1 (MQ=255)

                                                                                          aTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACa                                                           >  1:533616/1‑150 (MQ=255)

                                                                                              ggTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCATTTGGATCTGAAAATaa                                                       >  2:914735/1‑150 (MQ=255)

                                                                                               gTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAAc                                                      >  2:1201829/1‑150 (MQ=255)

                                                                                                  tAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTATATTCCGAGTACGGGCACCACTAAATCTTAACAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCtt                                                   >  1:594294/1‑150 (MQ=21)

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

                                                                                                     tGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCAc                                                <  2:1814511/150‑1 (MQ=255)

                                                                                                     tGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCAc                                                <  1:1128553/150‑1 (MQ=255)

                                                                                                     tGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCAc                                                <  1:439149/150‑1 (MQ=255)

                                                                                                        gcCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaa                                             >  1:1904327/1‑150 (MQ=255)

                                                                                                          cccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACACGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGAAAATAACCTTCAAGaaaa                                           >  2:231790/1‑150 (MQ=255)

                                                                                                               tGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTTATAATCCCCGTATCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGATCCCGCCCACACGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAAtt                                      <  2:1465104/150‑1 (MQ=25)

                                                                                                               tGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAAtt                                      >  2:1436471/1‑150 (MQ=255)

                                                                                                                  ccGGATTGATCAGTCGTTAGAGCAGTGGTTTGAAAATCCCCGTGTCCTTGTTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCAACAAGGCGGGTTTTTGCTTTTGGATCTTACAATAACCTTCACGAAAAAAATTAGc                                   <  1:1270708/150‑1 (MQ=37)

                                                                                                                    ggATAGGACAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCtt                                 <  1:1002225/150‑1 (MQ=255)

                                                                                                                     gATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTta                                <  2:1098116/150‑1 (MQ=255)

                                                                                                                         gCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTATGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGTTCTGACAATAACCTTCACGAAACAAATTAGCTTATaaa                            >  1:1946270/1‑150 (MQ=255)

                                                                                                                         gCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGCTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATaaa                            >  2:127590/1‑150 (MQ=255)

                                                                                                                          cTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAg                           <  1:1843079/150‑1 (MQ=255)

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

                                                                                                                                 ggTAGAGCAGGTGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGCCTgggg                    <  2:429598/150‑1 (MQ=255)

                                                                                                                                 ggTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTgggg                    <  1:1033070/150‑1 (MQ=255)

                                                                                                                                 ggTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTgggg                    >  1:857516/1‑150 (MQ=255)

                                                                                                                                  gTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCGGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTggggg                   <  2:1014860/150‑1 (MQ=255)

                                                                                                                                  gTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTggggg                   <  1:1349475/150‑1 (MQ=255)

                                                                                                                                     gacacggggATTTTCAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAAt                >  2:683299/6‑150 (MQ=255)

                                                                                                                                       gCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTa              <  2:1944406/150‑1 (MQ=255)

                                                                                                                                       cacggggATTTTCAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTa              >  2:2059885/4‑150 (MQ=255)

                                                                                                                                            ggATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAATAATTAGCTTATAAAGTCTGGGGGAATTACtctc         >  2:1175570/1‑150 (MQ=255)

                                                                                                                                              aTTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc       >  1:684555/1‑150 (MQ=255)

                                                                                                                                               ttGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCCCCACTAATTCTTAAGAACCCGCCCACAAGGCGGGGTTTTGCTTTTGGAGCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGTGAATTACTCTCgac      >  1:1129739/1‑148 (MQ=255)

                                                                                                                                                tGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGGCCGTGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa     <  1:148868/150‑1 (MQ=255)

                                                                                                                                                 gAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCAc    >  1:906507/1‑150 (MQ=255)

                                                                                                                                                 gAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCAc    >  2:560201/1‑150 (MQ=255)

                                                                                                                                                  aaaaTCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACg   >  1:710310/1‑150 (MQ=255)

TGCCTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGT  >  NC_000913/3110254‑3110550

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