BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,542,643	Δ1 bp	intergenic (+46/‑22)	b4316 → / → fimD	Pil/FimD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,542,640	0	G	.	99.1%	981.3 / ‑0.6	217	intergenic (+43/‑25)	b4316/fimD	Pil/FimD
Reads supporting (aligned to +/- strand): ref base G (2/0); new base . (114/101); total (116/101)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 240 total aligned reads displayed.

TAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATC  >  NC_000913/4542548‑4542728

tAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCg                                                                                          <  1:3989113/93‑1 (MQ=255)

 aaaTGATTATGGCGCACTTACACCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgg                                                                                         <  2:4659618/93‑1 (MQ=255)

   aTGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgggt                                                                                       >  1:569995/1‑92 (MQ=255)

   aTGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgggt                                                                                       >  1:3529178/1‑92 (MQ=255)

      aTTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgggtgat                                                                                    >  1:4290913/1‑89 (MQ=255)

      aTTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgggggat                                                                                    >  2:1637288/1‑90 (MQ=255)

       ttATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAtt                                                                                  <  1:2416525/93‑1 (MQ=255)

       ttATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAtt                                                                                  <  1:2246521/93‑1 (MQ=255)

          tGGTGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTa                                                                               <  1:544711/93‑1 (MQ=255)

          tGGCGCACTTGCCCCCAAAATGACGGGCGTAATGGAATAGCGCAGGGGGAATTTTTCGCTTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTa                                                                               <  2:501616/93‑1 (MQ=255)

           ggCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTaa                                                                              <  2:3375544/93‑1 (MQ=255)

            gcgcACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAg                                                                             <  1:21871/93‑1 (MQ=255)

             cgcACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGc                                                                            <  1:4510511/93‑1 (MQ=255)

             cgcACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGc                                                                            >  2:1086961/1‑93 (MQ=255)

             cgcACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGc                                                                            <  2:88339/93‑1 (MQ=255)

              gcACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGcc                                                                           <  1:4200047/93‑1 (MQ=255)

               cACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGG‑GATTTTAAGCCg                                                                          >  2:2842956/1‑93 (MQ=255)

                aCTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgg                                                                         >  1:763013/1‑93 (MQ=255)

                 cTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAACAAGAATTGACTGCC‑GGGTGATTTTAAGCCgga                                                                        >  2:2509005/1‑93 (MQ=255)

                 cTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgga                                                                        <  2:139173/93‑1 (MQ=255)

                 cTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCATAAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgga                                                                        <  2:4477523/93‑1 (MQ=255)

                   tACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggagg                                                                      >  1:3003671/1‑93 (MQ=255)

                    aCCCCCAAAATGACGGGCGTATTGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggagga                                                                     >  1:346988/1‑93 (MQ=255)

                    aCCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggagga                                                                     <  1:194270/93‑1 (MQ=255)

                     cccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaa                                                                    <  2:2238802/93‑1 (MQ=255)

                      ccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaat                                                                   <  1:432594/93‑1 (MQ=255)

                      ccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaat                                                                   >  2:1141380/1‑93 (MQ=255)

                      ccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaat                                                                   >  1:550601/1‑93 (MQ=255)

                       cccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaata                                                                  >  1:1658142/1‑93 (MQ=255)

                       cccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaata                                                                  <  2:3015971/93‑1 (MQ=255)

                        ccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaataa                                                                 <  2:614876/93‑1 (MQ=255)

                         cAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaataat                                                                >  1:3666123/1‑93 (MQ=255)

                          aaaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATg                                                               <  2:810204/93‑1 (MQ=255)

                           aaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGt                                                              >  2:2166261/1‑93 (MQ=255)

                           aaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGt                                                              >  2:3504552/1‑93 (MQ=255)

                           aaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGt                                                              <  1:3791014/93‑1 (MQ=255)

                             aTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGG‑GATTTTAAGCCGGAGGAATATTGtaa                                                            >  1:1612422/1‑91 (MQ=255)

                             aTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCa                                                            <  1:307246/93‑1 (MQ=255)

                               gACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCata                                                          <  1:900881/93‑1 (MQ=255)

                               gACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCata                                                          >  1:2445616/1‑93 (MQ=255)

                                 cGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATc                                                        <  1:321785/93‑1 (MQ=255)

                                  gggCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCt                                                       <  1:2939332/93‑1 (MQ=255)

                                   ggCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTg                                                      <  2:4176955/93‑1 (MQ=255)

                                   ggCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTg                                                      >  1:1949254/1‑93 (MQ=255)

                                   ggCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTg                                                      >  2:4438460/1‑93 (MQ=255)

                                    gCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGa                                                     <  1:2852561/93‑1 (MQ=255)

                                      gTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAt                                                   >  2:884500/1‑93 (MQ=255)

                                       tAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAtt                                                  >  2:1367818/1‑93 (MQ=255)

                                       tAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAtt                                                  >  1:1003005/1‑93 (MQ=255)

                                       tAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAAGAACGTCATATCTGAAtt                                                  >  1:4208780/1‑93 (MQ=255)

                                        aaTGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAttt                                                 >  1:2069076/1‑93 (MQ=255)

                                        aaTGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTAACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAtta                                                 >  1:4146214/1‑92 (MQ=255)

                                         aTGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGAC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTa                                                >  1:3689274/1‑93 (MQ=255)

                                          tGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTaa                                               >  2:2019564/1‑93 (MQ=255)

                                          tGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTaa                                               <  2:3452498/93‑1 (MQ=255)

                                           ggAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAg                                              <  1:2814517/93‑1 (MQ=255)

                                           ggAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAg                                              >  1:3506761/1‑93 (MQ=255)

                                             aaTAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGAc                                            >  2:1458151/1‑93 (MQ=255)

                                              aTAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACt                                           <  2:4323245/93‑1 (MQ=255)

                                              aTAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACt                                           >  1:2896671/1‑93 (MQ=255)

                                                aaCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACttt                                         <  1:1437254/93‑1 (MQ=255)

                                                   gCAGGGGCAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAcc                                      <  1:4168123/93‑1 (MQ=255)

                                                      gggggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGc                                   >  1:3368566/1‑93 (MQ=255)

                                                      gggggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGc                                   >  2:1755233/1‑93 (MQ=255)

                                                      gggggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGc                                   >  1:4403660/1‑93 (MQ=255)

                                                          gAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGaaa                               >  1:4437097/1‑93 (MQ=255)

                                                           aaTTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAac                              <  2:4253589/93‑1 (MQ=255)

                                                           aaTTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAac                              <  2:4019263/93‑1 (MQ=255)

                                                            aTTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAaca                             >  1:2681973/1‑93 (MQ=255)

                                                              ttttCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacaca                           <  1:1861633/93‑1 (MQ=255)

                                                                 tCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACaa                        <  2:2427003/93‑1 (MQ=255)

                                                                   gCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATg                      >  1:1140898/1‑93 (MQ=255)

                                                                    ccTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGc                     >  2:1817966/1‑93 (MQ=255)

                                                                     cTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCt                    >  2:3661264/1‑93 (MQ=255)

                                                                     cTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCt                    <  1:3869357/93‑1 (MQ=255)

                                                                     cTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCt                    >  1:3746060/1‑93 (MQ=255)

                                                                     cTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAAAGCt                    >  1:602212/1‑93 (MQ=255)

                                                                      tGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCtt                   >  1:1553518/1‑93 (MQ=255)

                                                                       gAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTg                  >  1:1372557/1‑93 (MQ=255)

                                                                         aTAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCa                <  2:1506724/93‑1 (MQ=255)

                                                                          tAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCat               >  2:3119344/1‑93 (MQ=255)

                                                                          tAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCat               >  2:3799230/1‑93 (MQ=255)

                                                                           aaaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCata              >  2:3247809/1‑93 (MQ=255)

                                                                           aaaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCata              <  1:3734446/93‑1 (MQ=255)

                                                                           aaaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCata              >  2:1995253/1‑93 (MQ=255)

                                                                            aaaaGAATTGACTGCCGGGG‑GATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCatat             >  2:3280716/1‑93 (MQ=255)

                                                                            aaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCatat             >  1:2287992/1‑93 (MQ=255)

                                                                             aaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATAtt            <  2:4677827/93‑1 (MQ=255)

                                                                             aaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATAtt            >  1:1027895/1‑93 (MQ=255)

                                                                               aGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCg          >  2:3646519/1‑93 (MQ=255)

                                                                                gAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGt         >  1:1262392/1‑93 (MQ=255)

                                                                                 aaTTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTa        <  1:1953487/93‑1 (MQ=255)

                                                                                 aaTTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTa        <  1:4093420/93‑1 (MQ=255)

                                                                                   ttGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAg      >  1:938385/1‑93 (MQ=255)

                                                                                    tGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGc     <  1:2087427/93‑1 (MQ=255)

                                                                                     gACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCa    >  1:644189/1‑93 (MQ=255)

                                                                                     gACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCa    <  2:4598278/93‑1 (MQ=255)

                                                                                      aCTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACCCACAATGCTTGCATATTCGTAAGCAt   <  2:2284221/93‑1 (MQ=255)

                                                                                       cTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATc  <  2:1762722/93‑1 (MQ=255)

                                                                                       cTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATc  <  1:3308346/93‑1 (MQ=255)

TAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATC  >  NC_000913/4542548‑4542728

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