BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,422,211	T→C	A53A (GCA→GCG)	yjfY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,422,211	0	T	C	99.6%	721.4 / ‑5.1	231	A53A (GCA→GCG)	yjfY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (114/116); total (115/116)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

Only 100 of 236 total aligned reads displayed.

ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTAAT  >  NC_000913/4422122‑4422300

 tatTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCtt                                                                                       >  1:4002648/1‑93 (MQ=255)

   ttGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCg                                                                                     >  2:1218986/1‑93 (MQ=255)

   ttGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCg                                                                                     >  2:2434520/1‑93 (MQ=255)

    tGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGc                                                                                    <  1:602416/93‑1 (MQ=255)

     gCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCt                                                                                   >  2:361047/1‑93 (MQ=255)

      ccGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGcgcttcgctt                                                                                  >  1:3461176/1‑93 (MQ=255)

       cGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCttt                                                                                 <  2:4352127/93‑1 (MQ=255)

       cGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCttt                                                                                 <  2:3322119/93‑1 (MQ=255)

        gTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTc                                                                                >  1:4519310/1‑93 (MQ=255)

        gTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTc                                                                                <  2:2790322/93‑1 (MQ=255)

         ttACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCa                                                                               >  1:1335141/1‑93 (MQ=255)

         ttACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCa                                                                               >  1:3814581/1‑93 (MQ=255)

         ttACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCa                                                                               <  2:314067/93‑1 (MQ=255)

          tACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAg                                                                              >  1:1833138/1‑93 (MQ=255)

           tcTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGt                                                                             <  1:3931293/92‑1 (MQ=255)

           aCTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGt                                                                             <  1:726012/93‑1 (MQ=255)

           aCTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGt                                                                             >  2:567148/1‑93 (MQ=255)

           aCTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGt                                                                             >  1:2092822/1‑93 (MQ=255)

           aCTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGt                                                                             <  2:514521/93‑1 (MQ=255)

                ccGTTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGa                                                                        <  2:78950/93‑1 (MQ=255)

                  ggTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGaaa                                                                      <  1:536981/93‑1 (MQ=255)

                  ggTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGaaa                                                                      <  2:2079706/93‑1 (MQ=255)

                   gTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAAt                                                                     <  2:4164488/93‑1 (MQ=255)

                    ttCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAAtt                                                                    >  2:2739485/1‑93 (MQ=255)

                    ttCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAAtt                                                                    >  2:46638/1‑93 (MQ=255)

                     tCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTa                                                                   <  1:1580190/93‑1 (MQ=255)

                      cccGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAt                                                                  <  2:2400306/93‑1 (MQ=255)

                      cccGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAt                                                                  <  1:4572575/93‑1 (MQ=255)

                      cccGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAt                                                                  >  2:1128412/1‑93 (MQ=255)

                      cccGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAt                                                                  <  2:2823185/93‑1 (MQ=255)

                       ccGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtg                                                                 <  2:3387406/93‑1 (MQ=255)

                       ccGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtg                                                                 <  2:160208/93‑1 (MQ=255)

                       ccGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtg                                                                 >  1:703153/1‑93 (MQ=255)

                        cGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtga                                                                <  2:1045252/93‑1 (MQ=255)

                         gCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgat                                                               >  2:4054973/1‑93 (MQ=255)

                          cATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgatt                                                              <  2:3449606/93‑1 (MQ=255)

                           aTCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattg                                                             >  1:3964724/1‑93 (MQ=255)

                           aTCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattg                                                             >  1:1548982/1‑93 (MQ=255)

                            tCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattga                                                            >  1:2671431/1‑93 (MQ=255)

                            tCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattga                                                            >  1:3012624/1‑93 (MQ=255)

                               gAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGatat                                                         >  1:3942786/1‑93 (MQ=255)

                                aaTTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGAtata                                                        <  1:1231993/93‑1 (MQ=255)

                                 aTTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATaa                                                       >  2:2351983/1‑93 (MQ=255)

                                   tACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt                                                     <  1:4116995/93‑1 (MQ=255)

                                    aCGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATc                                                    >  1:3415593/1‑93 (MQ=255)

                                    aCGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATc                                                    <  1:1443139/93‑1 (MQ=255)

                                     cGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCa                                                   <  1:778734/93‑1 (MQ=255)

                                       tGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACg                                                 >  2:386679/1‑93 (MQ=255)

                                        ggTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGc                                                <  1:371544/93‑1 (MQ=255)

                                            gtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTaa                                            >  1:300381/1‑93 (MQ=255)

                                              aCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGc                                          <  1:947313/93‑1 (MQ=255)

                                              aCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGc                                          <  1:4022122/93‑1 (MQ=255)

                                                gTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCtt                                        >  1:2189776/1‑93 (MQ=255)

                                                gTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCATAAGCtt                                        >  1:3562407/1‑93 (MQ=255)

                                                 ttGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCttt                                       <  2:2529497/93‑1 (MQ=255)

                                                   gcgcCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGc                                     <  2:1641795/93‑1 (MQ=255)

                                                      ccTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGcact                                  >  1:1489334/1‑92 (MQ=255)

                                                      ccTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGcaca                                  <  1:1176314/93‑1 (MQ=255)

                                                      ccTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGcaca                                  <  2:442810/93‑1 (MQ=255)

                                                        ttGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATc                                >  1:1000823/1‑93 (MQ=255)

                                                        ttGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATc                                >  1:889361/1‑93 (MQ=255)

                                                          gcgcATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCTGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCgc                              >  2:1454863/1‑92 (MQ=255)

                                                          gcgcATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGt                              <  2:3842265/93‑1 (MQ=255)

                                                           cgcATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTc                             <  2:3967299/93‑1 (MQ=255)

                                                           cgcATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTc                             <  1:2407393/93‑1 (MQ=255)

                                                             cATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCa                           >  1:1705095/1‑93 (MQ=255)

                                                             cATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCa                           >  2:4551503/1‑93 (MQ=255)

                                                              aTCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCat                          >  2:1634866/1‑93 (MQ=255)

                                                              aTCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCat                          >  1:2761149/1‑93 (MQ=255)

                                                              aTCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCat                          >  1:4160086/1‑93 (MQ=255)

                                                              aTCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCat                          >  1:22651/1‑93 (MQ=255)

                                                              aTCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCat                          >  2:1819039/1‑93 (MQ=255)

                                                                cTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATGACGCCTAAGCTTTGCACATCGTCCatat                        >  2:1979161/1‑93 (MQ=255)

                                                                cTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCatat                        >  1:2846834/1‑93 (MQ=255)

                                                                 tGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATAtt                       >  2:4284233/1‑93 (MQ=255)

                                                                  gTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATAttt                      <  1:2155520/93‑1 (MQ=255)

                                                                   tttCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTc                     >  1:3829631/1‑93 (MQ=255)

                                                                    ttcttcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCt                    <  1:3213327/93‑1 (MQ=255)

                                                                     tcttcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTg                   >  2:30094/1‑93 (MQ=255)

                                                                       ttcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGc                 >  2:1727849/1‑93 (MQ=255)

                                                                       ttcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGc                 <  2:3187904/93‑1 (MQ=255)

                                                                       ttcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGc                 >  1:4262061/1‑93 (MQ=255)

                                                                        tcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCTATATTTCTGGcc                <  1:2532530/93‑1 (MQ=255)

                                                                         cATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCt               >  2:2702028/1‑93 (MQ=255)

                                                                         cATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCt               <  1:3877324/93‑1 (MQ=255)

                                                                          aTTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtg              >  2:2908195/1‑93 (MQ=255)

                                                                           tttAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtgg             >  2:65675/1‑93 (MQ=255)

                                                                           tttAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtgg             >  2:1084097/1‑93 (MQ=255)

                                                                            ttAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtggt            >  1:4625748/1‑93 (MQ=255)

                                                                            ttAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtggt            >  1:3659951/1‑93 (MQ=255)

                                                                            ttAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtggt            <  2:2931876/93‑1 (MQ=255)

                                                                             tAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtggtg           <  2:3766201/93‑1 (MQ=255)

                                                                              aaTGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtggtgg          >  2:3504836/1‑93 (MQ=255)

                                                                               aTGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTggtggt         <  1:3681532/93‑1 (MQ=255)

                                                                                tGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGtt        >  2:3420484/1‑93 (MQ=255)

                                                                                   cTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGttatt     <  2:4223987/93‑1 (MQ=255)

                                                                                    tGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGttatta    >  2:1952496/1‑93 (MQ=255)

                                                                                    tGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGttatta    <  1:3251299/93‑1 (MQ=255)

                                                                                     gACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTaa   <  1:2893315/93‑1 (MQ=255)

                                                                                      aCGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTAAt  >  1:2871390/1‑93 (MQ=255)

ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTAAT  >  NC_000913/4422122‑4422300

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