BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,821,489	T→A	6.9%	I21N (ATT→AAT)	gmk →	guanylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,821,489	0	T	A	6.9%	406.5 / 14.7	145	I21N (ATT→AAT)	gmk	guanylate kinase
Reads supporting (aligned to +/- strand): ref base T (67/67); new base A (4/6); total (72/73)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

Only 100 of 146 total aligned reads displayed.

GTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCG  >  NC_000913/3821340‑3821637

gtgAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCTTATGTAGGCTTTATTTCGCTAATCACATACGTAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGACTGAt                                                                                                                                                      <  2:159181/150‑1 (MQ=255)

  gAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATc                                                                                                                                                    <  1:936552/150‑1 (MQ=255)

   aGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCa                                                                                                                                                   >  1:10010/1‑150 (MQ=255)

      aaGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGc                                                                                                                                                <  2:1463800/150‑1 (MQ=255)

         ggCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCttt                                                                                                                                             <  2:577535/150‑1 (MQ=255)

           cTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTttat                                                                                                                                           <  1:1145683/150‑1 (MQ=255)

           cTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCAGGCtggat                                                                                                                                           <  2:78315/150‑5 (MQ=255)

           cTACGTCACATGGCAGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCTGGCtggat                                                                                                                                           <  2:686979/150‑5 (MQ=255)

            tACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCAGGCtggatt                                                                                                                                          <  2:551583/150‑6 (MQ=255)

              cGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaa                                                                                                                                        <  1:647215/150‑1 (MQ=255)

              cGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaa                                                                                                                                        >  2:1423420/1‑150 (MQ=255)

               gTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaaa                                                                                                                                       >  1:2372/1‑150 (MQ=255)

                  acaTGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAAc                                                                                                                                    <  2:68549/150‑1 (MQ=255)

                  acaTGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAAc                                                                                                                                    >  1:1719670/1‑150 (MQ=255)

                          gcgcCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGt                                                                                                                            <  1:1157919/150‑1 (MQ=255)

                            gcCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTg                                                                                                                          <  2:1571965/150‑1 (MQ=255)

                               gtgtATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTAt                                                                                                                       <  2:1607536/150‑1 (MQ=255)

                               gtgtATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTAt                                                                                                                       <  1:796741/150‑1 (MQ=255)

                               gtgtATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTAt                                                                                                                       >  1:1131326/1‑150 (MQ=255)

                                 ttataATAAGCTCGTATGTAGGCTTTATTTCTCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGa                                                                                                                     <  2:433629/149‑1 (MQ=255)

                                  tataATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCATGATTAAGGCTTTATTAAAAACCCAACTGTTGTATGac                                                                                                                    >  1:758477/1‑150 (MQ=255)

                                     aataaGCTCGTATGTAGGCTTTATTTCGCTAATAACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTTCACCCAGTGGAGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAAAAGTTGTATCACAcc                                                                                                                 >  1:1260708/1‑150 (MQ=255)

                                           cTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGtt                                                                                                           >  1:858255/1‑150 (MQ=255)

                                           cTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGtt                                                                                                           >  1:679947/1‑150 (MQ=255)

                                           cTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGtt                                                                                                           >  1:1081345/1‑150 (MQ=255)

                                           cTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTATATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGAAGTATGACACCCAGGtt                                                                                                           >  1:1462779/1‑150 (MQ=255)

                                              gTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCt                                                                                                        <  1:1576140/150‑1 (MQ=255)

                                              gTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTCGTATGACACCCAGGTTTCt                                                                                                        <  2:1567953/150‑1 (MQ=255)

                                              gTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTCGTATGACACCCAGGTTTCt                                                                                                        <  2:1560993/150‑1 (MQ=255)

                                               tATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTg                                                                                                       <  1:546339/150‑1 (MQ=255)

                                                  gTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGttt                                                                                                    <  1:983612/150‑1 (MQ=255)

                                                   tAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTc                                                                                                   <  2:1162898/150‑1 (MQ=255)

                                                     ggCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcac                                                                                                 <  2:520102/150‑1 (MQ=255)

                                                     ggCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCAAACCGTTGTATGACACCCAGGTTTCTGTTTcac                                                                                                 >  1:783088/1‑150 (MQ=255)

                                                      gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTTTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCAGGCATTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTcaca                                                                                                >  1:1042009/1‑150 (MQ=255)

                                                      gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcaca                                                                                                >  1:1616019/1‑150 (MQ=255)

                                                       cTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcacac                                                                                               <  1:1251665/150‑1 (MQ=255)

                                                          tATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCAcacacc                                                                                            <  2:784890/150‑1 (MQ=255)

                                                           aTTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCAcacacca                                                                                           <  2:372796/150‑1 (MQ=255)

                                                                gCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTTTATGACACCCAGGTTTCTGTTTCACACACCACACGc                                                                                      >  1:333970/1‑150 (MQ=255)

                                                                    aTCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAAc                                                                                  >  2:1181376/1‑150 (MQ=255)

                                                                     tCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAAcc                                                                                 <  1:66052/150‑1 (MQ=255)

                                                                      cacaTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACcg                                                                                <  2:208128/150‑1 (MQ=255)

                                                                       acaTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACcgc                                                                               <  1:447263/150‑1 (MQ=255)

                                                                       acaTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACcgc                                                                               >  2:38740/1‑150 (MQ=255)

                                                                        caTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACcgcg                                                                              <  1:576431/150‑1 (MQ=255)

                                                                           aCGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTcc                                                                           <  1:346764/150‑1 (MQ=255)

                                                                             gAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACCCCACACGCCAACCGCGTCCTg                                                                         >  1:1074580/1‑150 (MQ=255)

                                                                               aaGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGt                                                                       >  2:160349/1‑150 (MQ=255)

                                                                                  aTACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGaa                                                                    >  2:641108/1‑150 (MQ=255)

                                                                                    aCTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGt                                                                  <  1:401659/150‑1 (MQ=255)

                                                                                     cTCATGGCTCAAGGCACGCTTTTTATTGTTTCTGCCCCCAGTGGCGCGGGTAAAACCAGCCTGATTCATCCTTTATTAAAACCCCAACCGTTGTATGCCACCCCGGTGTCGGTGTCACACCCCACCAGCCAACCGCGTCCTGGTGAAGTc                                                                 >  2:859167/1‑150 (MQ=255)

                                                                                     cTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGAGAAGTc                                                                 >  1:758030/1‑150 (MQ=255)

                                                                                       cATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAGTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACTCCAACCGCGTCCTGGTGAAGTccc                                                               >  1:1429842/1‑149 (MQ=255)

                                                                                         tGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCAGATTCAGGCTTTATAAAAAACCCAACCGTTGTATGAAACCCAGGTTTCTGTTTCACACCCCCCACGCCAACCGCGTCCTGGTGAAGTCCACg                                                             >  2:253482/1‑150 (MQ=255)

                                                                                          ggCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCCCGTCCTGGTGAAGTCCACgg                                                            >  2:856349/1‑150 (MQ=255)

                                                                                           gCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGt                                                           >  1:14137/1‑150 (MQ=255)

                                                                                                aGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACa                                                      >  1:1031170/1‑150 (MQ=255)

                                                                                                      gCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTAttt                                                >  2:145733/1‑150 (MQ=255)

                                                                                                       cTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCGTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTc                                               <  1:61283/150‑1 (MQ=255)

                                                                                                           atatTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTg                                           <  1:1514871/150‑1 (MQ=255)

                                                                                                              ttGTTTCTGCCCCCACTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTa                                        <  2:1227575/150‑1 (MQ=255)

                                                                                                               tGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCATCCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTaa                                       >  1:362428/1‑150 (MQ=255)

                                                                                                               tGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTaa                                       <  2:1243448/150‑1 (MQ=255)

                                                                                                               tGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTaa                                       <  2:641346/150‑1 (MQ=255)

                                                                                                               tGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTaa                                       <  2:688333/150‑1 (MQ=255)

                                                                                                               tGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTaa                                       <  1:1596272/150‑1 (MQ=255)

                                                                                                                gTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAAt                                      >  2:1050690/1‑150 (MQ=255)

                                                                                                                    cTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCatg                                  >  2:131985/1‑150 (MQ=255)

                                                                                                                      gCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCACCCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCatgat                                <  2:882171/150‑1 (MQ=255)

                                                                                                                       cccccAGTGGGGCGGGTAAATCCTGCCTGATTCAGGCTTTTTTAAAACCCCCCCCGTTGTATGACACCCAGGTTTCTGTTTCACAGACCACACGCCAACCGCGTCCGGGTGAAGTCCCCGGTGACCAGTATGTCTTTGTTAAACatgttg                               >  2:318543/1‑150 (MQ=255)

                                                                                                                        ccaaaGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCAACGGTGAACATTATTTCTTTGTTAATCAtgatga                              >  2:15032/5‑150 (MQ=255)

                                                                                                                           cAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAAtt                           <  2:596170/150‑1 (MQ=255)

                                                                                                                           cAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAAtt                           <  1:503736/150‑1 (MQ=255)

                                                                                                                           cAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTGTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGGAGTCCAGGGTGATCATTATTGCTGTGTTAATCTGGATGaagt                           >  2:317098/1‑148 (MQ=255)

                                                                                                                            aGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAAGATTATTTCTTTGTTAATCATGATGAAttt                          >  1:478224/1‑150 (MQ=255)

                                                                                                                             gTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAATTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTa                         >  1:929878/1‑150 (MQ=255)

                                                                                                                               ggCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTaaa                       >  2:1606761/1‑150 (MQ=255)

                                                                                                                               ggCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTaaa                       <  2:834130/150‑1 (MQ=255)

                                                                                                                                  gcgGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGaa                    <  2:10010/150‑1 (MQ=255)

                                                                                                                                      gTAAATCGATCCTGTTTGATTCTTTTTTAAAAACCCAACCGTTTTATGAAACCCAGGTTTCTGTTTCACACACCACACGCTAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCCTGATGAATTTAAAGAAATGa                <  1:859256/150‑1 (MQ=255)

                                                                                                                                       tAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGAt               <  1:636958/150‑1 (MQ=255)

                                                                                                                                       tAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACAAACCACACGCCAACCGCGTCCTTGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGAt               <  1:279730/150‑1 (MQ=255)

                                                                                                                                        aaaTCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGAtt              <  2:927068/150‑1 (MQ=255)

                                                                                                                                        aaaTCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGAtt              <  2:1067793/150‑1 (MQ=255)

                                                                                                                                         aataaaGCCTGAATCAGGCTTTATTAAAAACCCAACCGTTGTATCACACCCAGGTTTATGTTTCACACACAACACCCCACCAGCGACCAGGTGAAGTCAACGGTGAACATTATTTCTTTGTTAATAATGATGAATTTAAAGAAATGATTa             >  1:147423/6‑150 (MQ=255)

                                                                                                                                          aaccAGCCGGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACCCACCACACGCCAACCGCGTCCTGGAGAAGTCCACGGTGCACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAg            >  2:1119493/3‑150 (MQ=255)

                                                                                                                                            ccAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCa          <  2:1538814/150‑1 (MQ=255)

                                                                                                                                            ccAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGACCACGGTGAACATTATATAGTTGTTAATCATGATGAATTTAAAGAAATGATTAGCa          >  2:1066087/1‑150 (MQ=255)

                                                                                                                                             cAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCag         >  2:72508/1‑150 (MQ=255)

                                                                                                                                              aGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCaga        >  1:1423670/1‑150 (MQ=255)

                                                                                                                                               gCCTGATTCAGGCTTTATTAAAAACCCAACTGTTGTATGACACCCAGGTTTCTGTTTCACACTCCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCagag       <  2:758477/150‑1 (MQ=255)

                                                                                                                                               gCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGCCACCCAGGTTTCTGTTTCCCACACCACCCGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCagag       >  2:334336/1‑150 (MQ=255)

                                                                                                                                               gCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGCCCATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCagag       >  2:305502/1‑150 (MQ=255)

                                                                                                                                               gCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAGTGATTAGCagag       >  1:1183043/1‑150 (MQ=255)

                                                                                                                                               gCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCagag       <  1:160349/150‑1 (MQ=255)

                                                                                                                                               gCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCagag       >  1:650677/1‑150 (MQ=255)

                                                                                                                                                    aTTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCTTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCg  <  1:671831/150‑1 (MQ=255)

                                                                                                                                                    aTTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCg  >  2:968010/1‑150 (MQ=255)

                                                                                                                                                    aTTCACGCTTTATTAAAAACACAACCTTAGTATGACACACAGGTTTCTCATTCACACACCACACGCCAACCGCGACCTGGTGAAGTCCACGCTGAACATTATTTCTGTGTTAATCATGATGCATTTAAAGAAATGATTAGCAGAGATGCg  >  2:1125196/1‑150 (MQ=255)

GTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCG  >  NC_000913/3821340‑3821637

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