BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,537,076	C→G	5.4%	L54L (CTG→CTC)	pdxK ←	pyridoxal kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,537,076	0	C	G	5.4%	370.6 / 8.2	110	L54L (CTG→CTC)	pdxK	pyridoxal kinase I
Reads supporting (aligned to +/- strand): ref base C (50/54); new base G (4/2); total (54/56)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 113 total aligned reads displayed.

TTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAACAACAA  >  NC_000913/2536927‑2537225

ttGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTc                                                                                                                                                       <  1:742595/150‑1 (MQ=255)

ttGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTc                                                                                                                                                       <  1:506651/150‑1 (MQ=255)

   aTTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGc                                                                                                                                                    >  2:773820/1‑150 (MQ=255)

    tttGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCa                                                                                                                                                   >  1:1084517/1‑150 (MQ=255)

        cGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATAc                                                                                                                                               <  1:281322/150‑1 (MQ=255)

          aTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgc                                                                                                                                             <  1:577240/150‑3 (MQ=255)

          aTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgc                                                                                                                                             <  1:202910/150‑3 (MQ=255)

           tGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGt                                                                                                                                            >  1:441214/1‑150 (MQ=255)

             ccGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcg                                                                                                                                          >  2:1537868/1‑150 (MQ=255)

             ccGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcg                                                                                                                                          >  1:350796/1‑150 (MQ=255)

              cGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcgg                                                                                                                                         >  2:1280970/1‑150 (MQ=255)

              cGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcgg                                                                                                                                         <  1:1064246/150‑1 (MQ=255)

              cGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcg                                                                                                                                          >  1:536352/1‑149 (MQ=255)

              cGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcg                                                                                                                                          <  2:536352/149‑1 (MQ=255)

               gTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggc                                                                                                                                        >  1:1010348/1‑150 (MQ=255)

                  cccATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacg                                                                                                                                     >  2:515889/1‑150 (MQ=255)

                    cATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacggc                                                                                                                                   <  2:36325/150‑1 (MQ=255)

                    cATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacggc                                                                                                                                   >  1:1534492/1‑150 (MQ=255)

                        taGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAg                                                                                                                               <  1:967042/150‑1 (MQ=255)

                          gCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGAc                                                                                                                             >  2:406309/1‑150 (MQ=255)

                          gCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGAc                                                                                                                             >  1:379160/1‑150 (MQ=255)

                             cGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACAtt                                                                                                                          >  1:1002086/1‑150 (MQ=255)

                               tGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCa                                                                                                                        <  1:1400605/150‑1 (MQ=255)

                                 gTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAgg                                                                                                                      <  2:864037/150‑1 (MQ=255)

                                   tACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAg                                                                                                                                                     >  1:137859/1‑117 (MQ=255)

                                   tACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAg                                                                                                                                                     <  2:137859/117‑1 (MQ=255)

                                   tACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGcc                                                                                                                    <  1:832634/150‑1 (MQ=255)

                                    aCAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCg                                                                                                                   >  2:1066471/1‑150 (MQ=255)

                                       gCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTc                                                                                                                <  2:1084517/150‑1 (MQ=255)

                                       gCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTc                                                                                                                >  1:883437/1‑150 (MQ=255)

                                        cACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAAACATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCt                                                                                                               >  2:101767/1‑150 (MQ=255)

                                            aaGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACAGCAAAGACATTCAGGCCGTTCTGttt                                                                                                           <  2:350796/150‑1 (MQ=255)

                                              gTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGCCATTCAGGCCGTTCTGTTTGa                                                                                                         >  1:419091/1‑150 (MQ=255)

                                              gTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGa                                                                                                         >  1:801993/1‑150 (MQ=255)

                                               ttGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTgtt                                                                                                        <  2:152440/150‑3 (MQ=255)

                                                 ggCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAg                                                                                                      >  2:75694/1‑150 (MQ=255)

                                                  gcgcagcgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGc                                                                                                     <  2:1534492/150‑1 (MQ=255)

                                                  gcgcagcgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGc                                                                                                     <  2:261698/150‑1 (MQ=255)

                                                   cgcagcgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCa                                                                                                    <  1:1208157/150‑1 (MQ=255)

                                                      agcgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGATGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGc                                                                                                 <  2:1250815/150‑1 (MQ=255)

                                                           aTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGc                                                                                            <  2:971262/150‑1 (MQ=255)

                                                              aCGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAAt                                                                                         >  1:601868/1‑150 (MQ=255)

                                                                   ccTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACg                                                                                              >  1:613772/1‑140 (MQ=255)

                                                                   ccTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACg                                                                                              <  2:613772/140‑1 (MQ=255)

                                                                   ccTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCGGTTTGATAGCAGGCACGGCAATGCTGt                                                                                    >  2:273361/1‑150 (MQ=255)

                                                                    cTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGtt                                                                                   >  1:285611/1‑150 (MQ=255)

                                                                     tGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTg                                                                                  >  1:781089/1‑150 (MQ=255)

                                                                     tGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTg                                                                                  <  2:441214/150‑1 (MQ=255)

                                                                        aGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGccc                                                                               <  2:695044/150‑1 (MQ=255)

                                                                         gcgcACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCa                                                                              <  2:150053/150‑1 (MQ=255)

                                                                         gcgcACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCa                                                                              >  1:806606/1‑150 (MQ=255)

                                                                         gcgcACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCa                                                                              <  1:773820/150‑1 (MQ=255)

                                                                             acgcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCt                                                                          <  1:768891/150‑1 (MQ=255)

                                                                              cgcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTg                                                                         >  2:149178/1‑150 (MQ=255)

                                                                                caAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGcc                                                                       <  2:508037/150‑1 (MQ=255)

                                                                                 aaaTAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCg                                                                      <  1:1007639/150‑1 (MQ=255)

                                                                                 aaaTAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCg                                                                      <  1:81746/150‑1 (MQ=255)

                                                                                      gCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAac                                                                 <  1:75694/150‑1 (MQ=255)

                                                                                       ccGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAacc                                                                >  1:1187624/1‑150 (MQ=255)

                                                                                       ccGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAacc                                                                <  2:582283/150‑1 (MQ=255)

                                                                                        cGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATTATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAacca                                                               <  2:840929/150‑1 (MQ=255)

                                                                                          cTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGAGATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAaccacc                                                             <  2:1503089/150‑1 (MQ=255)

                                                                                          cTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAaccacc                                                             <  1:406309/150‑1 (MQ=255)

                                                                                           tAAACCATTCGTCCGGAATCGCACCACCGGAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCt                                                            >  2:971476/1‑150 (MQ=255)

                                                                                              aCCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGcc                                                         >  2:1050671/1‑149 (MQ=255)

                                                                                                 aTTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACt                                                      >  1:256902/1‑150 (MQ=255)

                                                                                                    cGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCa                                                   <  1:59993/150‑1 (MQ=255)

                                                                                                     gTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCAc                                                  >  2:945293/1‑150 (MQ=255)

                                                                                                     gTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCAc                                                  <  1:1280970/150‑1 (MQ=255)

                                                                                                      tCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACg                                                 >  1:562021/1‑150 (MQ=255)

                                                                                                         ggAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGcg                                              <  1:1074018/150‑1 (MQ=255)

                                                                                                         ggAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGcg                                              <  1:200839/150‑1 (MQ=255)

                                                                                                            aTCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGcgacg                                           >  2:232555/1‑150 (MQ=255)

                                                                                                             tCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGcgacga                                          <  1:881525/150‑1 (MQ=255)

                                                                                                              cGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGat                                         <  1:1066471/150‑1 (MQ=255)

                                                                                                                       gTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTg                                >  1:74239/1‑150 (MQ=255)

                                                                                                                        tAGAAAGTGTCATAATGCGTCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGTCGACGATATCCGCCTGc                               <  2:1117639/150‑1 (MQ=255)

                                                                                                                        tAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCGCGGCGACGATATCCGCCTGc                               >  1:351676/1‑150 (MQ=255)

                                                                                                                          gAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAg                             >  2:1025229/1‑150 (MQ=255)

                                                                                                                             aGTGTCATAATGCGGCGTATTGCTCAGCAATACCGGCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCGGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCGGCAGTGc                          <  2:192430/150‑1 (MQ=255)

                                                                                                                                gtCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATGGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCt                       >  1:1431734/1‑150 (MQ=255)

                                                                                                                                gtCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCt                       >  2:1214756/1‑150 (MQ=255)

                                                                                                                                 tCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGGTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGAGACGATATCCGCCTGCAGTGCCCTa                      <  2:427473/150‑1 (MQ=255)

                                                                                                                                  cATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAc                     <  2:1010348/150‑1 (MQ=255)

                                                                                                                                    tAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTActc                   >  1:1307407/1‑150 (MQ=255)

                                                                                                                                    tAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTActc                   <  2:562021/150‑1 (MQ=255)

                                                                                                                                       tGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTa                <  1:101767/150‑1 (MQ=255)

                                                                                                                                        gcggcgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTAt               <  2:683444/150‑1 (MQ=255)

                                                                                                                                         gacggTATTGCTGAGCAATACCGTCGGCACGTCAAAGACAATCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATc              >  2:1262511/5‑150 (MQ=255)

                                                                                                                                         gacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATc              >  1:767760/5‑150 (MQ=255)

                                                                                                                                         gacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATc              >  1:772975/5‑150 (MQ=255)

                                                                                                                                          ggcgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCg             <  1:1356172/150‑1 (MQ=255)

                                                                                                                                          ggcgTATTGCTCAGCAATACCGGCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCg             <  1:827660/150‑1 (MQ=255)

                                                                                                                                             gTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTa          <  1:1537868/150‑1 (MQ=255)

                                                                                                                                                ttGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAaac       <  1:536200/150‑1 (MQ=255)

                                                                                                                                                 tGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCAGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAaaca      <  2:443269/150‑1 (MQ=255)

                                                                                                                                                  gCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAaacaa     >  1:257283/1‑150 (MQ=255)

                                                                                                                                                   cTCAGCAATACCGTCGGCACGGCAATGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCGCGCTGCCGTAAACCCCCTGCGACTGCACGGCGAGGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAgacccc    >  2:1458789/1‑147 (MQ=255)

                                                                                                                                                    tCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAaacaaca   <  2:1300858/150‑1 (MQ=255)

                                                                                                                                                     cAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAaacaacaa  >  1:598588/1‑150 (MQ=255)

TTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAACAACAA  >  NC_000913/2536927‑2537225

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