BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,546,975	T→C	V200A (GTG→GCG)	fimH →	minor fimbrial subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,546,975	0	T	C	93.8%	1207.7 / 40.5	423	V200A (GTG→GCG)	fimH	minor fimbrial subunit
Reads supporting (aligned to +/- strand): ref base T (16/10); new base C (224/173); total (240/183)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

Only 100 of 440 total aligned reads displayed.

CGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTT  >  NC_000913/4546886‑4547064

 gCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCa                                                                                       >  1:7859584/1‑93 (MQ=255)

  ccAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCaa                                                                                      <  1:476377/93‑1 (MQ=255)

   cAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAAt                                                                                     <  2:541355/93‑1 (MQ=255)

   cAATAATGAGGTGGTGGTGCCTACTGCCGGCTGAGAAGTTTCTGCCCGTGATGTCAGCGTTACGCTGCCGGAGTAACCGTGTTCAgcgcacat                                                                                     >  2:1856173/1‑89 (MQ=255)

    aataatGATGTGGTGGTGCCTACTGGCGGCTGCGCTGTTTCTGCTCGTGATGTCACCGTTACCCTGCCGGACTACCCCGTTTCAGCGCAAAtt                                                                                    >  1:9859381/1‑93 (MQ=255)

    aataatGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGAAGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAAtt                                                                                    >  1:5353133/1‑93 (MQ=255)

     ataatGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTc                                                                                   <  1:1011241/93‑1 (MQ=255)

     ataatGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTc                                                                                   >  2:813835/1‑93 (MQ=255)

      taatGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTcc                                                                                  >  1:8127204/1‑93 (MQ=255)

        atgatgTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCctc                                                                                >  2:7396517/1‑93 (MQ=255)

         tgatgTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCctct                                                                               <  2:6649738/93‑1 (MQ=255)

         tgatgTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCctct                                                                               >  1:4467566/1‑93 (MQ=255)

           atgTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTCCTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTa                                                                             <  2:9657506/93‑1 (MQ=255)

               ggtggtgCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGt                                                                         >  1:7170620/1‑93 (MQ=255)

                 tggtgCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGttt                                                                       <  1:2578986/93‑1 (MQ=255)

                 tggtgCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGttt                                                                       >  2:8676725/1‑93 (MQ=255)

                  ggtgCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTa                                                                      <  2:2659499/93‑1 (MQ=255)

                     gCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATtg                                                                   >  1:1081090/1‑93 (MQ=255)

                     gCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATtg                                                                   <  2:4904881/93‑1 (MQ=255)

                       cTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATtgtg                                                                 >  2:1722766/1‑93 (MQ=255)

                        tACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGc                                                                <  2:4629355/93‑1 (MQ=255)

                         aCTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCg                                                               >  1:3597445/1‑93 (MQ=255)

                          cTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGa                                                              >  1:6235583/1‑93 (MQ=255)

                          cTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGa                                                              >  2:5274707/1‑93 (MQ=255)

                           tGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGaa                                                             >  1:8906072/1‑93 (MQ=255)

                           tGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGaa                                                             <  2:4101404/93‑1 (MQ=255)

                           tGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGaa                                                             >  1:10012426/1‑93 (MQ=255)

                            ggcggcTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGaaa                                                            <  1:2292402/93‑1 (MQ=255)

                            ggcggcTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGaaa                                                            >  2:8895902/1‑93 (MQ=255)

                               ggcTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGc                                                         <  1:5086682/93‑1 (MQ=255)

                               ggcTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGc                                                         >  1:9066274/1‑93 (MQ=255)

                                gcTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGcc                                                        >  1:3177692/1‑93 (MQ=255)

                                gcTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGcc                                                        >  2:485347/1‑93 (MQ=255)

                                gcTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGcc                                                        <  1:2611737/93‑1 (MQ=255)

                                gcTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGcc                                                        >  1:5129051/1‑93 (MQ=255)

                                 cTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCa                                                       <  2:9018977/93‑1 (MQ=255)

                                 cTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAACCCa                                                       >  1:3076845/1‑93 (MQ=255)

                                  tGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTCCCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCaa                                                      <  1:7432851/93‑1 (MQ=255)

                                  tGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCaa                                                      >  1:6890750/1‑93 (MQ=255)

                                  tGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCaa                                                      >  2:374108/1‑93 (MQ=255)

                                  tGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCaa                                                      >  1:8015102/1‑93 (MQ=255)

                                     gATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAAc                                                   >  2:8003370/1‑93 (MQ=255)

                                     gATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAAc                                                   <  2:8478664/93‑1 (MQ=255)

                                     gATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAAc                                                   <  2:4442828/93‑1 (MQ=255)

                                     gATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAAc                                                   >  2:1184110/1‑93 (MQ=255)

                                     gATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAAc                                                   >  1:3486121/1‑93 (MQ=255)

                                      aTGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAAcc                                                  >  1:5724658/1‑93 (MQ=255)

                                      aTGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAAcc                                                  >  2:3081910/1‑93 (MQ=255)

                                       tGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCt                                                 <  1:9639720/93‑1 (MQ=255)

                                       tGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCt                                                 >  2:3862487/1‑93 (MQ=255)

                                       tGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCt                                                 >  2:2851173/1‑93 (MQ=255)

                                       tGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCt                                                 <  2:1031124/93‑1 (MQ=255)

                                        gTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTg                                                >  1:5529368/1‑93 (MQ=255)

                                          ttCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTggg                                              <  2:2255546/93‑1 (MQ=255)

                                           tCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTgggg                                             <  2:2506893/93‑1 (MQ=255)

                                           tCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTgggg                                             >  2:5701848/1‑93 (MQ=255)

                                            cTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGt                                            <  2:1454685/93‑1 (MQ=255)

                                             tGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTa                                           <  1:24014/93‑1 (MQ=255)

                                                    gATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCCAAACCTGGGGTATTACctc                                    >  2:7655297/1‑93 (MQ=255)

                                                     aTGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCATACCGTTTATTGTGCGAAAAGCCAAAACCTGGGCTATTACctct                                   >  1:8401297/1‑93 (MQ=255)

                                                      tGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCtctc                                  <  1:9917504/93‑1 (MQ=255)

                                                       gTCACCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTcc                                 <  2:8290067/93‑1 (MQ=255)

                                                          aCCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGc                              >  1:2213706/1‑93 (MQ=255)

                                                          aCCGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGc                              >  1:4169076/1‑93 (MQ=255)

                                                           ccGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGca                             >  1:3987633/1‑93 (MQ=255)

                                                           ccGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGca                             <  2:3976450/93‑1 (MQ=255)

                                                            cGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGcac                            <  2:2755874/93‑1 (MQ=255)

                                                            cGTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGcac                            <  1:1184602/93‑1 (MQ=255)

                                                             gTTACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGccca                           >  2:6687894/1‑93 (MQ=255)

                                                               tACTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAAc                         >  2:9109761/1‑93 (MQ=255)

                                                                aCTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAAcc                        >  1:7661284/1‑93 (MQ=255)

                                                                aCTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAAcc                        >  2:9469666/1‑93 (MQ=255)

                                                                aCTCTGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAAcc                        >  1:8588371/1‑93 (MQ=255)

                                                                 ctctGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCg                       >  1:9482593/1‑93 (MQ=255)

                                                                  tctGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGc                      >  1:6531005/1‑93 (MQ=255)

                                                                   ctGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCa                     <  2:5898567/93‑1 (MQ=255)

                                                                   ctGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCa                     <  1:32803/93‑1 (MQ=255)

                                                                   ctGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCa                     <  2:4220642/93‑1 (MQ=255)

                                                                    tGCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGAATTACCTCTCCGGCACAACCGCAg                    >  2:7499579/1‑93 (MQ=255)

                                                                     gCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGa                   <  1:9825244/93‑1 (MQ=255)

                                                                     gCCGGACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGa                   <  2:2409328/93‑1 (MQ=255)

                                                                     gCCGGACTACCCTGGTGCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGa                   <  2:3237854/93‑1 (MQ=255)

                                                                      ccGGGCTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAACACCTGGGGTATTACCTCTCCGGCACAACCGCAGAt                  <  1:1790965/93‑1 (MQ=255)

                                                                      ccGGACTGCCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGAt                  <  2:7251242/93‑1 (MQ=255)

                                                                      ccGGACTACACTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGAt                  <  2:9279015/93‑1 (MQ=255)

                                                                        ggACTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGc                <  1:207081/93‑1 (MQ=255)

                                                                          aCTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCgg              >  2:6353528/1‑93 (MQ=255)

                                                                          aCTACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGCATTACCTCTCCGGCACAACCGCAGATGCgg              >  2:8775817/1‑93 (MQ=255)

                                                                           cttcccTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCggg             <  2:3433946/90‑1 (MQ=255)

                                                                            tACCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGc            >  2:4501637/1‑93 (MQ=255)

                                                                             aCCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCa           >  2:7141766/1‑93 (MQ=255)

                                                                             aCCCTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCa           >  1:5843081/1‑93 (MQ=255)

                                                                              cccTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCaa          <  1:5271859/93‑1 (MQ=255)

                                                                                cTGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACt        >  1:3751283/1‑93 (MQ=255)

                                                                                 tGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGTTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTc       <  1:1362194/93‑1 (MQ=255)

                                                                                 tGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTc       <  1:8752425/93‑1 (MQ=255)

                                                                                 tGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTc       >  1:6394927/1‑93 (MQ=255)

                                                                                 tGGTTCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTc       <  2:6111743/93‑1 (MQ=255)

                                                                                    ttCAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGAt    >  1:4644133/1‑93 (MQ=255)

                                                                                      cAGCGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGAttt  >  2:3627077/1‑93 (MQ=255)

CGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTT  >  NC_000913/4546886‑4547064

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