BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,866,583	A→G	S119S (AGT→AGC)	yidR ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,866,583	0	A	G	99.3%	951.0 / ‑3.1	277	S119S (AGT→AGC)	yidR	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (135/140); total (137/140)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

Only 100 of 284 total aligned reads displayed.

TTCACCGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTACTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAA  >  NC_000913/3866493‑3866672

ttCACCGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCt                                                                                         >  2:6191111/1‑93 (MQ=255)

 tCACCGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTc                                                                                        >  2:2011403/1‑93 (MQ=255)

 tCACCGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTc                                                                                        <  1:2349216/93‑1 (MQ=255)

 tCACCGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTc                                                                                        >  1:2674009/1‑93 (MQ=255)

    ccGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATc                                                                                     >  2:5174950/1‑93 (MQ=255)

     cGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCt                                                                                    >  2:2969153/1‑93 (MQ=255)

     cGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCt                                                                                    <  1:1638449/93‑1 (MQ=255)

     cGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCt                                                                                    <  2:3696872/93‑1 (MQ=255)

      gTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCtt                                                                                   >  1:2920220/1‑93 (MQ=255)

       ttCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTg                                                                                  >  1:3933822/1‑93 (MQ=255)

       ttCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTg                                                                                  <  1:3900266/93‑1 (MQ=255)

          gggCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTgccg                                                                               >  1:596957/1‑93 (MQ=255)

           ggCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTgccgc                                                                              >  1:3306504/1‑93 (MQ=255)

           ggCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTgccgc                                                                              >  2:6237656/1‑93 (MQ=255)

            gCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTgccgcc                                                                             <  2:4181534/93‑1 (MQ=255)

             cTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCt                                                                            <  1:1201450/93‑1 (MQ=255)

             cTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCt                                                                            <  1:4569108/93‑1 (MQ=255)

               aaaGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTc                                                                          <  1:3625556/93‑1 (MQ=255)

                aaGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTca                                                                         >  2:3641338/1‑93 (MQ=255)

                aaGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTca                                                                         <  2:163571/93‑1 (MQ=255)

                aaGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTca                                                                         <  2:1345776/93‑1 (MQ=255)

                 aGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTcaa                                                                        >  2:509518/1‑93 (MQ=255)

                  gACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTcaac                                                                       <  2:2274072/93‑1 (MQ=255)

                    cATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTcaacaa                                                                     <  1:5767802/93‑1 (MQ=255)

                    cATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTcaacaa                                                                     <  2:5550779/93‑1 (MQ=255)

                      tGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATc                                                                   >  1:5379299/1‑93 (MQ=255)

                       gCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCa                                                                  >  1:2619294/1‑93 (MQ=255)

                        cacaTGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCAc                                                                 <  2:1826001/93‑1 (MQ=255)

                         acaTGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACt                                                                <  1:3054765/93‑1 (MQ=255)

                           aTGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTcc                                                              >  1:651195/1‑93 (MQ=255)

                             ggCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCcgc                                                            <  2:5005922/93‑1 (MQ=255)

                             ggCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCcgc                                                            >  2:6296920/1‑93 (MQ=255)

                                tGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGcc                                                         <  1:6449181/93‑1 (MQ=255)

                                 gccgccGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCg                                                        >  1:6392124/1‑93 (MQ=255)

                                 gccgccGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCg                                                        >  1:3316402/1‑93 (MQ=255)

                                  ccgccgCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCga                                                       >  2:5755327/1‑93 (MQ=255)

                                  ccgccgCGCAGCACTCCTGGCATATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCga                                                       >  1:1473079/1‑93 (MQ=255)

                                   cgccgcGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCgat                                                      <  2:3136903/93‑1 (MQ=255)

                                   cgccgcGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCgat                                                      <  2:1001464/93‑1 (MQ=255)

                                    gccgcGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCgatg                                                     >  1:1117445/1‑93 (MQ=255)

                                     ccgcGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCgatga                                                    <  2:2540632/93‑1 (MQ=255)

                                     ccgcGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCgatga                                                    >  1:1862043/1‑93 (MQ=255)

                                        cgcAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATgg                                                 <  2:3982412/93‑1 (MQ=255)

                                        cgcAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATgg                                                 <  2:1170888/93‑1 (MQ=255)

                                        cgcAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATgg                                                 >  1:2997468/1‑93 (MQ=255)

                                         gcagcaCTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGa                                                <  2:5465781/93‑1 (MQ=255)

                                         gcagcaCTCCTGGCGTATACGGCGCGCTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGa                                                <  1:6035220/93‑1 (MQ=255)

                                          cagcaCTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGaa                                               >  1:2848760/1‑93 (MQ=255)

                                          cagcaCTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGaa                                               <  1:5867357/93‑1 (MQ=255)

                                           agcaCTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAg                                              >  2:4211250/1‑93 (MQ=255)

                                              aCTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCg                                           >  1:3712122/1‑93 (MQ=255)

                                                  cTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAAt                                       <  1:1350540/93‑1 (MQ=255)

                                                  cTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGCCGTAAt                                       >  2:1786386/1‑93 (MQ=255)

                                                  cTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTGGCTCATCTTGCCGCCTTCAACAATCGCTCCGCGCCGATGATGGAAGTCGTAAt                                       <  1:2566297/93‑1 (MQ=255)

                                                   tGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATg                                      >  1:6272160/1‑93 (MQ=255)

                                                   tGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATg                                      >  2:3423872/1‑93 (MQ=255)

                                                    ggCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGc                                     <  2:4612520/93‑1 (MQ=255)

                                                     gCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGcc                                    <  2:4413162/93‑1 (MQ=255)

                                                      cGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCa                                   <  2:5923735/93‑1 (MQ=255)

                                                      cGTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCa                                   >  2:3294201/1‑93 (MQ=255)

                                                       gTATACGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCAt                                  >  1:4241621/1‑93 (MQ=255)

                                                        tataCGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATg                                 >  2:904309/1‑93 (MQ=255)

                                                        tataCGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATg                                 >  2:5586197/1‑93 (MQ=255)

                                                        tataCGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATg                                 <  1:4983850/93‑1 (MQ=255)

                                                        tataCGGCGCGGTAATATCCATCGCATCAAGGTTGCCCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATg                                 <  1:4376365/93‑1 (MQ=255)

                                                         ataCGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGt                                >  2:6453973/1‑93 (MQ=255)

                                                          taCGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGtt                               >  1:836294/1‑93 (MQ=255)

                                                          taCGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACACCGCGCCGATGATGAAAGTCGTAATGCCATGtt                               <  2:4389721/93‑1 (MQ=255)

                                                           aCGGCGCGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGttt                              >  2:4497338/1‑93 (MQ=255)

                                                              gcgcGGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTtcat                           >  1:2413477/1‑93 (MQ=255)

                                                                gcgGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTtcatca                         <  2:1281162/93‑1 (MQ=255)

                                                                 cgGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAg                        >  2:2644745/1‑93 (MQ=255)

                                                                 cgGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAg                        >  2:210362/1‑93 (MQ=255)

                                                                 cgGTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAg                        >  2:4701111/1‑93 (MQ=255)

                                                                   gTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGa                      >  1:5601735/1‑93 (MQ=255)

                                                                   gTAATATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGa                      >  2:726397/1‑93 (MQ=255)

                                                                     aaTATCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGAtt                    <  2:4854723/93‑1 (MQ=255)

                                                                      atatCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGAttt                   >  1:6161122/1‑93 (MQ=255)

                                                                      atatCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGAttt                   <  2:5984548/93‑1 (MQ=255)

                                                                      atatCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGAttt                   <  2:3247684/93‑1 (MQ=255)

                                                                         tCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTcc                <  2:2456470/93‑1 (MQ=255)

                                                                         tCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTcc                >  1:4312885/1‑93 (MQ=255)

                                                                         tCCATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTcc                <  2:1395622/93‑1 (MQ=255)

                                                                          ccATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCg               >  1:3828522/1‑93 (MQ=255)

                                                                           cATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCGCTCCCCGCCGATGATGAAAGTCGTAAGGCCATCTTTCATCAGGATTTTCCgc              >  1:5244056/1‑92 (MQ=255)

                                                                           cATCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCgg              >  1:1830816/1‑93 (MQ=255)

                                                                            aTCGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCggg             >  1:186931/1‑93 (MQ=255)

                                                                              cggaTCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGcc           >  2:5283558/4‑93 (MQ=255)

                                                                              cGCATCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGcc           <  1:709299/93‑1 (MQ=255)

                                                                               cgaTCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCg          <  2:1921196/91‑1 (MQ=255)

                                                                                 aTCAAGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTg        >  1:5127918/1‑93 (MQ=255)

                                                                                    aaGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGAt     <  1:4076568/93‑1 (MQ=255)

                                                                                    aaGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGAt     <  1:6113833/93‑1 (MQ=255)

                                                                                    aaGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGAt     <  1:4512432/93‑1 (MQ=255)

                                                                                     aGGTTGCTCATCTTGCCGCCTTCAACAATCACTCTGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATg    >  1:6361566/1‑93 (MQ=255)

                                                                                     aGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATg    >  1:2153843/1‑93 (MQ=255)

                                                                                     aGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATg    <  2:1668031/93‑1 (MQ=255)

                                                                                     aGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATg    >  2:6062407/1‑93 (MQ=255)

                                                                                     aGGTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATg    <  2:2782031/93‑1 (MQ=255)

                                                                                       gTTGCTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGaa  <  2:3384543/93‑1 (MQ=255)

TTCACCGTTCGGGCTAAAGACATGCACATGGCTGCCGCCGCGCAGCACTCCTGGCGTATACGGCGCGGTAATATCCATCGCATCAAGGTTACTCATCTTGCCGCCTTCAACAATCACTCCGCGCCGATGATGGAAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAA  >  NC_000913/3866493‑3866672

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