BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,854,698	T→C	Y445C (TAC→TGC)	yidJ ←	putative sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,854,698	0	T	C	99.5%	1324.4 / ‑3.9	384	Y445C (TAC→TGC)	yidJ	putative sulfatase
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (178/204); total (180/204)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 395 total aligned reads displayed.

ACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCG  >  NC_000913/3854608‑3854787

aCTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCaa                                                                                         <  2:7143441/93‑1 (MQ=255)

 cTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAAc                                                                                        <  2:1905210/93‑1 (MQ=255)

 cTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAAc                                                                                        <  2:573412/93‑1 (MQ=255)

   ggCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTg                                                                                      >  1:509036/1‑93 (MQ=255)

     cGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCg                                                                                    <  2:302470/93‑1 (MQ=255)

      gAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCgg                                                                                   >  2:8205993/1‑93 (MQ=255)

       aaTAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGa                                                                                  >  1:5047305/1‑93 (MQ=255)

       aaTAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGa                                                                                  <  2:1731747/93‑1 (MQ=255)

        aTAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGaa                                                                                 <  1:7248259/93‑1 (MQ=255)

        aTAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGaa                                                                                 >  2:3633544/1‑93 (MQ=255)

        aTAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGaa                                                                                 >  2:4913082/1‑93 (MQ=255)

         tAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAAc                                                                                <  1:7309616/93‑1 (MQ=255)

          aGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACg                                                                               <  1:1835654/93‑1 (MQ=255)

          aGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACg                                                                               >  2:525783/1‑93 (MQ=255)

            ccATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGa                                                                             <  1:113800/93‑1 (MQ=255)

            ccATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGa                                                                             <  2:768775/93‑1 (MQ=255)

            ccATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGa                                                                             <  1:131255/93‑1 (MQ=255)

              aTCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATc                                                                           >  2:6064861/1‑93 (MQ=255)

              aTCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATc                                                                           >  2:6936433/1‑93 (MQ=255)

                 ttGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATcgcg                                                                        >  1:2207642/1‑93 (MQ=255)

                 ttGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATcgcg                                                                        <  1:5609669/93‑1 (MQ=255)

                  tgtgGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGa                                                                       >  2:6567676/1‑93 (MQ=255)

                   gtggacgtggacGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGaa                                                                      <  1:3872222/93‑1 (MQ=255)

                   gtggacgtggacGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGaa                                                                      <  1:4933365/93‑1 (MQ=255)

                   gtggacgtggacGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGaa                                                                      <  1:1419583/93‑1 (MQ=255)

                    gggacgtggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAAt                                                                     <  1:3063097/92‑1 (MQ=255)

                     ggacgtggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAAtt                                                                    >  2:2199005/1‑93 (MQ=255)

                          tggacgAAACGCCCCCATCCAGCGCGGCCGTGCTTCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTc                                                               <  1:8523761/93‑1 (MQ=255)

                          tggacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTc                                                               >  2:121791/1‑93 (MQ=255)

                            gacgAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCa                                                             <  1:2735478/93‑1 (MQ=255)

                              cgAAACTCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACGCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATg                                                           <  1:1575307/93‑1 (MQ=255)

                               gAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGt                                                          <  2:3856761/93‑1 (MQ=255)

                                aaaCGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTa                                                         >  1:756243/1‑93 (MQ=255)

                                aaaCGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGAGCATTTTGTCCATGTa                                                         >  2:3898335/1‑93 (MQ=255)

                                 aaCGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTaa                                                        <  2:1380209/93‑1 (MQ=255)

                                  aCGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAAt                                                       <  1:6633600/93‑1 (MQ=255)

                                   cGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATc                                                      >  2:2947128/1‑93 (MQ=255)

                                     cccccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCa                                                    >  1:5403257/1‑93 (MQ=255)

                                      ccccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCaa                                                   <  1:2478218/93‑1 (MQ=255)

                                      ccccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCaa                                                   >  1:4364840/1‑93 (MQ=255)

                                      ccccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCaa                                                   >  1:1967430/1‑93 (MQ=255)

                                       cccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCCAt                                                  >  1:6001448/1‑93 (MQ=255)

                                       cccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAAt                                                  >  2:7367926/1‑93 (MQ=255)

                                       cccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAAt                                                  <  1:8663787/93‑1 (MQ=255)

                                        ccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCCATa                                                 >  2:8347947/1‑93 (MQ=255)

                                        ccATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATa                                                 >  1:6597038/1‑93 (MQ=255)

                                         cATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATaa                                                <  2:6321540/93‑1 (MQ=255)

                                         cATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATaa                                                <  1:5655989/93‑1 (MQ=255)

                                          aTCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAg                                               >  1:6126613/1‑93 (MQ=255)

                                           tCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAgg                                              <  1:4061434/93‑1 (MQ=255)

                                           tCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAgg                                              >  1:6223155/1‑93 (MQ=255)

                                           tCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAgg                                              >  1:7148120/1‑93 (MQ=255)

                                             cAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCg                                            <  2:4883789/93‑1 (MQ=255)

                                             cAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCg                                            <  1:4656530/93‑1 (MQ=255)

                                              aGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGt                                           >  1:3611975/1‑93 (MQ=255)

                                               gcgcGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTc                                          >  1:1029310/1‑93 (MQ=255)

                                               gcgcGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGGCCATGTTATCCAATAACGCGTc                                          >  2:2668493/1‑93 (MQ=255)

                                                cgcgGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCa                                         >  2:8580255/1‑93 (MQ=255)

                                                cgcgGCCGTGCATCTTTACGCCACGGACCCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCa                                         <  2:2622127/93‑1 (MQ=255)

                                                 gcgGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCTAATTTTGTCCATGGAATCCAAAAAGGCGTCAt                                        >  2:6670510/1‑93 (MQ=255)

                                                 gcgGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCAt                                        <  2:4676980/93‑1 (MQ=255)

                                                   ggCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGc                                      >  1:4763273/1‑93 (MQ=255)

                                                    gCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCa                                     <  1:2942063/93‑1 (MQ=255)

                                                     ccGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAt                                    >  2:6179915/1‑93 (MQ=255)

                                                      cGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAtt                                   <  1:7014291/93‑1 (MQ=255)

                                                      cGTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAtt                                   <  1:8262792/93‑1 (MQ=255)

                                                       gTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAttt                                  <  2:272351/93‑1 (MQ=255)

                                                       gTGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAttt                                  >  2:3701385/1‑93 (MQ=255)

                                                        tGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAtttt                                 >  1:1659012/1‑93 (MQ=255)

                                                        tGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAtttt                                 >  2:1846536/1‑93 (MQ=255)

                                                        tGCATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAtttt                                 <  1:5868090/93‑1 (MQ=255)

                                                          cATCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgt                               >  1:7972445/1‑92 (MQ=255)

                                                            tCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgctg                             >  2:3502463/1‑93 (MQ=255)

                                                            tCTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgctg                             >  1:681560/1‑93 (MQ=255)

                                                             cTTTACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgctgc                            <  2:4487240/93‑1 (MQ=255)

                                                              tttACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCg                           >  1:3796897/1‑93 (MQ=255)

                                                                tACGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGaa                         <  1:8047736/93‑1 (MQ=255)

                                                                 aCGCCACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAAc                        <  2:8601335/93‑1 (MQ=255)

                                                                    ccACGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTc                     <  1:6464996/93‑1 (MQ=255)

                                                                      aCGGACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTg                   >  2:1464080/1‑93 (MQ=255)

                                                                        ggACGCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCa                 <  2:5934164/93‑1 (MQ=255)

                                                                           cccAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCaaaa              <  1:8137717/91‑1 (MQ=255)

                                                                            gCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAAc             >  1:6363776/1‑93 (MQ=255)

                                                                            gCAGACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAAc             <  2:7596783/93‑1 (MQ=255)

                                                                               gACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGGATTTTGCTGCGAACGTCTGCAAAACGGa          >  1:7753371/1‑93 (MQ=255)

                                                                               gACTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGa          <  1:5786767/93‑1 (MQ=255)

                                                                                aCTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACggct         >  2:255975/1‑91 (MQ=255)

                                                                                aCTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGat         <  2:3448585/93‑1 (MQ=255)

                                                                                aCTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGat         >  2:7850507/1‑93 (MQ=255)

                                                                                aCTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGat         <  1:945971/93‑1 (MQ=255)

                                                                                aCTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGat         >  1:4886367/1‑93 (MQ=255)

                                                                                 cTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGata        <  1:2388632/93‑1 (MQ=255)

                                                                                 cTCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGata        <  2:900040/93‑1 (MQ=255)

                                                                                  tCCATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGCACGTCTGCAAAACGGatat       >  1:6513222/1‑93 (MQ=255)

                                                                                   ccATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatc      <  2:5480917/93‑1 (MQ=255)

                                                                                   ccATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatc      <  2:7552798/93‑1 (MQ=255)

                                                                                    cATTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatca     >  1:3140701/1‑93 (MQ=255)

                                                                                     aTTGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatcat    >  2:6825142/1‑93 (MQ=255)

                                                                                      ttGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatcatc   <  2:4265297/93‑1 (MQ=255)

                                                                                       tGGCAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCg  <  2:4147034/93‑1 (MQ=255)

ACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCG  >  NC_000913/3854608‑3854787

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