BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,001,447	(G)_7→8	coding (418/1071 nt)	b0941 →	FimH‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,001,440	1	.	G	98.2%	333.7 / ‑0.7	108	coding (411/1071 nt)	b0941	FimH‑like protein
Reads supporting (aligned to +/- strand): ref base . (1/1); new base G (55/51); total (56/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 128 total aligned reads displayed.

TGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑GGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCC  >  NC_000913/1001348‑1001528

tGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGcaca                                                                                             >  2:2830012/1‑93 (MQ=255)

tGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGcaca                                                                                             >  2:635442/1‑93 (MQ=255)

tGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGcaca                                                                                             >  2:735405/1‑93 (MQ=255)

 ggTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑g                                                                                           <  1:938295/93‑1 (MQ=255)

 ggTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑g                                                                                           <  1:3057070/93‑1 (MQ=255)

  gTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑gg                                                                                          >  1:269417/1‑93 (MQ=255)

  gTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑gg                                                                                          >  1:778003/1‑93 (MQ=255)

    aCAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑gggg                                                                                        <  1:1321164/93‑1 (MQ=255)

      aaaaTTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑gggg‑g‑g                                                                                    <  1:3156169/93‑1 (MQ=255)

        aaTTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑gggg‑g‑ggg                                                                                  >  1:3541754/1‑92 (MQ=255)

        aaTTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑gggg‑g‑ggg                                                                                  >  1:609973/1‑92 (MQ=255)

             aTGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATAcc                                                                              <  1:2600300/93‑1 (MQ=255)

             aTGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATAcc                                                                              <  1:3931727/93‑1 (MQ=255)

             aTGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATAcc                                                                              <  2:1601551/93‑1 (MQ=255)

                 cTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGAt                                                                          >  2:3503970/1‑93 (MQ=255)

                  tATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATg                                                                         >  1:3918251/1‑93 (MQ=255)

                  tATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATg                                                                         >  2:388864/1‑93 (MQ=255)

                   aTCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATgg                                                                        >  1:2371852/1‑93 (MQ=255)

                   aTCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATgg                                                                        >  2:3444415/1‑93 (MQ=255)

                   aTCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGCTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATgg                                                                        <  1:3296226/93‑1 (MQ=255)

                    tCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGc                                                                       >  2:2576384/1‑93 (MQ=255)

                      ttGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCat                                                                     >  1:27496/1‑93 (MQ=255)

                       tGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCata                                                                    <  2:409360/93‑1 (MQ=255)

                        gccgccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCatat                                                                   <  2:2608699/93‑1 (MQ=255)

                          cgccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATAttt                                                                 >  2:2103042/1‑93 (MQ=255)

                          cgccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATAttt                                                                 <  1:3837643/93‑1 (MQ=255)

                            ccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGc                                                               <  2:2458809/93‑1 (MQ=255)

                              aaaGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGa                                                             <  2:2129653/93‑1 (MQ=255)

                              aaaGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGa                                                             >  2:1642623/1‑93 (MQ=255)

                               aatttttATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCAGAGGGGG‑G‑GGATACCCGATGGCATATTTGCGAc                                                            <  2:4170552/90‑1 (MQ=255)

                               aaGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCCCA‑GGGGAG‑GGATAACCGATGGCATATTTGCGaa                                                            >  2:3377943/1‑92 (MQ=255)

                               aaGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGTGAc                                                            >  1:3733936/1‑93 (MQ=255)

                                 gTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGAccc                                                          >  1:2161780/1‑93 (MQ=255)

                                 gTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGAccc                                                          <  1:4089006/93‑1 (MQ=255)

                                  ttttATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCt                                                         <  1:476443/93‑1 (MQ=255)

                                   tttATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTg                                                        <  1:3535316/93‑1 (MQ=255)

                                    ttATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTGTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTgg                                                       <  2:493154/93‑1 (MQ=255)

                                     tATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGc                                                      >  1:1948841/1‑93 (MQ=255)

                                     tATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGc                                                      >  2:178075/1‑93 (MQ=255)

                                          cAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCa                                                 >  1:76319/1‑93 (MQ=255)

                                          cAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCa                                                 >  2:1679156/1‑93 (MQ=255)

                                           aGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCaa                                                <  2:4062587/93‑1 (MQ=255)

                                            ggTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAAc                                               >  2:4031167/1‑93 (MQ=255)

                                            ggTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAAc                                               <  2:2886497/93‑1 (MQ=255)

                                            ggTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAAc                                               >  2:1312584/1‑93 (MQ=255)

                                            ggTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAAc                                               <  1:1090405/93‑1 (MQ=255)

                                            ggTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAAc                                               >  1:4089935/1‑93 (MQ=255)

                                             gTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑GGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGt                                             <  1:4373991/93‑1 (MQ=255)

                                              ttAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGt                                             <  2:4643728/93‑1 (MQ=255)

                                              ttAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGt                                             >  1:3318478/1‑93 (MQ=255)

                                               tAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGtt                                            <  1:2475540/93‑1 (MQ=255)

                                                aaaGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGttt                                           <  1:3467608/93‑1 (MQ=255)

                                                 aaGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCAAAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTa                                          >  1:2365578/1‑93 (MQ=255)

                                                     cAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTg                                      <  2:3765315/93‑1 (MQ=255)

                                                     cAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTg                                      <  2:2360273/93‑1 (MQ=255)

                                                     cAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTg                                      <  1:1472730/93‑1 (MQ=255)

                                                     cAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTg                                      <  1:1219223/93‑1 (MQ=255)

                                                      aGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTgg                                     <  2:2269107/93‑1 (MQ=255)

                                                      aGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTgg                                     <  1:2040016/93‑1 (MQ=255)

                                                        ttCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGc                                   <  1:1741759/93‑1 (MQ=255)

                                                           tcCTACTGCGGTTTCTGTCGAAAACGTGGGTACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCTgg                                >  2:4686696/1‑93 (MQ=255)

                                                           tcCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAgg                                >  1:834434/1‑93 (MQ=255)

                                                              tACTGCTGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAggggg                             <  1:231487/93‑1 (MQ=255)

                                                               aCTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGc                            >  1:2234857/1‑93 (MQ=255)

                                                               aCTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGc                            >  2:4405273/1‑93 (MQ=255)

                                                               aCTGCGGTTCCTTTCGAAAACGTGGGCACA‑GGGG‑G‑GG‑TACCCGATGGGATATTTTTGACCCTCGCGGTCTACCTTTGTGGGGGGAggggggca                          >  2:540139/1‑90 (MQ=255)

                                                                cTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCa                           >  1:1526946/1‑93 (MQ=255)

                                                                 tGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCag                          >  1:1422418/1‑92 (MQ=255)

                                                                 tGCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCaa                          >  1:4407784/1‑93 (MQ=255)

                                                                  gCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAg                         >  1:3116297/1‑93 (MQ=255)

                                                                  gCGGTTCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAg                         <  1:1575613/93‑1 (MQ=255)

                                                                  gCGGTTCCTTTCGAAAACGTGGGAACAGGGGG‑G‑GGAAACCCGATGGCATATTTGCGACCCTGGCGGCCAACGTTTAGGTGGGCAGGGGGCAAg                         >  1:1879801/1‑93 (MQ=255)

                                                                      ttCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGt                     >  1:912184/1‑93 (MQ=255)

                                                                      ttCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGt                     >  1:266997/1‑93 (MQ=255)

                                                                       tCCTTTCGAAAACGTGGGCAGAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTa                    <  1:895308/93‑1 (MQ=255)

                                                                       tCCTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTa                    >  2:4260620/1‑93 (MQ=255)

                                                                        ccTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTaa                   <  2:823066/93‑1 (MQ=255)

                                                                        ccTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTaa                   <  1:4598580/93‑1 (MQ=255)

                                                                        ccTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTaa                   <  2:2654108/93‑1 (MQ=255)

                                                                         cTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTact                  >  2:3250107/1‑91 (MQ=255)

                                                                         cTTTCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAAt                  <  2:2133514/93‑1 (MQ=255)

                                                                          tttCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATa                 <  2:2919421/93‑1 (MQ=255)

                                                                           ttCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAg                <  2:3109086/93‑1 (MQ=255)

                                                                           ttCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAg                <  2:1398820/93‑1 (MQ=255)

                                                                           ttCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGCAATAg                >  2:3400933/1‑93 (MQ=255)

                                                                            tCGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGc               >  2:3429186/1‑93 (MQ=255)

                                                                             cGAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCg              <  1:471033/93‑1 (MQ=255)

                                                                              gAAAACGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCgg             >  1:3446717/1‑93 (MQ=255)

                                                                               aaaaCGTGGGCACAGGGGG‑G‑GGATACTCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGt            <  2:3795163/93‑1 (MQ=255)

                                                                               aaaaCGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGt            <  2:408904/93‑1 (MQ=255)

                                                                               aaaaCGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGt            >  2:229809/1‑93 (MQ=255)

                                                                               aaaaCGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGt            <  1:3729102/93‑1 (MQ=255)

                                                                               aaaaCGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGt            >  1:1324379/1‑93 (MQ=255)

                                                                                aaaCGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTa           >  2:3378144/1‑93 (MQ=255)

                                                                                   cGTGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCt        <  2:1120466/93‑1 (MQ=255)

                                                                                     tGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCttt      >  1:489384/1‑93 (MQ=255)

                                                                                     tGGGCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCttt      <  1:3067911/93‑1 (MQ=255)

                                                                                        gCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTcc   <  2:2141606/93‑1 (MQ=255)

                                                                                        gCACAGGGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTcc   <  1:3283582/93‑1 (MQ=255)

                                                                                         caca‑GGGG‑GAGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTccc  <  2:172853/93‑1 (MQ=255)

TGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑GGGG‑G‑GGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCC  >  NC_000913/1001348‑1001528

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