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.6%	1042.1 / ‑3.6	226	intergenic (+43/‑25)	b4316/fimD	Pil/FimD
Reads supporting (aligned to +/- strand): ref base G (1/0); new base . (110/115); total (111/115)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 256 total aligned reads displayed.

TAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATC  >  NC_000913/4542548‑4542728

tAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCGGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCg                                                                                          >  2:5172089/1‑93 (MQ=255)

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

tAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCg                                                                                          >  2:2165026/1‑93 (MQ=255)

 aaaTGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgg                                                                                         <  2:4848608/93‑1 (MQ=255)

 aaaTGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgg                                                                                         >  2:5351582/1‑93 (MQ=255)

  aaTGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCggg                                                                                        <  1:2057689/93‑1 (MQ=255)

   aTGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgggt                                                                                       <  2:700606/93‑2 (MQ=255)

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

    tGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgggtg                                                                                      >  1:3837485/1‑91 (MQ=255)

    tGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgggtg                                                                                      <  1:4635671/93‑3 (MQ=255)

    tGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgggtg                                                                                      <  2:599417/93‑3 (MQ=255)

    tGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgggtg                                                                                      <  2:2559348/93‑3 (MQ=255)

      aTTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGTAATAACGCAGGGGGAATTTTTCGCTTGAATAAAAAGAATTGACTGCCgggtgat                                                                                    <  1:3104441/93‑5 (MQ=255)

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

       ttATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAtt                                                                                  >  2:6535801/1‑93 (MQ=255)

       ttATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAtt                                                                                  >  2:4390251/1‑93 (MQ=255)

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

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

           ggCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTaa                                                                              >  2:2131546/1‑93 (MQ=255)

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

              gcACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGcc                                                                           >  2:5786456/1‑93 (MQ=255)

              gcACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGcc                                                                           >  1:1715571/1‑93 (MQ=255)

               cACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCg                                                                          <  1:4754442/93‑1 (MQ=255)

               cACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCg                                                                          >  2:1262339/1‑93 (MQ=255)

               cACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCg                                                                          >  1:6379610/1‑93 (MQ=255)

                aCTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgg                                                                         <  2:4161046/93‑1 (MQ=255)

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

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

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

                 cTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgga                                                                        >  1:349561/1‑93 (MQ=255)

                  ttACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggag                                                                       <  1:583616/93‑1 (MQ=255)

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

                    aCCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggagga                                                                     >  1:2605622/1‑93 (MQ=255)

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

                     cccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaa                                                                    >  2:3956349/1‑93 (MQ=255)

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

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

                         cAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaataat                                                                <  1:6331300/93‑1 (MQ=255)

                          aaaaTGACGGGCGTAATGGAATAACGCGGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATg                                                               <  1:202002/93‑1 (MQ=255)

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

                            aaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGCATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTc                                                             >  2:970415/1‑93 (MQ=255)

                            aaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTc                                                             >  1:944916/1‑93 (MQ=255)

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

                                aCGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCatat                                                         >  2:237694/1‑93 (MQ=255)

                                     cGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGaa                                                    >  2:4634639/1‑93 (MQ=255)

                                     cGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGaa                                                    <  1:5142827/93‑1 (MQ=255)

                                      gTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAt                                                   <  1:1077442/93‑1 (MQ=255)

                                      gTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAt                                                   <  2:6266734/93‑1 (MQ=255)

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

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

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

                                          tGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAAGGTCATATCTGAATTTaa                                               <  1:6232735/93‑1 (MQ=255)

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

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

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

                                               tAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACtt                                          <  2:3857493/93‑1 (MQ=255)

                                                 aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTa                                        <  1:1225280/93‑1 (MQ=255)

                                                 aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTa                                        <  2:1051491/93‑1 (MQ=255)

                                                  cGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAc                                       <  2:5834142/93‑1 (MQ=255)

                                                   gCAGGGGGAATCTTTCGCCTGAATAAAAAGAAGTGACTGCC‑GGTTGATTTTAAGCCGGAGTAATAATGTCATATCTGAATTTAAGACTTTAcc                                      <  1:4565781/93‑1 (MQ=255)

                                                    cAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCa                                     >  1:3192531/1‑93 (MQ=255)

                                                    cAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCa                                     <  2:2996769/93‑1 (MQ=255)

                                                    cAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCa                                     >  1:6338156/1‑93 (MQ=255)

                                                     aGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAg                                    >  1:1393100/1‑93 (MQ=255)

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

                                                      gggggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAACCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGc                                   <  2:3564629/93‑1 (MQ=255)

                                                       ggggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCg                                  >  2:855352/1‑93 (MQ=255)

                                                        gggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGa                                 >  2:3681251/1‑93 (MQ=255)

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

                                                            aTTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAaca                             <  2:4773650/93‑1 (MQ=255)

                                                              ttttCGCCTGAATAAAAAGAATTGACTGCC‑AGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacaca                           <  2:1862237/93‑1 (MQ=255)

                                                               tttCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacacac                          <  2:3778610/93‑1 (MQ=255)

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

                                                                  cGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAAt                       <  2:3240559/93‑1 (MQ=255)

                                                                  cGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAAt                       <  2:6729501/93‑1 (MQ=255)

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

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

                                                                    ccTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATCTCTGAATTTAAGACTTTACCAGCGAAACACACAATgt                     >  2:920130/1‑92 (MQ=255)

                                                                    ccTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGc                     <  2:2897327/93‑1 (MQ=255)

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

                                                                        aaTAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATACCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGc                 <  1:5438147/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                            aaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCatat             <  1:4537091/93‑1 (MQ=255)

                                                                             gaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATAtt            >  1:6438692/2‑93 (MQ=255)

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

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

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

                                                                                  aTTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTaa       <  1:560822/93‑1 (MQ=255)

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

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

                                                                                    tGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGc     >  1:2219072/1‑93 (MQ=255)

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

                                                                                      aCTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCAt   <  1:2479840/93‑1 (MQ=255)

                                                                                       cTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATc  >  1:5718068/1‑93 (MQ=255)

TAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATC  >  NC_000913/4542548‑4542728

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