BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,477,321	Δ1 bp	coding (12/1182 nt)	yjgN →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,477,321	0	A	.	99.3%	686.5 / ‑3.5	151	coding (12/1182 nt)	yjgN	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base . (106/44); minor base C (0/1); total (106/45)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 151 total aligned reads displayed.

AACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATAGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGATTATAACG  >  NC_000913/4477232‑4477408

aaCATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTcc                                                                                    >  1:4497103/1‑93 (MQ=255)

 aCATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCg                                                                                   <  1:3308268/93‑1 (MQ=255)

 aCATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCg                                                                                   >  1:2164940/1‑93 (MQ=255)

  cATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGt                                                                                  >  2:4243728/1‑93 (MQ=255)

   aTTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTc                                                                                 >  2:3931695/1‑93 (MQ=255)

        aCTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCAtt                                                                            <  2:153485/93‑1 (MQ=255)

          tgtctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCg                                                                          >  1:3399788/1‑93 (MQ=255)

          tgtctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCg                                                                          >  1:4008752/1‑93 (MQ=255)

          tgtctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCg                                                                          >  2:4040055/1‑93 (MQ=255)

           gtctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgt                                                                         >  2:533399/1‑93 (MQ=255)

            tctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtt                                                                        >  1:3448963/1‑93 (MQ=255)

            tctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtt                                                                        <  2:834821/93‑1 (MQ=255)

            tctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtt                                                                        <  2:2180265/93‑1 (MQ=255)

             ctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgttt                                                                       >  1:2923801/1‑93 (MQ=255)

              tgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttg                                                                      <  1:1733575/93‑1 (MQ=255)

              tgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttg                                                                      >  1:4247175/1‑93 (MQ=255)

              tgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttg                                                                      <  2:1127590/93‑1 (MQ=255)

                tatatataAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttgtt                                                                    >  2:3035776/1‑93 (MQ=255)

                tatatataAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttgtt                                                                    >  2:3923633/1‑93 (MQ=255)

                tatatataAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttgtt                                                                    >  1:1129438/1‑93 (MQ=255)

                 atatataAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtctgttt                                                                   >  1:767671/1‑93 (MQ=255)

                  tatataAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGtttt                                                                  >  2:2091019/1‑93 (MQ=255)

                   atataAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGttttt                                                                 >  2:1872013/1‑93 (MQ=255)

                    tataAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTc                                                                <  1:1529620/93‑1 (MQ=255)

                    tataAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTc                                                                >  1:4376735/1‑93 (MQ=255)

                     ataAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCa                                                               >  2:620636/1‑93 (MQ=255)

                     ataAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCa                                                               >  2:4085441/1‑93 (MQ=255)

                      taAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCAt                                                              >  2:361162/1‑93 (MQ=255)

                      taAGTTATACCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCAt                                                              >  2:3642897/1‑93 (MQ=255)

                       aaGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATg                                                             >  2:931246/1‑93 (MQ=255)

                        aGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATgg                                                            <  2:3494238/93‑1 (MQ=255)

                         gTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGt                                                           >  1:3829305/1‑93 (MQ=255)

                           tataTCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTAc                                                         >  2:1291576/1‑93 (MQ=255)

                            atatCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACa                                                        >  1:1986259/1‑93 (MQ=255)

                             tatCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAg                                                       >  2:1132365/1‑93 (MQ=255)

                             tatCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAg                                                       >  2:1787003/1‑93 (MQ=255)

                               tCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAggg                                                     >  1:3330409/1‑92 (MQ=255)

                               tCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGt                                                     >  1:3187971/1‑93 (MQ=255)

                                cAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTg                                                    >  1:4303481/1‑93 (MQ=255)

                                cAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTg                                                    >  2:3779121/1‑93 (MQ=255)

                                cAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTg                                                    >  2:2875177/1‑93 (MQ=255)

                                 aaTGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgg                                                   >  2:3731882/1‑92 (MQ=255)

                                  aTGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgag                                                  >  1:299203/1‑93 (MQ=255)

                                  aTGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgag                                                  <  1:1249837/93‑1 (MQ=255)

                                   tGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgaga                                                 <  1:3553308/93‑1 (MQ=255)

                                     gATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgagaga                                               >  2:4147912/1‑93 (MQ=255)

                                      aTTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGat                                              <  2:4122961/93‑1 (MQ=255)

                                       ttttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGata                                             >  1:915479/1‑93 (MQ=255)

                                       ttttttAA‑TTCAGGGAATTCTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGata                                             >  2:4547802/1‑93 (MQ=255)

                                        tttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGatat                                            >  1:3453149/1‑93 (MQ=255)

                                        tttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGatat                                            <  2:4329539/93‑1 (MQ=255)

                                             taattCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTc                                        <  2:279985/90‑1 (MQ=255)

                                             taattCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTc                                        <  2:2106209/90‑1 (MQ=255)

                                             taattCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTc                                        >  2:1615272/4‑93 (MQ=255)

                                              aattCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCt                                       >  1:3980603/3‑93 (MQ=255)

                                               attCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCtt                                      >  1:1414049/2‑93 (MQ=255)

                                               attCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCtt                                      <  1:3351852/92‑1 (MQ=255)

                                                ttCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTa                                     >  1:2912800/1‑93 (MQ=255)

                                                  cAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTAtt                                   <  1:1062303/93‑1 (MQ=255)

                                                   aGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTAttt                                  >  1:1878977/1‑93 (MQ=255)

                                                   aGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTAttt                                  >  2:2612270/1‑93 (MQ=255)

                                                   aGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGA‑CGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTAttt                                  <  2:784996/93‑1 (MQ=255)

                                                      gAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTtgtg                               >  1:3561898/1‑93 (MQ=255)

                                                       aaTTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTtgtgt                              <  2:1082278/93‑1 (MQ=255)

                                                         tttttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTgtgg                            >  1:481952/1‑93 (MQ=255)

                                                         tttttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTgtgg                            >  2:3675494/1‑93 (MQ=255)

                                                          ttttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTgtggt                           <  2:2183530/93‑1 (MQ=255)

                                                            ttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGa                         >  2:3424852/1‑93 (MQ=255)

                                                            ttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGa                         >  1:2081918/1‑93 (MQ=255)

                                                             tATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGaa                        <  2:1577842/93‑1 (MQ=255)

                                                             tATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGaa                        <  1:1999449/93‑1 (MQ=255)

                                                              aTGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAAt                       <  2:188648/93‑1 (MQ=255)

                                                               tGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAAtg                      <  1:3148375/93‑1 (MQ=255)

                                                               tGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAAtg                      >  2:1950892/1‑93 (MQ=255)

                                                               tGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAAtg                      >  1:2710915/1‑93 (MQ=255)

                                                                 gCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAAtgtg                    >  2:1713299/1‑93 (MQ=255)

                                                                  cTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATgtgt                   <  1:3014715/93‑1 (MQ=255)

                                                                  cTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATgtgt                   <  1:1967724/93‑1 (MQ=255)

                                                                    cAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGtgttg                 >  2:2684042/1‑93 (MQ=255)

                                                                     aaGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGtgttgt                >  1:697024/1‑93 (MQ=255)

                                                                     aaGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGtgttgt                <  2:616614/93‑1 (MQ=255)

                                                                      aGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTgttgtt               <  1:965477/93‑1 (MQ=255)

                                                                      aGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTgttgtt               <  1:1651344/93‑1 (MQ=255)

                                                                        ttattaATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTaa             <  1:2170599/93‑1 (MQ=255)

                                                                         tattaATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAAc            >  1:4511946/1‑93 (MQ=255)

                                                                          attaATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACg           >  2:334591/1‑93 (MQ=255)

                                                                          attaATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACg           >  1:1211539/1‑93 (MQ=255)

                                                                          attaATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACg           >  1:4426510/1‑93 (MQ=255)

                                                                           ttaATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGa          >  1:2814836/1‑93 (MQ=255)

                                                                           ttaATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGa          >  1:885636/1‑93 (MQ=255)

                                                                             aaTGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGAtt        >  1:904805/1‑93 (MQ=255)

                                                                             aaTGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGAtt        >  1:4163942/1‑93 (MQ=255)

                                                                             aaTGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGAtt        >  2:1761882/1‑93 (MQ=255)

                                                                              aTGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGATta       <  1:565456/93‑1 (MQ=255)

                                                                               tGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGATtat      >  1:1574244/1‑93 (MQ=255)

                                                                               tGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGATtat      >  1:3791806/1‑93 (MQ=255)

                                                                                  aaTGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGATTATAAc   >  1:1569459/1‑93 (MQ=255)

                                                                                  aaTGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGATTATAAc   <  1:3067767/93‑1 (MQ=255)

                                                                                   aTGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGATTATAACg  >  1:3918419/1‑93 (MQ=255)

                                                                                   aTGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGATTATAACg  >  2:3109206/1‑93 (MQ=255)

AACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATAGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGTTGTTAACGATTATAACG  >  NC_000913/4477232‑4477408

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