BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	824,957	C→A	5.2%	L110I (CTT→ATT)	ybhQ →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	824,957	0	C	A	5.2%	391.7 / 7.5	124	L110I (CTT→ATT)	ybhQ	inner membrane protein
Reads supporting (aligned to +/- strand): ref base C (56/61); new base A (5/2); total (61/63)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

Only 100 of 126 total aligned reads displayed.

ACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATC  >  NC_000913/824811‑825106

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)

          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)

             cTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGt                                                                                                                                       <  1:1358859/150‑1 (MQ=255)

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

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

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

                      gggggCAGACGTGCTGTAAGAGCTGACCACGACCTTGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCGGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCttt                                                                                                                              <  2:341303/147‑1 (MQ=255)

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

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

                       gTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTg                                                                                                                             >  1:499031/1‑150 (MQ=255)

                        tGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGc                                                                                                                            <  2:1746237/150‑1 (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:1644336/1‑150 (MQ=255)

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

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

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

                                       aaGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATaaga                                                                                                             >  2:427299/1‑150 (MQ=255)

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

                                             tGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTAc                                                                                                       >  2:490468/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                                                                                             >  2:1517993/1‑150 (MQ=255)

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

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

                                                        tGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACg                                                                                            <  1:1489480/150‑1 (MQ=255)

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

                                                         ggTATTTTGGCGGAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGt                                                                                           >  1:857881/1‑150 (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)

                                                                        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)

                                                                                gtgCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATAACCTTAAGTCTAAACATCGCGTACTCAGCTGa                                                                    >  1:1210768/1‑150 (MQ=255)

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

                                                                                    tgtgGCTGTTGTCCCGCGTTCTTGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGACCGGCGGTAGTGTCTTTGCTGGCGAAAGATAACAAGATACATCACCTTACGTCTTAACATAGCGTACGCCGCTGACCCt                                                                >  2:442349/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)

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

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

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

                                                                                              gtCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCg                                                      <  1:307995/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)

                                                                                                                     ttttGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACaaa                               <  1:1817309/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)

                                                                                                                              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:608982/2‑150 (MQ=255)

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

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

                                                                                                                                           ggCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGaa         >  1:1024686/1‑149 (MQ=255)

                                                                                                                                              aaTCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCg      <  1:1094124/150‑1 (MQ=255)

                                                                                                                                              aaTCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCg      <  1:1128498/150‑1 (MQ=255)

                                                                                                                                                tCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCa    <  1:1453471/150‑1 (MQ=255)

                                                                                                                                                 ccTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCAt   >  1:1007995/1‑150 (MQ=255)

                                                                                                                                                  cTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATc  >  1:1560099/1‑150 (MQ=255)

ACCCCGAGCAACGCTGGCGTGAGGTGGCAGACGTGCTGGAAGAGCTGACCACGACCTGGTATTTTGGCGCAGCGCTGATTGTGCTGTGGCTGTTGTCCCGCGTTCTGGAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATC  >  NC_000913/824811‑825106

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