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	92.6%	557.5 / 36.8	190	coding (411/1071 nt)	b0941	FimH‑like protein
Reads supporting (aligned to +/- strand): ref base . (8/6); new base G (100/76); total (108/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 222 total aligned reads displayed.

TGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCC  >  NC_000913/1001348‑1001527

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

tGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGcaca                                                                                             <  2:6009669/93‑1 (MQ=255)

  gTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑g                                                                                       <  1:2425150/93‑1 (MQ=255)

  gTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑g                                                                                       <  2:461576/93‑1 (MQ=255)

  gTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑g                                                                                       <  2:2560024/93‑1 (MQ=255)

   tACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑gg                                                                                      <  1:6892831/93‑1 (MQ=255)

    aCAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑ggg                                                                                     >  1:6430761/1‑93 (MQ=255)

    aCAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑ggg                                                                                     <  2:3885970/93‑1 (MQ=255)

     cAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑gggg                                                                                    <  2:3046322/93‑1 (MQ=255)

      aaaaTTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑ggggg                                                                                   >  2:1797854/1‑93 (MQ=255)

      aaaaTTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑ggggg                                                                                   >  2:5519691/1‑93 (MQ=255)

      aaaaTTAATGACTATCTTGCCGCACAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTCCCTTTGGCAAACGGGGGCACA‑‑‑g‑ggggg                                                                                   >  1:5805302/1‑93 (MQ=255)

       aaaTTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑gggggg                                                                                  <  2:6095034/93‑1 (MQ=255)

        aaTTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑ggggggg                                                                                 >  1:1817365/1‑92 (MQ=255)

         aTTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑gggggggt                                                                                >  1:6245254/1‑91 (MQ=255)

         aTTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑ggggggga                                                                                <  1:4922164/93‑3 (MQ=255)

           tAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑g‑gggggggata                                                                              <  1:3559312/93‑5 (MQ=255)

                aCTATCTTGCCGCGAAAGTTTTATTGCAGGTTAAAGGCAGTACTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGa                                                                          <  2:6870752/93‑1 (MQ=255)

                 cTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGAt                                                                         <  2:6669876/93‑1 (MQ=255)

                  tATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATg                                                                        >  1:1329049/1‑93 (MQ=255)

                  tATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATg                                                                        >  1:4365276/1‑93 (MQ=255)

                      ttGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCat                                                                    >  1:7262611/1‑93 (MQ=255)

                       tGCCGCCAAAGTTTTATTGCAGTTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGGTACCCGATGGCata                                                                   <  1:1862840/93‑1 (MQ=255)

                          cgccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGG‑G‑GGGGGGATACCCGATGGCATAtt                                                                 <  2:7589153/93‑1 (MQ=255)

                          cgccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGGTACCCGATGGCATAttt                                                                >  1:4340903/1‑93 (MQ=255)

                          cgccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATAttt                                                                <  1:3230482/93‑1 (MQ=255)

                           gccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTg                                                               <  2:6185735/93‑1 (MQ=255)

                           gccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTg                                                               <  1:4780861/93‑1 (MQ=255)

                            ccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGc                                                              <  1:2753108/93‑1 (MQ=255)

                            ccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGc                                                              <  2:7827409/93‑1 (MQ=255)

                            ccAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGc                                                              >  2:4990098/1‑93 (MQ=255)

                             cAAAGTTTTATTGCAGGGTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATAtttttg                                                             >  2:2748317/1‑90 (MQ=255)

                              aaaGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGa                                                            <  2:5360411/93‑1 (MQ=255)

                              aaaGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCCGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGa                                                            >  1:4801747/1‑93 (MQ=255)

                               aaGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGAc                                                           <  1:6478403/93‑1 (MQ=255)

                                aGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGAcc                                                          >  1:1166216/1‑93 (MQ=255)

                                 gTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATCTTTGAGAccc                                                         >  2:903418/1‑93 (MQ=255)

                                 gTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGAccc                                                         <  2:5733810/93‑1 (MQ=255)

                                  ttttATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCt                                                        <  2:6561136/93‑1 (MQ=255)

                                   tttATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTg                                                       >  2:538780/1‑93 (MQ=255)

                                   tttATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTg                                                       >  2:4787695/1‑93 (MQ=255)

                                    ttATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGG‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTg                                                       >  1:3862301/1‑93 (MQ=255)

                                         gCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTc                                                 <  1:1604530/93‑1 (MQ=255)

                                         gCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTc                                                 <  2:3628240/93‑1 (MQ=255)

                                            ggTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAAc                                              >  1:5098943/1‑93 (MQ=255)

                                              ttAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGt                                            >  2:7674092/1‑93 (MQ=255)

                                              ttAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑GTGGGGGGGTACCAGATCGCATATTTGCGACCCCGCCGCTCagctt                                            >  1:5937468/1‑89 (MQ=255)

                                                aaaGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTa                                         >  2:6468603/1‑93 (MQ=255)

                                                 aaGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTa                                         >  2:38235/1‑93 (MQ=255)

                                                  aGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAg                                        <  1:3189305/93‑1 (MQ=255)

                                                   ggCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGG‑GGGGGGGTACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTa                                         >  1:7435792/1‑93 (MQ=255)

                                                   ggCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAgg                                       <  1:6205506/93‑1 (MQ=255)

                                                     cAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTg                                     <  1:194449/93‑1 (MQ=255)

                                                      aGTTCTCCTACTTCGGTTCCTTTCGTAAACGTGGGCACAG‑‑G‑GGGGGGATACTCGATGGCTTATTTGCGACCCTGGCGGTCATCGCTTAGGTgg                                    >  2:7035529/1‑93 (MQ=255)

                                                      aGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTgg                                    <  2:6784477/93‑1 (MQ=255)

                                                      aGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTgg                                    <  2:7309679/93‑1 (MQ=255)

                                                      aGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTgg                                    <  2:3259557/93‑1 (MQ=255)

                                                       gTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTggg                                   >  1:2287624/1‑93 (MQ=255)

                                                       gTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTggg                                   >  2:4442317/1‑93 (MQ=255)

                                                       gTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTggg                                   >  1:5711268/1‑93 (MQ=255)

                                                       gTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTggg                                   >  1:5250722/1‑93 (MQ=255)

                                                       gTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTggg                                   >  2:205044/1‑93 (MQ=255)

                                                        ttCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGc                                  >  2:318072/1‑93 (MQ=255)

                                                        ttCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGGGGGcg                                 >  2:7388350/1‑92 (MQ=255)

                                                           tcCTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAgg                               >  1:234620/1‑93 (MQ=255)

                                                            ccTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAggg                              >  2:5228323/1‑93 (MQ=255)

                                                            ccTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGCTTAAGGGGGCAggg                              >  2:1487928/1‑93 (MQ=255)

                                                             cTACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAgggg                             >  1:4173497/1‑93 (MQ=255)

                                                              tACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAggggg                            >  2:3770027/1‑93 (MQ=255)

                                                              tACTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAggggg                            >  1:5554732/1‑93 (MQ=255)

                                                               tcTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGc                           <  1:3436828/92‑1 (MQ=255)

                                                               aCTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGc                           >  2:2615116/1‑93 (MQ=255)

                                                                cTGCGGTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCa                          >  2:1981812/1‑93 (MQ=255)

                                                                 tGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAg                        >  1:384054/1‑93 (MQ=255)

                                                                    ggTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCg                      >  2:6307700/1‑93 (MQ=255)

                                                                    ggTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCg                      <  2:3450849/93‑1 (MQ=255)

                                                                     gTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCgg                     >  2:3694202/1‑93 (MQ=255)

                                                                     gTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCgg                     <  1:5771540/93‑1 (MQ=255)

                                                                     gTTCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCgg                     >  1:7540407/1‑93 (MQ=255)

                                                                       tCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTa                   <  2:5065036/93‑1 (MQ=255)

                                                                       tCCTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTa                   >  1:653300/1‑93 (MQ=255)

                                                                         cTTTCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTaaa                 >  2:4574660/1‑92 (MQ=255)

                                                                          tttCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATa                >  1:6494045/1‑93 (MQ=255)

                                                                           ttCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAg               >  1:4163875/1‑93 (MQ=255)

                                                                            tCGAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGc              <  2:1632342/93‑1 (MQ=255)

                                                                            tCGAAAACGTGGGCACA‑‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCg             <  2:2132718/93‑1 (MQ=255)

                                                                              gAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCgg            <  2:4022871/93‑1 (MQ=255)

                                                                              gAAAACGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCgg            >  1:4439629/1‑93 (MQ=255)

                                                                               aaaaCGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGt           <  1:988265/93‑1 (MQ=255)

                                                                                aaaCGTGGGCACAG‑‑G‑GGGGGGATCCCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTa          <  2:4359561/93‑1 (MQ=255)

                                                                                aaaCGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTa          <  1:3602931/93‑1 (MQ=255)

                                                                                 aaCGTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGAGGTAg         >  1:468005/1‑93 (MQ=255)

                                                                                   cGTGGGCACA‑‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCtt      >  1:3374832/1‑93 (MQ=255)

                                                                                    gTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCtt      >  1:78704/1‑93 (MQ=255)

                                                                                    gTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCtt      >  1:6164348/1‑93 (MQ=255)

                                                                                    gTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCtt      >  1:358184/1‑93 (MQ=255)

                                                                                    gTGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCtt      <  1:1207780/93‑1 (MQ=255)

                                                                                     tGGGCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCttt     >  1:2685288/1‑93 (MQ=255)

                                                                                       ggCACAG‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTc   <  2:4097026/93‑1 (MQ=255)

                                                                                        gCACAG‑‑G‑GGGGGGATACCCGATGGCTTATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTcc  >  1:5776259/1‑93 (MQ=255)

TGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACA‑‑‑G‑GGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCC  >  NC_000913/1001348‑1001527

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