BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,906,231	C→A	6.3%	L182I (CTC→ATC)	mntP →	Mn(2(+)) exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,906,231	0	C	A	6.3%	529.7 / 12.9	157	L182I (CTC→ATC)	mntP	Mn(2(+)) exporter
Reads supporting (aligned to +/- strand): ref base C (84/63); new base A (5/5); total (89/68)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 160 total aligned reads displayed.

ACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGCGACG  >  NC_000913/1906082‑1906380

    tATCTCGACCGTATTGGCTATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATTGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGTAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTg                                                                                                                                                   <  2:1070087/150‑1 (MQ=255)

       cgcgACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAGGATCTGGAc                                                                                                                                                <  1:256114/150‑1 (MQ=255)

        gcgACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACg                                                                                                                                               >  2:155260/1‑150 (MQ=255)

        gcgACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAGGATCTGGACg                                                                                                                                               <  1:1379494/150‑1 (MQ=255)

         cgACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAGGATCTGGACGc                                                                                                                                              <  1:975209/150‑1 (MQ=255)

          gACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAGGATCTGGACGcc                                                                                                                                             <  2:2050894/150‑2 (MQ=255)

           aCCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCAc                                                                                                                                            >  2:2172924/1‑150 (MQ=255)

                      cATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTta                                                                                                                                 <  1:1017505/150‑1 (MQ=255)

                      cATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTta                                                                                                                                 >  2:1225862/1‑150 (MQ=255)

                      cATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTta                                                                                                                                 >  1:1662569/1‑150 (MQ=255)

                         tGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaata                                                                                                                              <  1:1184116/150‑1 (MQ=255)

                           gTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAg                                                                                                                            <  2:677263/150‑1 (MQ=255)

                           gTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAg                                                                                                                            >  2:968125/1‑150 (MQ=255)

                           gTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAg                                                                                                                            <  2:1649995/150‑1 (MQ=255)

                           gTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAg                                                                                                                            >  1:821211/1‑150 (MQ=255)

                              gtgCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAAc                                                                                                                         <  1:1272147/150‑1 (MQ=255)

                               tgCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACg                                                                                                                        >  2:1468124/1‑150 (MQ=255)

                                gCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTTGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGc                                                                                                                       >  2:1433020/1‑150 (MQ=255)

                                gCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGc                                                                                                                       >  2:1908169/1‑150 (MQ=255)

                                gCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTATCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGc                                                                                                                       >  2:1342650/1‑150 (MQ=255)

                                  aaCCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTg                                                                                                                     <  1:1259521/150‑1 (MQ=255)

                                   aCCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGc                                                                                                                    >  1:1046190/1‑150 (MQ=255)

                                     cTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCTATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCa                                                                                                                  >  2:1822601/1‑150 (MQ=255)

                                     cTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCa                                                                                                                  >  1:623970/1‑150 (MQ=255)

                                     cTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATACTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCa                                                                                                                  >  1:2095421/1‑150 (MQ=255)

                                      ttGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCag                                                                                                                 <  1:699712/150‑1 (MQ=255)

                                               tCAACATTAGGGATGATGGTTGGACGCTTTATCGGCTCAAATATTGGGAAAAAACCGGAAATTCTCGGAGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACCGTAAATAAGAACGCTGCCAGAGGTGAata                                                                                                        >  2:465356/1‑150 (MQ=255)

                                                 aaCATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATAtt                                                                                                      >  1:1530320/1‑150 (MQ=255)

                                                  aCATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTa                                                                                                     >  2:195141/1‑150 (MQ=255)

                                                     ttAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAg                                                                                                  <  1:2191903/150‑1 (MQ=255)

                                                     ttAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAg                                                                                                  <  2:1016367/150‑1 (MQ=255)

                                                      tAGGGATGATGGTTGGTCGCTTTATCAGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGt                                                                                                 >  2:2351581/1‑150 (MQ=255)

                                                          gatgatGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGt                                                                                             >  2:2404120/1‑150 (MQ=255)

                                                            tgatgGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCt                                                                                           >  1:1453112/1‑150 (MQ=255)

                                                             gatgGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTg                                                                                          <  2:1417786/150‑1 (MQ=255)

                                                               tggttggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGc                                                                                        <  1:2032822/150‑1 (MQ=255)

                                                               tggttggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGc                                                                                        <  1:1928174/150‑1 (MQ=255)

                                                               tggttggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGc                                                                                        <  1:407508/150‑1 (MQ=255)

                                                               tgggtggtCGCTTTTTCGGCTCAATTATTGGGAAAAACCCGTAAATTCTCGGCGGGGGGGTGGTTATGGGCTTCGGGGTCCATATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGc                                                                                        <  1:602832/150‑1 (MQ=255)

                                                                 gttggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAg                                                                                      >  2:2420440/1‑150 (MQ=255)

                                                                  ttggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGt                                                                                     <  1:155260/150‑1 (MQ=255)

                                                                    ggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCg                                                                                   <  2:1662569/150‑1 (MQ=255)

                                                                     gtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGa                                                                                  >  2:206673/1‑150 (MQ=255)

                                                                      tCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGCTCGGCATCGGCGTCCAGATCCTCTGGCCGCACTTCCACGGTTAATAAGAACGCTGCCAGCGGTGAATATTAAAGTCCGTCTGGCAGTCGac                                                                                 >  2:87516/1‑149 (MQ=255)

                                                                        ggtttATCGGCTCAATTATTGGGAAAAAAGCGGATATTCTCGGCGGGCTGGTGCTGATGGGCATCGGCGTCCAGATCCTCTGGACGCACTTGGTCGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAaca                                                                               <  1:465356/148‑1 (MQ=255)

                                                                        gCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTTATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAaca                                                                               <  1:474075/150‑1 (MQ=255)

                                                                         cTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAac                                                                                >  1:405729/1‑148 (MQ=255)

                                                                         cTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAac                                                                                <  2:405729/148‑1 (MQ=255)

                                                                           ttATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACActt                                                                            >  1:433007/1‑150 (MQ=255)

                                                                               cGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCttt                                                                        <  1:546736/150‑1 (MQ=255)

                                                                                ggCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCtttt                                                                       >  1:32906/1‑150 (MQ=255)

                                                                                ggCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCtttt                                                                       >  1:712264/1‑150 (MQ=255)

                                                                                 ggtCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTtg                                                                      >  2:833418/3‑150 (MQ=255)

                                                                                    cAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTtgctg                                                                   >  2:625808/1‑150 (MQ=255)

                                                                                        tattGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCCGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAg                                                               >  1:1023190/1‑150 (MQ=255)

                                                                                         attGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAgt                                                              >  1:136625/1‑150 (MQ=255)

                                                                                         attGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAgt                                                              >  1:1858287/1‑150 (MQ=255)

                                                                                           tGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTACGTCTGGCAGTCGAACACTTCTTTTGCTGCCAgtgt                                                            >  2:1618857/1‑150 (MQ=255)

                                                                                             ggAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGttt                                                          >  1:2472755/1‑150 (MQ=255)

                                                                                              gAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTg                                                         <  1:784420/150‑1 (MQ=255)

                                                                                                aaaaaGCGGAAATTCTCGGAGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGcc                                                       >  2:2316871/1‑150 (MQ=255)

                                                                                                  aaaGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAg                                                     <  1:195141/150‑1 (MQ=255)

                                                                                                  aaaGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAg                                                     >  1:1800837/1‑150 (MQ=255)

                                                                                                    aGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGAc                                                   <  1:335298/150‑1 (MQ=255)

                                                                                                      cGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACtt                                                 >  1:2408414/1‑150 (MQ=255)

                                                                                                       ggAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCCGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTc                                                >  2:2386019/1‑150 (MQ=255)

                                                                                                            ttCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCt                                           >  1:1908651/1‑150 (MQ=255)

                                                                                                               tcGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCTGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCg                                        <  1:1604853/150‑1 (MQ=255)

                                                                                                               tcGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCg                                        <  2:831519/150‑1 (MQ=255)

                                                                                                               tcGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCg                                        >  2:1008216/1‑150 (MQ=255)

                                                                                                               tcGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGAAGTCTGGcaaac                                        <  2:1329803/150‑5 (MQ=255)

                                                                                                               tcGGCGGGCTGGTGCGGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGAAGTCTGGcaaac                                        <  2:2097819/150‑5 (MQ=255)

                                                                                                                         ggTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGc                              <  1:2104581/150‑1 (MQ=255)

                                                                                                                         ggTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGc                              >  1:286778/1‑150 (MQ=255)

                                                                                                                           tGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCaa                            >  2:1546759/1‑150 (MQ=255)

                                                                                                                             cTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAAc                          <  2:712264/150‑1 (MQ=255)

                                                                                                                               gATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACgg                        >  2:561481/1‑150 (MQ=255)

                                                                                                                               gATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACgg                        >  1:2258108/1‑150 (MQ=255)

                                                                                                                               gATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACgg                        >  2:31033/1‑150 (MQ=255)

                                                                                                                                 tcggcatcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCg                      <  2:789229/150‑1 (MQ=255)

                                                                                                                                  cggcatcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGt                     >  1:1957649/1‑150 (MQ=255)

                                                                                                                                   ggcatcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAACCGGCggc                    >  1:2085897/1‑148 (MQ=255)

                                                                                                                                   ggcatcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTc                    >  2:301664/1‑150 (MQ=255)

                                                                                                                                         tgcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc              <  2:1065459/149‑1 (MQ=255)

                                                                                                                                         tgcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc              >  1:1065459/2‑150 (MQ=255)

                                                                                                                                         ggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc              <  1:537446/150‑1 (MQ=255)

                                                                                                                                         ggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc              >  2:1217038/1‑150 (MQ=255)

                                                                                                                                          gcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCa             >  2:1923124/1‑150 (MQ=255)

                                                                                                                                          gcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCa             >  1:1766113/1‑150 (MQ=255)

                                                                                                                                           cGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCaa            <  2:1046190/150‑1 (MQ=255)

                                                                                                                                            gTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCaat           <  1:2000776/150‑1 (MQ=255)

                                                                                                                                            gTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCaat           <  2:1268301/150‑1 (MQ=255)

                                                                                                                                            gTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCaat           <  2:875178/150‑1 (MQ=255)

                                                                                                                                             tCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCaata          <  1:1908169/150‑1 (MQ=255)

                                                                                                                                                aGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCaataaag       >  2:1828131/1‑148 (MQ=255)

                                                                                                                                                aGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATg       >  1:2467175/1‑150 (MQ=255)

                                                                                                                                                 gATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGc      <  1:910909/150‑1 (MQ=255)

                                                                                                                                                 gATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGACAGGCAAACGGCGTCATCTGCATTAATGc      <  2:697186/150‑1 (MQ=255)

                                                                                                                                                  aTCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGCg     <  2:2124457/150‑1 (MQ=255)

                                                                                                                                                  aTCCTCTGGACGCACTTCCACGGTTAAAAAGAACGCTGCCAGAGGTGAATATTAAATTCCGTATGGCAGTCGAACAATTCTTTTCCTGCCAGTGTATGACAGACTTATGGCTTCGCACTCCAGGCAAACGCCGTCATATGGAATAATGcc     >  1:2484987/1‑149 (MQ=255)

ACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGCGACG  >  NC_000913/1906082‑1906380

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