BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	824,949	T→C	5.6%	L107P (CTG→CCG)	ybhQ →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	824,949	0	T	C	5.6%	362.0 / 9.1	125	L107P (CTG→CCG)	ybhQ	inner membrane protein
Reads supporting (aligned to +/- strand): ref base T (59/57); new base C (5/2); total (64/61)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

Only 100 of 126 total aligned reads displayed.

TTTCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGA  >  NC_000913/824801‑825098

tttCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTg                                                                                                                                                      >  1:755189/1‑150 (MQ=255)

 ttCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTgg                                                                                                                                                     >  1:1004849/1‑150 (MQ=255)

 ttCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTgg                                                                                                                                                     >  1:385822/1‑150 (MQ=255)

  tCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGc                                                                                                                                                    >  1:139620/1‑150 (MQ=255)

   cAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCGGGACGAGCTGACCACGACCTGGGATTTTGGCGCAGCGCTGATGGTGCTGTGGCTGTTGTCCCGCGTGCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCa                                                                                                                                                   <  1:1044117/150‑1 (MQ=255)

    aGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCaa                                                                                                                                                  >  2:518813/1‑150 (MQ=255)

          aCCCCGAGCAACGCTGGCGTGAGGTGGCAGCCGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTg                                                                                                                                            >  2:1128498/1‑150 (MQ=255)

          aCCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTg                                                                                                                                            <  1:1677246/150‑1 (MQ=255)

          aCCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTg                                                                                                                                            <  2:435526/150‑1 (MQ=255)

          aCCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTg                                                                                                                                            <  2:239786/150‑1 (MQ=255)

            cccGAGCAACGCTGGCGTGAGTGGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGcc                                                                                                                                          <  2:145417/150‑1 (MQ=255)

             ccGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGccg                                                                                                                                         >  2:95533/1‑150 (MQ=255)

              cGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGccgg                                                                                                                                        >  2:155959/1‑150 (MQ=255)

                  cAACGCTGGCGTGAGGCGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGccggcccg                                                                                                                                    >  2:1449025/1‑150 (MQ=255)

                   aaCGCTGGCGTGAGGGGGCAGACGTGCGGGTAGAGCTGACCACGACCTGGTATTATGGCGCAGCGCGGATTGTGCTGTGGCTGTTGTCCCGCGGGCTGGTAAACAACTTTTTGCTGGCAATTGCAGGGCCGGCAAGGATTGcctgccctg                                                                                                                                   <  1:230359/150‑3 (MQ=255)

                    ccgcgGGCGGGAGGTGGCAGACGTGCTGGAAGCGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGAGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGccggcccggc                                                                                                                                  <  1:923904/149‑1 (MQ=255)

                    aCGCTGGCGTGAGGTGGCAGACGCGCTGGACGAGCGGCCCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGGTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGccggcccggc                                                                                                                                  <  1:1143805/150‑1 (MQ=255)

                      gcgggCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCCCGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTGCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCcggcgg                                                                                                                                <  1:884891/147‑1 (MQ=255)

                      gCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCcggcgg                                                                                                                                <  1:1769490/150‑1 (MQ=255)

                      gCTGGCGGGAGGTGGCAGACGAGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTAGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCcggcgg                                                                                                                                <  1:633270/150‑1 (MQ=255)

                       cTGGCGTGAGGGGGCAGACGTGCTGGAAGAGCGGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCGGTTGTCCCTCGTTCTGGAAAACAACTTTATGCTGGCAATTGCAGGGCCGGCAAGGATTGccagccctgcaat                                                                                                                               <  1:1149382/150‑4 (MQ=255)

                         ggCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAg                                                                                                                             >  2:638940/1‑150 (MQ=255)

                           cGTGAGGTGGCAGACGTGCTGGAAGAGCTGCCCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAgtg                                                                                                                           >  2:430198/1‑150 (MQ=255)

                           cGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAgtg                                                                                                                           >  1:1408468/1‑150 (MQ=255)

                                 ggggCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGAGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTg                                                                                                                     <  2:204948/148‑1 (MQ=255)

                                 gTGGCAGACGTGCTGGAAGAGCTGACCACGACCTTGTATTTTGGAGCAGCGCTTATTGTGCTGTTGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCTGCCCGGCGGTAGTGTCtttt                                                                                                                     <  2:975714/150‑2 (MQ=255)

                                  tGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGc                                                                                                                    <  2:1746237/150‑1 (MQ=255)

                                   ggCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCt                                                                                                                   >  1:387559/1‑150 (MQ=255)

                                     cAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTgg                                                                                                                 >  1:1020342/1‑150 (MQ=255)

                                     cAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTgg                                                                                                                 >  1:1525029/1‑150 (MQ=255)

                                      aGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGc                                                                                                                <  2:547441/150‑1 (MQ=255)

                                          gTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGaaa                                                                                                            <  2:917642/150‑1 (MQ=255)

                                            gatGGAAGAGCTGACAACTACCTGGTATTTTTGAGCAGCGATGATTGTGCTGTGGCTGTTGTCACGCATTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGa                                                                                                          >  2:53201/3‑150 (MQ=255)

                                             cTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGAt                                                                                                         >  1:452464/1‑150 (MQ=255)

                                             cTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGAt                                                                                                         <  2:1297309/150‑1 (MQ=255)

                                                gAAGATCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCAGTGGCTGTTGTCCCGCGTTCCGGAAAACAACTTTTTGCTGGCAATTGCAGGGCGGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATaag                                                                                                      <  2:370620/150‑1 (MQ=255)

                                                 aaGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATaaga                                                                                                     <  2:435940/150‑1 (MQ=255)

                                                     gCTGACCACGACGTGGTATTTTGGCGCAGCGCTGTTTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGtt                                                                                                 <  1:1454301/150‑1 (MQ=255)

                                                     gCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGtt                                                                                                 >  2:564112/1‑150 (MQ=255)

                                                       tGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTAc                                                                                               >  2:490468/1‑150 (MQ=255)

                                                        gACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACa                                                                                              >  2:422808/1‑150 (MQ=255)

                                                         aCCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCTTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCTAAAGATTAGAAGTTACAt                                                                                             >  1:158267/1‑150 (MQ=255)

                                                          ccACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATc                                                                                            >  2:486835/1‑150 (MQ=255)

                                                          ccACGACCTGGTATTTTGGCGCAGCGCTGATAGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATc                                                                                            >  2:1448106/1‑150 (MQ=255)

                                                           cACGACCTGGTATTTTGGCTCAGCGCTGATTGTGATGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCa                                                                                           <  2:1390490/150‑1 (MQ=255)

                                                           cACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCa                                                                                           >  1:1212284/1‑150 (MQ=255)

                                                              gACCTGGTATTTTGGCGGAGCGCTGATTGTGCTGGGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTAGCTGGCAATTGCAGGGCTGGCAATCCTTGCGGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCAGCt                                                                                        <  2:1491943/150‑1 (MQ=255)

                                                                 cTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTAc                                                                                     >  1:1089010/1‑150 (MQ=255)

                                                                 cTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTAc                                                                                     >  2:1453471/1‑150 (MQ=255)

                                                                 cTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTAc                                                                                     >  2:1517993/1‑150 (MQ=255)

                                                                  tGGTATTTTGGCGCAGCGCTGATTGTCCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACg                                                                                    <  2:116244/150‑1 (MQ=255)

                                                                   ggTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGt                                                                                   >  1:1314727/1‑150 (MQ=255)

                                                                    gTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTc                                                                                  <  2:767927/150‑1 (MQ=255)

                                                                        tttGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTaaa                                                                              <  2:1790473/150‑1 (MQ=255)

                                                                             cgcAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATcg                                                                         >  1:1335901/1‑150 (MQ=255)

                                                                             cgcAGCGCTGATTGTGCGGTGGCGGGTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATcg                                                                         <  2:1577345/150‑1 (MQ=255)

                                                                              gcAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATcgc                                                                        <  2:268653/150‑1 (MQ=255)

                                                                               cAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATcgcg                                                                       <  1:1322570/150‑1 (MQ=255)

                                                                                  cgcTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTAc                                                                    >  2:640013/1‑150 (MQ=255)

                                                                                   gcTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTATTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTAcg                                                                   >  2:852733/1‑150 (MQ=255)

                                                                                       aTTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTAcgccgc                                                               <  2:829762/150‑1 (MQ=255)

                                                                                       aTTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTAcgccga                                                               >  2:865872/1‑149 (MQ=255)

                                                                                            gCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGAcc                                                          >  2:354338/1‑150 (MQ=255)

                                                                                              tgtgGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCt                                                        >  1:886064/1‑150 (MQ=255)

                                                                                              tgtgGCTGTAGTCCCGCGTTCTGGATAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATATGAAGTTACATCACCTTTCGTCAAAACATCGAGTACGCAGCTGACCCt                                                        >  1:880229/1‑150 (MQ=255)

                                                                                               gtgGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTg                                                       >  2:1824890/1‑150 (MQ=255)

                                                                                                tgGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTTGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGt                                                      <  2:788812/150‑1 (MQ=255)

                                                                                                tgGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCTAATGCAGGGCTGGCAGTCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATTAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGt                                                      <  2:1500478/150‑1 (MQ=255)

                                                                                                tgGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGt                                                      <  1:1824890/150‑1 (MQ=255)

                                                                                                  gCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCg                                                    <  2:363523/150‑1 (MQ=255)

                                                                                                   cTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGt                                                   <  1:1780746/150‑1 (MQ=255)

                                                                                                    tgttgtCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTg                                                  <  1:1721731/150‑1 (MQ=255)

                                                                                                       tgtCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGcc                                               <  1:1038093/150‑1 (MQ=255)

                                                                                                         gcccGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGAAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATTGCGTACGCCGCTGACCCTGTCGTGGCCGt                                             <  2:1429707/149‑1 (MQ=255)

                                                                                                          cccGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGtt                                            >  2:122290/1‑150 (MQ=255)

                                                                                                             gcgTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATc                                         >  2:691751/1‑150 (MQ=255)

                                                                                                             gcgTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCATTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCTAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCTCGTACGCCGCTGCCCCTGTCGTGGCCGatata                                         >  2:1816384/1‑149 (MQ=255)

                                                                                                              cgTTCTGGTATACAAATTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGGCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCa                                        <  2:754641/150‑1 (MQ=255)

                                                                                                              cgTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTTGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCa                                        <  1:1184522/150‑1 (MQ=255)

                                                                                                                  cTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAg                                    <  1:376174/150‑1 (MQ=255)

                                                                                                                   tGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGt                                   >  2:904513/1‑150 (MQ=255)

                                                                                                                    ggAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTa                                  >  2:221878/1‑150 (MQ=255)

                                                                                                                          caacTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGcc                            >  2:1232782/1‑150 (MQ=255)

                                                                                                                              tttttGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACaa                        <  1:1522091/150‑1 (MQ=255)

                                                                                                                                tttGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAAc                      <  2:1089010/150‑1 (MQ=255)

                                                                                                                                    cTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTagt                  >  2:556385/1‑148 (MQ=255)

                                                                                                                                    cTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCCTTATCACCAGTAGCGGACACAAACTAtt                  >  1:199811/1‑150 (MQ=255)

                                                                                                                                     tGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTc                 >  2:1506910/1‑150 (MQ=255)

                                                                                                                                      ggCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTcc                <  2:568415/150‑1 (MQ=255)

                                                                                                                                      ggCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGAGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTcc                <  1:155959/150‑1 (MQ=255)

                                                                                                                                        cAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGTGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCcac              <  2:156684/150‑1 (MQ=255)

                                                                                                                                        cAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCcac              >  2:112522/1‑150 (MQ=255)

                                                                                                                                         aaTTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCcacc             >  2:1642622/1‑149 (MQ=255)

                                                                                                                                          aTTGCAGTGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCcacac            <  2:139620/150‑1 (MQ=255)

                                                                                                                                           ccgccgggCCGGCAAGGATTGCCAGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCcccacc           >  1:564689/6‑149 (MQ=255)

                                                                                                                                              cAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAata        >  1:1458878/1‑150 (MQ=255)

                                                                                                                                               cgggCCGGCAAGGATTGCCAGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAatat       >  2:763702/2‑150 (MQ=255)

                                                                                                                                               cgggCCGGCAAGGATTGCCAGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAatat       >  2:608982/2‑150 (MQ=255)

                                                                                                                                                  gCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCa    <  1:73158/150‑1 (MQ=255)

                                                                                                                                                  gCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCa    <  2:1298857/150‑1 (MQ=255)

TTTCAGCGCCACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGA  >  NC_000913/824801‑825098

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