BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,759,194	C→T	A208T (GCA→ACA)	ynhC ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,759,194	0	C	T	99.3%	473.6 / ‑5.1	135	A208T (GCA→ACA)	ynhC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (58/76); total (58/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

Only 100 of 142 total aligned reads displayed.

CTCAGCCAGCAACAAATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGCGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCACCTTCCG  >  NC_000913/1759105‑1759283

cTCAGCCAGCAACAAATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGAc                                                                                        >  2:2765649/1‑93 (MQ=255)

  cagccagcAACAAATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGt                                                                                      <  2:756637/93‑1 (MQ=255)

  cagccagcAACAAATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGt                                                                                      <  2:321882/93‑1 (MQ=255)

   agccagcAACAAATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGtt                                                                                     <  1:972866/93‑1 (MQ=255)

      cagcAACAAATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGTCATTCGCTGTGACGTTGAt                                                                                  <  1:3673225/93‑1 (MQ=255)

      cagcAACAAATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGAt                                                                                  >  1:1102920/1‑93 (MQ=255)

       agcAACAAATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATa                                                                                 >  2:748018/1‑93 (MQ=255)

       agcAACAAATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATa                                                                                 >  2:3480246/1‑93 (MQ=255)

        gcAACAAATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAg                                                                                <  1:735093/93‑1 (MQ=255)

        gcAACAAATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAg                                                                                <  2:4608771/93‑1 (MQ=255)

         caacaaATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGt                                                                               <  1:1673168/93‑1 (MQ=255)

          aacaaATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTg                                                                              <  1:1188877/93‑1 (MQ=255)

            caaATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGaa                                                                            >  2:896201/1‑93 (MQ=255)

            caaATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGaa                                                                            <  2:4372293/93‑1 (MQ=255)

              actCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAAcc                                                                          <  1:2647243/91‑1 (MQ=255)

                tCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGt                                                                        <  1:2540556/93‑1 (MQ=255)

                   ttATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGcc                                                                     <  1:541680/93‑1 (MQ=255)

                   ttATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGcc                                                                     <  1:3447622/93‑1 (MQ=255)

                    tATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGccc                                                                    <  2:247338/93‑1 (MQ=255)

                    tATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGccc                                                                    <  1:4477554/93‑1 (MQ=255)

                     aTGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGcccc                                                                   <  2:3432863/93‑1 (MQ=255)

                      tGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCg                                                                  >  2:2614181/1‑93 (MQ=255)

                       gAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCgg                                                                 <  1:1898447/93‑1 (MQ=255)

                         gCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTa                                                               <  2:57342/93‑1 (MQ=255)

                           aaaGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTG‑CCCGGTaaaa                                                            <  1:3034700/93‑1 (MQ=255)

                            aaGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTaaaa                                                            <  1:3303773/93‑1 (MQ=255)

                            aaGGGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTaaaa                                                            <  1:1074001/93‑1 (MQ=255)

                             aGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAAt                                                           <  2:4478755/93‑1 (MQ=255)

                              gtggtgACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGATCCGTGCCCCGGTAAAATg                                                          >  1:591929/1‑93 (MQ=255)

                               tggtgACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGa                                                         <  2:1710488/93‑1 (MQ=255)

                               tggtgACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGa                                                         <  1:2588878/93‑1 (MQ=255)

                                ggtgACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGAc                                                        <  2:463295/93‑1 (MQ=255)

                                  tgACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTACAATGACGa                                                      >  2:3256305/1‑93 (MQ=255)

                                   gACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAg                                                     <  1:3694587/93‑1 (MQ=255)

                                         gCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGcatcat                                               >  1:3210957/1‑93 (MQ=255)

                                         gCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGcatcat                                               >  2:3493676/1‑93 (MQ=255)

                                           gggTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGCGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTc                                             <  2:2153977/93‑1 (MQ=255)

                                             tctttCATACGCCAGCTTGATTTGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTATAATGACGAGCATCATTCAg                                           >  2:4196213/3‑93 (MQ=255)

                                             gTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAg                                           >  1:4560178/1‑93 (MQ=255)

                                              ttttCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAgg                                          >  1:861803/1‑93 (MQ=255)

                                               tttCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGc                                         <  1:2556836/93‑1 (MQ=255)

                                                 tCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTg                                       >  2:2335384/1‑93 (MQ=255)

                                                  cAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGa                                      <  2:4132682/93‑1 (MQ=255)

                                                   aaaCGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGAc                                     <  1:1426411/93‑1 (MQ=255)

                                                       gCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACaaaa                                 >  2:2781965/1‑93 (MQ=255)

                                                        ccAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAAt                                <  2:1920394/93‑1 (MQ=255)

                                                         cAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATg                               >  1:1228233/1‑93 (MQ=255)

                                                         cAGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATg                               >  1:675595/1‑93 (MQ=255)

                                                          aGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGt                              <  2:1954120/93‑1 (MQ=255)

                                                          aGCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGt                              >  1:2511595/1‑93 (MQ=255)

                                                           gCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGtt                             <  1:1011149/93‑1 (MQ=255)

                                                           gCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGtt                             <  2:995665/93‑1 (MQ=255)

                                                           gCTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGtt                             >  1:389974/1‑93 (MQ=255)

                                                            cTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTc                            <  1:4393936/93‑1 (MQ=255)

                                                            cTTGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTG‑CCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCg                           <  1:2133386/93‑1 (MQ=255)

                                                             ttGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCg                           <  1:1520410/93‑1 (MQ=255)

                                                             ttGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCg                           >  2:4302210/1‑93 (MQ=255)

                                                             ttGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCg                           >  1:3762112/1‑93 (MQ=255)

                                                             ttGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCg                           <  1:1520411/93‑1 (MQ=255)

                                                              tGATATTGTGCAAGTGGGCATTCGCCGTGACGTTGATAGTGAACCGTGCCCCTGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGa                          <  1:1808321/93‑1 (MQ=255)

                                                              tGATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGa                          <  1:2384503/93‑1 (MQ=255)

                                                               gATATGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGAt                         <  1:444222/93‑1 (MQ=255)

                                                                atatGCTGGAAGTGGTCATCCGCTGTAACGTTGATAGTGATCCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATc                        <  2:1677024/93‑1 (MQ=255)

                                                                atatGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATc                        >  1:487939/1‑93 (MQ=255)

                                                                 tatGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCa                       >  2:1803871/1‑93 (MQ=255)

                                                                 tatGCTGCAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCa                       >  2:882632/1‑93 (MQ=255)

                                                                  atGCTGCAAGTGGCCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCAc                      >  1:2792579/1‑93 (MQ=255)

                                                                   tgctgcAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCAcc                     >  1:3040876/1‑93 (MQ=255)

                                                                   tgctgcAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCAcc                     <  2:2648238/93‑1 (MQ=255)

                                                                   tgctgcAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCAcc                     <  1:2155531/93‑1 (MQ=255)

                                                                   tgctgcAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCAcc                     <  2:1037450/93‑1 (MQ=255)

                                                                    gctgcAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCg                    <  1:4594269/93‑1 (MQ=255)

                                                                     ctgcAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATTTTCGATCACCGt                   <  1:3833780/93‑1 (MQ=255)

                                                                     ctgcAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGt                   <  2:4017404/93‑1 (MQ=255)

                                                                     ctgcAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGt                   <  1:890531/93‑1 (MQ=255)

                                                                      tgcAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGtt                  <  2:3539712/93‑1 (MQ=255)

                                                                      tgcAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGtt                  <  2:4663627/93‑1 (MQ=255)

                                                                      tgcAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGtt                  <  1:2787096/93‑1 (MQ=255)

                                                                      tgcAAGTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGtt                  <  2:924877/93‑1 (MQ=255)

                                                                           gTGGGCATTCGCTGTGACGTTGATAGTGAACCGTTCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTc             <  1:581360/93‑1 (MQ=255)

                                                                           gTGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTc             >  2:2858731/1‑93 (MQ=255)

                                                                            tGGGCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCg            >  2:747991/1‑93 (MQ=255)

                                                                             gggCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCgg           >  1:4347722/1‑93 (MQ=255)

                                                                              ggCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGc          <  2:3909813/93‑1 (MQ=255)

                                                                              ggCATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGc          >  2:1332334/1‑93 (MQ=255)

                                                                                cATTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCAc        <  1:2003291/93‑1 (MQ=255)

                                                                                 aTTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCAcc       >  1:3101790/1‑93 (MQ=255)

                                                                                 aTTCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCAcc       <  2:3634408/93‑1 (MQ=255)

                                                                                  ttCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTCCGATCACCGTTGCTTCGGCACCt      <  2:4102481/93‑1 (MQ=255)

                                                                                   tCGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCACCtt     <  2:3578492/93‑1 (MQ=255)

                                                                                    cGCTGTTACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCACCTTc    <  1:407813/93‑1 (MQ=255)

                                                                                    cGCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCACCTTc    <  1:546353/93‑1 (MQ=255)

                                                                                    cGCTGTGACGTTGATAATGAACCGTGACCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCACCTTc    <  1:3526068/93‑1 (MQ=255)

                                                                                     gCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTTGATCACCGTTGCTTCGGCACCTTcc   >  1:244499/1‑93 (MQ=255)

                                                                                     gCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCACCTTcc   <  1:1708601/93‑1 (MQ=255)

                                                                                     gCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCACCTTcc   <  2:3683477/93‑1 (MQ=255)

                                                                                     gCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCACCTTcc   >  2:3358379/1‑93 (MQ=255)

                                                                                     gCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCACCTTcc   <  2:2110914/93‑1 (MQ=255)

                                                                                     gCTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCACCTTcc   <  1:60947/93‑1 (MQ=255)

                                                                                      cTGTGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCACCTTCCg  >  2:699907/1‑93 (MQ=255)

CTCAGCCAGCAACAAATCGTTATGAGCAAAGTGGTGACTGAGCGGGTTTTCAAACGCCAGCTTGATATGCTGCAAGTGGGCATTCGCTGCGACGTTGATAGTGAACCGTGCCCCGGTAAAATGACGAGCATCATTCAGGCTGACAAAATGTTCGATCACCGTTGCTTCGGCACCTTCCG  >  NC_000913/1759105‑1759283

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