BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	858,424	A→T	5.3%	L213Q (CTG→CAG) ‡	ybiU ←	DUF1479 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	858,424	0	A	T	5.3%	310.7 / 3.5	104	L213Q (CTG→CAG) ‡	ybiU	DUF1479 family protein
Reads supporting (aligned to +/- strand): ref base A (50/48); new base T (2/4); total (52/52)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

AGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGA  >  NC_000913/858275‑858573

aGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagca                                                                                                                                                       >  2:132748/1‑150 (MQ=255)

aGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagca                                                                                                                                                       <  1:206922/150‑1 (MQ=255)

 gCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcag                                                                                                                                                      <  2:523338/150‑1 (MQ=255)

  cccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcagc                                                                                                                                                     >  2:465984/1‑150 (MQ=255)

      ggAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGcagccagc                                                                                                                                                 >  1:581482/1‑150 (MQ=255)

          tGTCCGAAACACTTAACATATGGTCCTGTTGTCAACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGCCTCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGAAACGCTGGAAGCAGCCAGAGTTc                                                                                                                                             >  2:1614086/1‑150 (MQ=255)

           gTCCTAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTcc                                                                                                                                            <  2:1175634/150‑1 (MQ=255)

           gTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTcc                                                                                                                                            >  2:1636771/1‑150 (MQ=255)

           gTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTcc                                                                                                                                            <  2:1559153/150‑1 (MQ=255)

             ccGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAg                                                                                                                                          <  1:276501/150‑1 (MQ=255)

              cGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGt                                                                                                                                         >  2:759434/1‑150 (MQ=255)

              cGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGAGCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGtatcagc                                                                                                                                         <  2:1598999/150‑7 (MQ=255)

               gAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGtatcagcg                                                                                                                                        <  2:226866/150‑8 (MQ=255)

               gAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGtatcagcg                                                                                                                                        <  1:1152952/150‑8 (MQ=255)

                 aaCACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGcc                                                                                                                                      <  2:1044081/150‑1 (MQ=255)

                   cacGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGcccc                                                                                                                                    >  1:1485730/1‑150 (MQ=255)

                    ccGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGAGCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCg                                                                                                                                   <  1:1168514/149‑1 (MQ=255)

                       gAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAg                                                                                                                                >  1:1469955/1‑150 (MQ=255)

                          cATTTGGTGGTGTTGTCCACAGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCg                                                                                                                             >  1:357448/1‑150 (MQ=255)

                           aTTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCgg                                                                                                                            >  1:1183950/1‑150 (MQ=255)

                            tttGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGt                                                                                                                           <  1:59591/150‑1 (MQ=255)

                             ttGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTa                                                                                                                          <  1:1138310/150‑1 (MQ=255)

                             ttGGTGGTGTTGTCCACCGTGTACTCTTAAACTTCCGTACGATGTGCCGCATGACAGGGATCATTTTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGATGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCTGAGTCGGTa                                                                                                                          <  1:1261654/150‑1 (MQ=255)

                             ttGGTGGGGTTGTCCACCGTGTACTCTTCATCTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTa                                                                                                                          <  1:1695342/150‑1 (MQ=255)

                              tggtggTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTAt                                                                                                                         >  1:1067686/1‑150 (MQ=255)

                               ggtggtGTTGTCGACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATg                                                                                                                        <  1:1137841/150‑1 (MQ=255)

                                 tggtgTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATgcg                                                                                                                      >  2:1027619/1‑150 (MQ=255)

                                        gtCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAg                                                                                                               >  1:454049/1‑150 (MQ=255)

                                          ccACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGAc                                                                                                             >  2:222221/1‑150 (MQ=255)

                                           cACCGTGTACTGTTCAACTTCCGTACGATTTGCCGCATGCCAGGGATCATATTGAGCCAGATTGCCATTAAAGACTTTGGCGAAAACGCGCTGATTCGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGAcc                                                                                                            <  1:1737672/150‑1 (MQ=255)

                                           cACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGAcc                                                                                                            <  2:207784/150‑1 (MQ=255)

                                             ccGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCtt                                                                                                          >  2:1499304/1‑150 (MQ=255)

                                                 gtACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGa                                                                                                      <  1:486530/150‑1 (MQ=255)

                                                  tACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAg                                                                                                     <  1:1155692/150‑1 (MQ=255)

                                                     tctTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTc                                                                                                  <  2:1339213/150‑1 (MQ=255)

                                                      ctTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCg                                                                                                 >  2:1257859/1‑150 (MQ=255)

                                                       ttCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTACGGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGt                                                                                                <  1:924529/150‑1 (MQ=255)

                                                           aCTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCg                                                                                            >  1:1428672/1‑150 (MQ=255)

                                                           aCTTCCGTACGATGTGCCGCATGCCAGGGATCATAATGCGCCAGATTGCAATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTTCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCg                                                                                            <  2:1379255/150‑1 (MQ=255)

                                                             ttCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCggg                                                                                          >  1:1734094/1‑150 (MQ=255)

                                                             ttCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCggg                                                                                          <  1:1551812/150‑1 (MQ=255)

                                                                   aCGATGTGCCGCCTGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGAcg                                                                                    <  2:1701947/150‑1 (MQ=255)

                                                                   aCGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGAcg                                                                                    >  2:1212155/1‑150 (MQ=255)

                                                                     gATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGGCGGTATGCGCTCCAAGACCTTTGGAGGGCGTTCCGGGCGGAcggc                                                                                  <  2:1095829/150‑1 (MQ=255)

                                                                         tgCCGCATGCCAGGGATCATATTGCGCCAGATTGACATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGGCGTTCCGGGCGGAcggcatcg                                                                              <  2:907376/150‑5 (MQ=255)

                                                                         ggCCGCACGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGAcggcggcg                                                                              <  2:320507/149‑1 (MQ=255)

                                                                           ccGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGATGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGa                                                                            >  1:945627/1‑150 (MQ=255)

                                                                           ccGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGa                                                                            >  1:799307/1‑150 (MQ=255)

                                                                               aTGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACg                                                                        <  2:1433659/150‑1 (MQ=255)

                                                                                  ccAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTc                                                                     >  2:1553644/1‑150 (MQ=255)

                                                                                   cAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCa                                                                    >  2:392904/1‑150 (MQ=255)

                                                                                     gggATCATATTGCGCCAGATTGCCCTTAAAGACGTTGGCGAACACGCGCTGATACGCTGGAAGCAGCCAGCGGTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAgg                                                                  <  1:234037/150‑1 (MQ=255)

                                                                                      ggATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAggg                                                                 <  2:357448/150‑1 (MQ=255)

                                                                                      ggATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAggg                                                                 >  2:752572/1‑150 (MQ=255)

                                                                                       gATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGt                                                                <  1:384907/150‑1 (MQ=255)

                                                                                        aTCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCTGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTa                                                               <  2:933673/150‑1 (MQ=255)

                                                                                        aTCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTa                                                               >  2:1530549/1‑150 (MQ=255)

                                                                                           atatTGCGACAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGCAAGAt                                                            >  2:1312742/1‑150 (MQ=255)

                                                                                            tatTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATg                                                           >  1:472634/1‑150 (MQ=255)

                                                                                             atTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGa                                                          >  2:1251704/1‑150 (MQ=255)

                                                                                             atTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGa                                                          >  2:1480581/1‑150 (MQ=255)

                                                                                                 cgcCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCt                                                      >  2:422922/1‑150 (MQ=255)

                                                                                                  gcCAGATTGCCATTAAAGACGTTGGCGACAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTc                                                     <  2:241148/150‑1 (MQ=255)

                                                                                                  gcCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCTGACGGCGGCGGATACGGTCAGGGTAGATGACGCTc                                                     >  1:1626453/1‑150 (MQ=255)

                                                                                                   ccAGATTGCCATTAACGACGTGGGCGAAAACGCGCGGATACGCTGGAAGCAGCCAGCGTTCCAGGGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGGCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCGca                                                    <  1:1146456/150‑1 (MQ=255)

                                                                                                   ccAGATTGCCATTAAAGCCGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTca                                                    <  1:1396123/150‑1 (MQ=255)

                                                                                                    aaGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCGCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGGCAGGGTAGATGACGCTcac                                                   <  2:865694/149‑1 (MQ=255)

                                                                                                        ttGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGCTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATcc                                               >  1:464046/1‑150 (MQ=255)

                                                                                                           ccATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCggg                                            <  1:427610/150‑1 (MQ=255)

                                                                                                           ccATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCggg                                            >  2:609628/1‑150 (MQ=255)

                                                                                                            cATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGt                                           >  2:1035517/1‑150 (MQ=255)

                                                                                                             aTTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGtt                                          >  1:1222405/1‑150 (MQ=255)

                                                                                                              ttAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGGCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTa                                         <  1:759434/150‑1 (MQ=255)

                                                                                                                 aaGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCGCATCCGGGTTacac                                      >  1:908130/1‑150 (MQ=255)

                                                                                                                 aaGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAAc                                      >  2:1380382/1‑150 (MQ=255)

                                                                                                                      gTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCCTTCCAGTGCCCCGGAGTCGGTATGCACTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTg                                 >  2:795418/1‑150 (MQ=255)

                                                                                                                       ttGGCGAAAACGCGCTGATACGCGGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGc                                <  1:1678478/150‑1 (MQ=255)

                                                                                                                         ggCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCtt                              >  1:123456/1‑150 (MQ=255)

                                                                                                                          gCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCttt                             >  1:823018/1‑150 (MQ=255)

                                                                                                                           cGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTc                            >  1:1114659/1‑150 (MQ=255)

                                                                                                                           cGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTc                            <  2:1183950/150‑1 (MQ=255)

                                                                                                                           cGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGCACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTc                            >  2:1077907/1‑150 (MQ=255)

                                                                                                                             aaaaCGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCa                          >  1:1761730/1‑150 (MQ=255)

                                                                                                                              aaaCGCGCAGATACGCTGGAAGCAGCCAGCGGTCCAGGGCCCCGGCGGCGGGATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCAt                         <  1:318850/150‑1 (MQ=255)

                                                                                                                                aCGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCa                       <  1:1416127/150‑1 (MQ=255)

                                                                                                                                 cgcgCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCAc                      <  2:1439246/150‑1 (MQ=255)

                                                                                                                                 cgcgCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCAc                      >  2:1518475/1‑150 (MQ=255)

                                                                                                                                   cgcTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACtt                    <  2:594855/150‑1 (MQ=255)

                                                                                                                                      tggaaCGCTGGCTGCTTCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCa                 >  1:254104/5‑150 (MQ=255)

                                                                                                                                       ggaaCGCTGGCTGCTTCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCaa                >  1:1296014/4‑150 (MQ=255)

                                                                                                                                        aTACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGGTCCGGGCGGCCGGCGGCGGAAACGGTCAGGGTAGATGACGCTCACATCCGGGTGACACCTTTGCTTTCCATCACTTTCaaa               >  1:709491/1‑150 (MQ=255)

                                                                                                                                          aCGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATg             <  1:1212155/150‑1 (MQ=255)

                                                                                                                                            gCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTc           >  1:1709710/1‑150 (MQ=255)

                                                                                                                                             cTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTcc          <  2:1818472/150‑1 (MQ=255)

                                                                                                                                               ggAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCcac        >  1:1712902/1‑150 (MQ=255)

                                                                                                                                               ggAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCAATCACTTTCAAATGTcccc        >  2:221875/1‑148 (MQ=255)

                                                                                                                                                 aaGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAg      <  2:1769110/150‑1 (MQ=255)

                                                                                                                                                    cagccagcGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACg   >  1:1258615/1‑150 (MQ=255)

                                                                                                                                                    cagccagcGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACg   <  2:419959/150‑1 (MQ=255)

                                                                                                                                                     agccagcGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGa  >  2:778029/1‑150 (MQ=255)

AGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGA  >  NC_000913/858275‑858573

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