BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	945,097	T→C	S2P (TCC→CCC)	ycaD →	putative transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	945,097	0	T	C	99.5%	591.5 / NA	189	S2P (TCC→CCC)	ycaD	putative transport
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (86/102); minor base A (1/0); total (87/102)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

GGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTT  >  NC_000913/945007‑945186

gggATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGccc                                                                                         >  1:1507096/1‑93 (MQ=255)

gggATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACCGTTTTAATTCATGccc                                                                                         >  2:7538199/1‑93 (MQ=255)

 ggATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCa                                                                                        <  1:1888984/93‑1 (MQ=255)

 ggATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCa                                                                                        >  1:6982871/1‑93 (MQ=255)

  gATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCAc                                                                                       >  2:6143703/1‑93 (MQ=255)

  gATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCAc                                                                                       >  1:2654400/1‑93 (MQ=255)

  gATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCAc                                                                                       >  2:5101971/1‑93 (MQ=255)

    tatCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGt                                                                                     <  1:5286708/93‑1 (MQ=255)

     atCGCGTATTTTCTCCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTGCCTCACGGTTTTAATTCATGCCCACGta                                                                                    <  1:3223562/93‑1 (MQ=255)

     atCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGATTTAATTCATGCCCACGta                                                                                    >  1:4558650/1‑93 (MQ=255)

       cgcgTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGtata                                                                                  <  1:4360415/93‑1 (MQ=255)

           tATTTTCACCCATTGTCAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCa                                                                              <  1:4972599/93‑1 (MQ=255)

           tATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATTCCCACGTATACCCa                                                                              <  1:1998766/93‑1 (MQ=255)

           tATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCa                                                                              <  2:331839/93‑1 (MQ=255)

           tATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCa                                                                              <  1:2433057/93‑1 (MQ=255)

           tATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCa                                                                              >  2:5567068/1‑93 (MQ=255)

            aTTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAg                                                                             >  2:6546542/1‑93 (MQ=255)

             ttttCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGc                                                                            <  1:965057/93‑1 (MQ=255)

             ttttCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGc                                                                            >  1:3619385/1‑93 (MQ=255)

              tttCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGcc                                                                           <  2:2732305/93‑1 (MQ=255)

                tCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTg                                                                         <  1:3721319/93‑1 (MQ=255)

                  aCCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTc                                                                       <  2:3019806/93‑1 (MQ=255)

                   cccATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCa                                                                      <  2:1092683/93‑1 (MQ=255)

                   cccATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCa                                                                      >  1:62124/1‑93 (MQ=255)

                     cATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATg                                                                    <  1:6470143/93‑1 (MQ=255)

                     cATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATg                                                                    >  1:4515549/1‑93 (MQ=255)

                      aTTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGc                                                                   <  2:2554093/93‑1 (MQ=255)

                       ttGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTTACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCt                                                                  >  2:3160577/1‑93 (MQ=255)

                       ttGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCt                                                                  >  1:3543084/1‑93 (MQ=255)

                       ttGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCt                                                                  >  1:3439859/1‑93 (MQ=255)

                        tGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCtt                                                                 <  2:219045/93‑1 (MQ=255)

                        tGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCtt                                                                 <  1:4672112/93‑1 (MQ=255)

                        tGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCtt                                                                 >  2:289207/1‑93 (MQ=255)

                         taCAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATTCCCACGTATACCCAGCCTGTCATGCttt                                                                <  1:5859949/92‑1 (MQ=255)

                         gACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATTCCCACGTATACCCAGCCTGTCATGCttt                                                                <  2:1778304/93‑1 (MQ=255)

                         gACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCttt                                                                <  1:1143741/93‑1 (MQ=255)

                          aCAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCtttt                                                               >  1:4785745/1‑93 (MQ=255)

                           cAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCAAGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTg                                                              >  2:2242770/1‑93 (MQ=255)

                            aaTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGc                                                             >  1:5022881/1‑93 (MQ=255)

                            aaTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGc                                                             >  1:2883689/1‑93 (MQ=255)

                             aTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCt                                                            <  2:863957/93‑1 (MQ=255)

                                gttttGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTc                                                         <  2:2593864/92‑1 (MQ=255)

                                 ttttGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCt                                                        <  1:4140725/93‑1 (MQ=255)

                                 ttttGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCt                                                        <  2:6950736/93‑1 (MQ=255)

                                 ttttGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCt                                                        <  2:7467333/93‑1 (MQ=255)

                                  tttGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTg                                                       <  1:1139177/93‑1 (MQ=255)

                                   ttGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTgg                                                      <  2:4974143/93‑1 (MQ=255)

                                     ggCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGcc                                                    <  1:3448786/93‑1 (MQ=255)

                                     ggCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGcc                                                    <  2:7882384/93‑1 (MQ=255)

                                        ggTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGc                                                 <  2:7767814/93‑1 (MQ=255)

                                         gTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCt                                                >  1:7603953/1‑93 (MQ=255)

                                         gTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCt                                                >  1:7834207/1‑93 (MQ=255)

                                         gTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCt                                                <  1:6344322/93‑1 (MQ=255)

                                           ggCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCttt                                              >  1:1577561/1‑93 (MQ=255)

                                           ggCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCttt                                              >  2:3656752/1‑93 (MQ=255)

                                            gCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCtttt                                             <  2:1939300/93‑1 (MQ=255)

                                             cATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTttg                                            <  2:2834603/93‑1 (MQ=255)

                                             cATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTttg                                            <  2:5400505/93‑1 (MQ=255)

                                              atgatgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTttgt                                           >  2:5667075/1‑93 (MQ=255)

                                              atgatgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTttgt                                           >  1:2777929/1‑93 (MQ=255)

                                               tgatgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTttgtt                                          <  1:7228863/93‑1 (MQ=255)

                                                gatgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTttgttg                                         >  1:3861842/1‑93 (MQ=255)

                                                  tgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGGTGTCTGGCCTGCTTTTGTTGat                                       >  2:7397107/1‑92 (MQ=255)

                                                  tgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGAc                                       >  1:586818/1‑93 (MQ=255)

                                                    cgcATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGActc                                     <  2:5216590/93‑1 (MQ=255)

                                                    cgcATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGActc                                     <  2:5216591/93‑1 (MQ=255)

                                                     gcATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGActct                                    >  2:5457978/1‑93 (MQ=255)

                                                      cATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTg                                   <  2:1554376/93‑1 (MQ=255)

                                                      cATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTg                                   <  1:5467256/93‑1 (MQ=255)

                                                        tGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGc                                 >  1:4474167/1‑93 (MQ=255)

                                                         gAAATTTGAACTTCCTCACGATTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCg                                >  2:7615947/1‑93 (MQ=255)

                                                          aaaTTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGa                               <  2:663296/93‑1 (MQ=255)

                                                          aaaTTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGa                               <  1:3992483/93‑1 (MQ=255)

                                                          aaaTTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGa                               <  2:5139222/93‑1 (MQ=255)

                                                            aTTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGAtt                             <  2:4793885/93‑1 (MQ=255)

                                                             tttGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTg                            >  2:6384966/1‑93 (MQ=255)

                                                              ttGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGc                           <  1:7163696/93‑1 (MQ=255)

                                                                 aaCTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGgt                        <  2:4864719/93‑1 (MQ=255)

                                                                  aCTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGgtg                       >  1:1944164/1‑93 (MQ=255)

                                                                  aCTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGgtg                       <  1:7449480/93‑1 (MQ=255)

                                                                  aCTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGgtg                       >  1:6896964/1‑93 (MQ=255)

                                                                   cTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGgtgt                      <  1:2243866/93‑1 (MQ=255)

                                                                    ttCCTCACGGTTTTAATTCATGCCCACGTATACCAAGCCTGTCATGCTTTTGCTGTATGGCCTGCTTTTGTTGAATCTGGCGATTGCGGTGtt                     >  1:3748879/1‑93 (MQ=255)

                                                                      ccTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTaa                   <  1:5975365/93‑1 (MQ=255)

                                                                       cTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTaaa                  >  1:3186143/1‑93 (MQ=255)

                                                                         cACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATa                >  2:275892/1‑93 (MQ=255)

                                                                         cACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATa                <  2:6845892/93‑1 (MQ=255)

                                                                          aCGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATac               <  1:5103768/93‑1 (MQ=255)

                                                                          aCGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATac               >  1:6641157/1‑93 (MQ=255)

                                                                          aCGGTGTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATac               <  1:7810582/93‑1 (MQ=255)

                                                                           cGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATaca              <  2:4075946/93‑1 (MQ=255)

                                                                           cGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATaca              <  1:7075197/93‑1 (MQ=255)

                                                                           cGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATaca              <  1:6258832/93‑1 (MQ=255)

                                                                            ggTTTTAATTTCTGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATacac             <  2:6388728/93‑1 (MQ=255)

                                                                               tttAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCg          <  1:1253887/93‑1 (MQ=255)

                                                                               tttAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCg          <  1:6359970/93‑1 (MQ=255)

                                                                                ttAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGt         >  1:6728996/1‑93 (MQ=255)

                                                                                  aaTTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGc       >  2:3480985/1‑93 (MQ=255)

                                                                                   aTTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGcc      >  2:4926841/1‑93 (MQ=255)

                                                                                       aTGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCtt  <  1:3368930/93‑1 (MQ=255)

GGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTT  >  NC_000913/945007‑945186

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`