BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,090,226	T→C	5.2%	C158R (TGC→CGC)	yggI →	conserved protein YggI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,090,226	0	T	C	5.2%	545.9 / 12.9	173	C158R (TGC→CGC)	yggI	conserved protein YggI
Reads supporting (aligned to +/- strand): ref base T (79/85); new base C (5/4); total (84/89)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

Only 100 of 175 total aligned reads displayed.

GTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGT  >  NC_000913/3090077‑3090375

  gtTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGGTTACCGTTCGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGc                                                                                                                                                     <  1:1729871/150‑1 (MQ=255)

  gtTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGc                                                                                                                                                     <  2:2077393/150‑1 (MQ=255)

  gtTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTCTGTTCACTGc                                                                                                                                                     >  2:940090/1‑150 (MQ=255)

    tgccgcccgTCGGACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGAATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCgg                                                                                                                                                   <  2:1559884/148‑1 (MQ=255)

    tCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCgg                                                                                                                                                   <  1:925611/150‑1 (MQ=255)

     cccgcccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCCACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCCGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGACTATCGCTGTGTTAACTGCGGt                                                                                                                                                  >  1:674601/1‑150 (MQ=255)

     cccgcccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGt                                                                                                                                                  <  1:2404906/150‑1 (MQ=255)

     cccgcccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGt                                                                                                                                                  >  2:1287862/1‑150 (MQ=255)

     cccgcccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGt                                                                                                                                                  >  1:1724262/1‑150 (MQ=255)

       cgcccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGa                                                                                                                                                <  2:1235204/150‑1 (MQ=255)

          ccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACa                                                                                                                                             >  2:158134/1‑150 (MQ=255)

           cgtcgtACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACCGCAGTGAACac                                                                                                                                            <  1:107144/150‑2 (MQ=255)

             tcgtACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCt                                                                                                                                          <  2:253459/150‑1 (MQ=255)

              cgtACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACCGCAGTGAACacagc                                                                                                                                         <  2:2211929/150‑5 (MQ=255)

              cgtACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCACCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACCGCAGTGAACacagc                                                                                                                                         <  2:1169485/150‑5 (MQ=255)

               gtACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTgg                                                                                                                                        >  2:851275/1‑150 (MQ=255)

                  cGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCCTCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTg                                                                                                                                     >  2:2354928/1‑150 (MQ=255)

                    cATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCg                                                                                                                                   >  1:687843/1‑150 (MQ=255)

                    cATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCg                                                                                                                                   <  1:1101474/150‑1 (MQ=255)

                       cAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGaaa                                                                                                                                <  1:508280/150‑1 (MQ=255)

                        aGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAAt                                                                                                                               >  1:490793/1‑150 (MQ=255)

                         gTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATa                                                                                                                              <  1:1546956/150‑1 (MQ=255)

                         gTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCACAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATa                                                                                                                              <  2:653967/150‑1 (MQ=255)

                           tCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAAc                                                                                                                            >  2:1814605/1‑150 (MQ=255)

                           tCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAAc                                                                                                                            >  2:2471359/1‑150 (MQ=255)

                            cGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAAcc                                                                                                                           >  2:829725/1‑150 (MQ=255)

                              aaCTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCAt                                                                                                                         >  2:1535060/1‑150 (MQ=255)

                               aCTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATc                                                                                                                        <  2:1516433/150‑1 (MQ=255)

                                  gCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGa                                                                                                                     <  1:211975/150‑1 (MQ=255)

                                  gCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGa                                                                                                                     >  2:652092/1‑150 (MQ=255)

                                  gCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGa                                                                                                                     >  1:1830606/1‑150 (MQ=255)

                                  gCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGa                                                                                                                     >  1:1839319/1‑150 (MQ=255)

                                  gCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGCAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGa                                                                                                                     >  1:995398/1‑150 (MQ=255)

                                       ccGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTAt                                                                                                                >  1:478591/1‑150 (MQ=255)

                                        cGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATc                                                                                                               <  2:758810/150‑1 (MQ=255)

                                              tCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAAc                                                                                                         <  2:452519/150‑1 (MQ=255)

                                                gCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACtt                                                                                                       >  1:724760/1‑150 (MQ=255)

                                                   acacCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTcc                                                                                                    <  1:1894756/150‑1 (MQ=255)

                                                     acCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTg                                                                                                  <  2:808546/150‑1 (MQ=255)

                                                      ccTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGa                                                                                                 >  1:214935/1‑150 (MQ=255)

                                                      ccTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGa                                                                                                 >  2:894161/1‑150 (MQ=255)

                                                      ccTTCCACTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGa                                                                                                 >  2:1500232/1‑150 (MQ=255)

                                                            ccTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGc                                                                                           <  2:1722867/150‑1 (MQ=255)

                                                            ccTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGc                                                                                           <  2:481983/150‑1 (MQ=255)

                                                             cTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCt                                                                                          <  2:324962/150‑1 (MQ=255)

                                                             cTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCt                                                                                          <  2:1214767/150‑1 (MQ=255)

                                                             cTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCt                                                                                          <  2:1136652/150‑1 (MQ=255)

                                                              tACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGACCTTTCCTGATCTGGCTg                                                                                         >  1:570605/1‑150 (MQ=255)

                                                               aCCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGa                                                                                        <  1:2214658/150‑1 (MQ=255)

                                                                  gCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGGACTTTCCTGATCTGGCTGATTg                                                                                     <  2:1152278/150‑1 (MQ=255)

                                                                   cTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGc                                                                                    <  1:1287862/150‑1 (MQ=255)

                                                                     gCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCAt                                                                                  >  2:1451029/1‑150 (MQ=255)

                                                                       aaGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATAc                                                                                >  1:2022565/1‑150 (MQ=255)

                                                                       aaGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATAc                                                                                <  1:581928/150‑1 (MQ=255)

                                                                         gTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCa                                                                              <  1:1794834/150‑1 (MQ=255)

                                                                         gTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCa                                                                              <  2:242734/150‑1 (MQ=255)

                                                                          tGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaa                                                                             <  1:1535060/150‑1 (MQ=255)

                                                                          tGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaa                                                                             >  1:2469737/1‑150 (MQ=255)

                                                                          tGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaa                                                                             <  2:1947902/150‑1 (MQ=255)

                                                                           gCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTGCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaaa                                                                            <  2:1101554/150‑1 (MQ=255)

                                                                            ccAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaaaa                                                                           >  2:1672809/1‑150 (MQ=255)

                                                                             cAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAAc                                                                          <  1:862785/150‑1 (MQ=255)

                                                                                gAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGc                                                                       <  2:1115986/150‑1 (MQ=255)

                                                                                        gCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTa                                                               <  2:1829310/150‑1 (MQ=255)

                                                                                        gCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTa                                                               <  1:1410639/150‑1 (MQ=255)

                                                                                         cTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTAc                                                              <  1:731734/150‑1 (MQ=255)

                                                                                         cTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTAc                                                              <  2:2007216/150‑1 (MQ=255)

                                                                                             ccGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg                                                          <  1:266049/150‑1 (MQ=255)

                                                                                                 aCGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCATGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCTCTACTTTTCttg                                                      <  2:1253068/150‑3 (MQ=255)

                                                                                                  cgccgcCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc                                                     >  2:91308/1‑150 (MQ=255)

                                                                                                  cgccgcCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTACATACCAAAACAGCTTTCGCTACGTTGCTGGc                                                     >  1:1048807/1‑150 (MQ=255)

                                                                                                       ccATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGtt                                                <  2:2399320/150‑1 (MQ=255)

                                                                                                        cATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGttt                                               <  1:1597401/150‑1 (MQ=255)

                                                                                                                cgTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACgg                                       <  2:899270/150‑1 (MQ=255)

                                                                                                                cgTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACgg                                       <  1:851275/150‑1 (MQ=255)

                                                                                                                cgTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACgg                                       <  1:467756/150‑1 (MQ=255)

                                                                                                                 gTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGa                                      >  2:96985/1‑150 (MQ=255)

                                                                                                                       cGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGt                                <  2:1428959/150‑1 (MQ=255)

                                                                                                                       cGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGt                                >  1:177463/1‑150 (MQ=255)

                                                                                                                        gTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTg                               >  1:1489903/1‑150 (MQ=255)

                                                                                                                          ggCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGAt                             <  2:497091/150‑1 (MQ=255)

                                                                                                                          ggCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGAt                             >  2:448683/1‑150 (MQ=255)

                                                                                                                            cGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTCGCTCGTTTTAACACGGAGTAAGTCATGt                           >  2:590518/1‑150 (MQ=255)

                                                                                                                             gAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTa                          >  1:1626856/1‑150 (MQ=255)

                                                                                                                                 ccGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGTTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGt                      <  1:1141149/150‑1 (MQ=255)

                                                                                                                                  cGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGtt                     >  1:2347138/1‑150 (MQ=255)

                                                                                                                                       aTCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACATGGAGTAAGTGATGTACCGTTATTTg                >  1:2313196/1‑150 (MQ=255)

                                                                                                                                       aTCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTg                <  1:91308/150‑1 (MQ=255)

                                                                                                                                          gCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTa                                   >  2:1157547/1‑128 (MQ=255)

                                                                                                                                          gCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCt             >  1:2095816/1‑150 (MQ=255)

                                                                                                                                          gCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCt             >  1:342105/1‑150 (MQ=255)

                                                                                                                                          gCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCt             >  2:1544571/1‑150 (MQ=255)

                                                                                                                                            tgtgTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTAt           >  2:1887480/1‑150 (MQ=255)

                                                                                                                                            gctgtTCACCGCAGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTAt           >  1:1395323/3‑150 (MQ=255)

                                                                                                                                            gctgtTCAACGCAGTGAAAAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTAt           >  2:230420/3‑150 (MQ=255)

                                                                                                                                                ttCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATtgct       <  2:1063051/150‑1 (MQ=255)

                                                                                                                                                  cACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc                                                   <  2:276261/104‑1 (MQ=255)

                                                                                                                                                  cACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc                                                   >  1:276261/1‑104 (MQ=255)

                                                                                                                                                  cACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATtgctgc     <  1:409074/150‑1 (MQ=255)

                                                                                                                                                     tGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCggt  <  2:724760/150‑1 (MQ=255)

GTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGT  >  NC_000913/3090077‑3090375

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