BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,544,598	A→G	15.1%	N99D (AAT→GAT)	fimC →	periplasmic chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,598	0	A	G	15.1%	432.5 / 64.3	190	N99D (AAT→GAT)	fimC	periplasmic chaperone
Reads supporting (aligned to +/- strand): ref base A (81/80); new base G (14/15); total (95/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.74e-01

Only 100 of 201 total aligned reads displayed.

CCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAAT  >  NC_000913/4544500‑4544697

ccTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGaa                                                                                                    <  1:3045784/100‑1 (MQ=255)

ccTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGaa                                                                                                    >  1:2736993/1‑100 (MQ=255)

  tATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATa                                                                                                  >  2:937445/1‑100 (MQ=255)

  tATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATa                                                                                                  <  2:1659107/100‑1 (MQ=255)

   aTTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATAc                                                                                                 >  2:1553985/1‑100 (MQ=255)

   aTTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATAc                                                                                                 >  2:2087575/1‑100 (MQ=255)

    tttAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATAcc                                                                                                <  1:3391028/100‑1 (MQ=255)

      tAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCtt                                                                                              >  2:3114029/1‑100 (MQ=255)

      tAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCtt                                                                                              >  1:988680/1‑100 (MQ=255)

       aaTTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAgagagatcggaa                                                                                             >  2:3274721/1‑92 (MQ=255)

        aTTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTAc                                                                                            >  2:2412794/1‑100 (MQ=255)

           cAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTa                                                                                         >  1:2337762/1‑100 (MQ=255)

             aTCATGGGTGGAAAGTGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTAtt                                                                                       <  1:1343830/100‑1 (MQ=255)

               cATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCt                                                                                     >  2:2212859/1‑100 (MQ=255)

               cATGGATGAAAAAAGCCGATGATGTAAAGAACGCTCGTTTTATTGAGACGCCTCCTCTGTTTGAAAAAAAAAAAAAAAAAGAGAATACCTTACGTATTCt                                                                                     <  2:3178714/100‑1 (MQ=255)

                atgaggagAAAATGCCGATGGTGTAAAAGAAGGTCGTTTTATCGTGACGCCTCCTCTGTTTAAAAAAAAAAGAAAAAAAGAGAATACCTTACGTATTCt                                                                                     <  2:196224/92‑1 (MQ=255)

                 tGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTg                                                                                   <  1:772859/100‑1 (MQ=255)

                  gggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACGTAAGGTATTCtctt                                                                                  <  1:490830/100‑4 (MQ=255)

                   ggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGAt                                                                                 >  1:2094089/1‑100 (MQ=255)

                   ggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGAt                                                                                 >  1:2903304/1‑100 (MQ=255)

                   ggTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGAt                                                                                 <  2:1128019/100‑1 (MQ=255)

                    gTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGCAGGGAAAAAAAGAGAATACCTTACGTATTCTTGa                                                                                  >  2:1695633/1‑98 (MQ=255)

                    gTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATg                                                                                <  1:1412133/100‑1 (MQ=255)

                    gTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATg                                                                                >  1:3931513/1‑100 (MQ=255)

                    gTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATg                                                                                >  2:3716955/1‑100 (MQ=255)

                    gTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATg                                                                                >  2:4134910/1‑100 (MQ=255)

                      ggAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACGTAAGGTATTCtctttttt                                                                              <  2:3002230/100‑8 (MQ=38)

                      ggAAAATGCCGATGGTGTAAAGGATGGTCGTGTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGc                                                                               <  1:3802256/99‑1 (MQ=255)

                       gAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGca                                                                              <  1:1973314/99‑1 (MQ=255)

                       gAAAATGCCGATGATGTAAAGGATGGTCGTTTTATCGTGAAGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACGTAAGGTATTCtctttttt                                                                              <  2:393855/99‑8 (MQ=37)

                        aaaaTGCCGATGGTGTAAAGGATGGTCGTTTTATCTTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACGTAAGGTATTCtctttttttc                                                                            <  1:2923698/100‑10 (MQ=37)

                        aaaaTGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGcaac                                                                            <  2:1074119/100‑1 (MQ=255)

                           aTGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACaaa                                                                         >  2:618571/1‑100 (MQ=255)

                           aTGCCGATGGAATAAAGGATGGTCGTTATATCGTGACGCCTCATCTGTTAGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGcaacaa                                                                          >  2:269007/1‑99 (MQ=255)

                               cGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAAc                                                                     >  1:1875548/1‑100 (MQ=255)

                               cGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCGTACGTATTCTTGATGCAACAAATAAc                                                                     <  2:4162365/100‑1 (MQ=255)

                               cGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTAAAAAAAAAAAAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAAc                                                                     <  2:1109250/100‑1 (MQ=255)

                                gATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAAAAAATAAc                                                                     >  2:1460527/1‑99 (MQ=255)

                                gATGGTGGAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCAATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAAcc                                                                    <  2:40629/100‑1 (MQ=255)

                                gATGGAGTAAAGGAAGGTCGTATTATCGTGACGCCTCCTTTGTTTGAAAAAAAAAAAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAAc                                                                     <  2:3836046/99‑1 (MQ=255)

                                  tGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCaa                                                                  <  1:690388/100‑1 (MQ=255)

                                   ggTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAAt                                                                 >  1:3835244/1‑100 (MQ=255)

                                    gtgtAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAAtt                                                                >  1:538329/1‑100 (MQ=255)

                                    gtgtAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAAtt                                                                >  1:3365944/1‑100 (MQ=255)

                                    gtgtAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTGGCGATGCAGGGAAAAATAGAGAATACCTTACGTATTCTTGATGCaacaaataacaaat                                                                 >  2:2918883/1‑99 (MQ=255)

                                    gtgtAAAGAATGGTCGTTTTATCGTGAAGCCTCCTTTGTTTGCGAAAAAAAAAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAAtt                                                                <  2:2110795/100‑1 (MQ=255)

                                     tgtAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTg                                                               >  2:87023/1‑100 (MQ=255)

                                     tgtAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTg                                                               >  2:2139330/1‑100 (MQ=255)

                                      gtAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGc                                                              >  1:3630693/1‑100 (MQ=255)

                                        aaaGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGcc                                                             >  1:3080560/1‑99 (MQ=255)

                                        aaaGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCca                                                            >  2:148929/1‑100 (MQ=255)

                                         aaGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCcac                                                           >  1:2160312/1‑100 (MQ=255)

                                          aGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCcaca                                                          <  1:1373263/100‑1 (MQ=255)

                                                gTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCg                                                    >  2:567650/1‑100 (MQ=255)

                                                gTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCg                                                    >  1:2399368/1‑100 (MQ=255)

                                                 tCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCgg                                                   <  1:3993923/100‑1 (MQ=255)

                                                 tCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCgg                                                   <  1:1135110/100‑1 (MQ=255)

                                                 tCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAAAGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCgg                                                   <  2:3599271/100‑1 (MQ=255)

                                                 tCGTTTTATCGTGACGCCTCCTCTGTTTGCGATAAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCgg                                                   <  2:694942/100‑1 (MQ=255)

                                                     tttATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGaaa                                               >  2:1266737/1‑100 (MQ=255)

                                                     tttATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGaaa                                               >  2:582265/1‑100 (MQ=255)

                                                      ttATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAg                                              <  1:591355/100‑1 (MQ=255)

                                                        aTCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGtt                                            >  2:2003085/1‑100 (MQ=255)

                                                         tCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGttt                                           <  1:1054032/100‑1 (MQ=255)

                                                           gTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAt                                         <  1:1383736/100‑1 (MQ=255)

                                                           gTGACGCCCCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAt                                         <  1:2525070/100‑1 (MQ=255)

                                                            tagcGCCTTCTCTGTTTGCGTTGAAGGTAAAAAAAGAAAAGACCTTACGTATTCGTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGtt                                            <  2:3305922/93‑1 (MQ=255)

                                                            tGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAtt                                        <  1:1749878/100‑1 (MQ=255)

                                                            tGACGCCTCCTCTGTGTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAtt                                        <  1:4115644/100‑1 (MQ=255)

                                                               cGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTg                                     >  2:3594131/1‑100 (MQ=255)

                                                               cGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTg                                     >  1:351125/1‑100 (MQ=255)

                                                                 cctcctCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGa                                   <  1:1360373/100‑1 (MQ=255)

                                                                 cctcctCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGa                                   >  1:2592993/1‑100 (MQ=255)

                                                                  ctcctcTGTTTGAGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTgg                                    <  2:1839159/98‑1 (MQ=255)

                                                                  cgcctgTGTTTGAGATGAAGAGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                                  <  1:3332399/98‑1 (MQ=255)

                                                                         gTTTGCGATGAAGGGAAAAAAAGAGAATACGTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTAGATGAACGtt                           >  2:2733408/1‑100 (MQ=255)

                                                                         gTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGtt                           >  1:562141/1‑100 (MQ=255)

                                                                          tttGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTa                          >  2:4009562/1‑100 (MQ=255)

                                                                              cGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGc                      <  1:2261570/100‑1 (MQ=255)

                                                                               gATGAAGGGAAAAAAAGAGGCTACCTTCCGCATTCTTGATGCACCAAATCAACAATTGCCCCAGGCCCGGGCAAGGTTATTCTGGATGAACGTTaaa                        <  2:1860511/97‑1 (MQ=255)

                                                                               gATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCg                     <  1:3905925/100‑1 (MQ=255)

                                                                                aTGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGAAGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGa                    >  2:2576602/1‑100 (MQ=255)

                                                                                  gAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGAtt                  <  1:1006025/100‑1 (MQ=255)

                                                                                     gggAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCg               >  1:1636219/1‑100 (MQ=255)

                                                                                      ggAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGt              >  1:1400550/1‑100 (MQ=255)

                                                                                       gAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTc             <  1:315433/100‑1 (MQ=255)

                                                                                       gAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTc             >  1:3265311/1‑100 (MQ=255)

                                                                                       gAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTc             >  2:722585/1‑100 (MQ=255)

                                                                                       gAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTc             <  1:1956681/100‑1 (MQ=255)

                                                                                        aaaaaaaGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCa            <  1:2646096/100‑1 (MQ=255)

                                                                                         aaaaaaGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCaa           >  2:3195270/1‑100 (MQ=255)

                                                                                          aaaaaGAGAATACCTTACGTATTCTTGATGCAACAAATAACAAATTGCCACAGGACCGGGAAAGTTTACTCTGGATGAACGTTAAAGCGATTCCGTAAAt          >  2:513710/1‑100 (MQ=255)

                                                                                           cataaagaATACGTAAGGTATTCTTGATGCAACAACTAAACAATTGCCACAGGACGAGGAAAGTTCATTCTGGATGAAAGTGAAAGCGATTCAGTCAATg         >  2:628486/6‑100 (MQ=255)

                                                                                           aaaaGAGAATACCTTACGTATTCTTGATGCAACAAATAACAAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATg         >  2:2753977/1‑100 (MQ=255)

                                                                                            aaaGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATgg        >  1:1663639/1‑100 (MQ=255)

                                                                                             tcaagaATACGTAAGGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGa       >  2:2030082/4‑100 (MQ=255)

                                                                                              caagaATACGTAAGGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGAt      >  2:2649232/3‑100 (MQ=255)

                                                                                              caagaATACGTAAGGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGAt      >  2:2304929/3‑100 (MQ=255)

                                                                                               gagaATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATa     <  1:1227915/100‑1 (MQ=255)

                                                                                                  aaTACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAAt  <  2:656855/100‑1 (MQ=255)

CCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAAT  >  NC_000913/4544500‑4544697

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