BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,544,598	A→G	47.2%	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	47.2%	25.2 / 253.4	184	N99D (AAT→GAT)	fimC	periplasmic chaperone
Reads supporting (aligned to +/- strand): ref base A (48/49); new base G (45/42); total (93/91)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 187 total aligned reads displayed.

ACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGAT  >  NC_000913/4544499‑4544693

aCCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAgaga                                                                                                 <  2:2167254/100‑1 (MQ=255)

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

  cTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAAt                                                                                               >  2:554709/1‑100 (MQ=255)

   tATTTAATTTAATCATGGGTGGAAAATGCCGATGGTGTAAAAGATGGTCGTTTTATAGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGAt                                                                                               >  1:3614020/1‑99 (MQ=255)

   tATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATa                                                                                              >  2:1745571/1‑100 (MQ=255)

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

        aactCAATCATGGGTGGAAAATGCAGATGGTGTAAAGGATGGTCGTTTTATAGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTa                                                                                         <  2:649945/97‑1 (MQ=255)

        aaTTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTa                                                                                         >  2:3284638/1‑100 (MQ=255)

         aTTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTAc                                                                                        <  2:3823990/100‑1 (MQ=255)

         aTTCAATAATGGGTGAAAAAAGCCGATAGTGTAAAGGATAGACATTATATCGTGACGCCTCCTCTGATTAAAAAAAAAAAAAAAAAAGAGAATACCTTAc                                                                                        <  1:3542531/100‑1 (MQ=255)

          ttCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACg                                                                                       >  2:2065719/1‑100 (MQ=255)

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

             aagcATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTAGTTTTATCGTGACGCCTTCTTTGTTTGAGAAAAAAGAAACAAAAGAGGATACCTTACGTa                                                                                     <  1:2606415/96‑1 (MQ=255)

             aaTCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTa                                                                                     <  1:1887180/99‑1 (MQ=255)

              aTCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTAtt                                                                                   >  1:441858/1‑100 (MQ=255)

              aTCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTAtt                                                                                   >  1:3047596/1‑100 (MQ=255)

                cATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCt                                                                                 >  1:1304229/1‑100 (MQ=255)

                cATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCt                                                                                 <  2:1933374/100‑1 (MQ=255)

                 aTGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCtt                                                                                <  2:1784797/100‑1 (MQ=255)

                 aTGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCtt                                                                                >  2:746076/1‑100 (MQ=255)

                 aTGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCtt                                                                                >  1:148152/1‑100 (MQ=255)

                 aTGGGTGGAAAACGCCGATGGTGTAAAGGAAGATAGTTTTATCGTGACGCCTCCTCTGTTTGAGAAAAAAAAAAAAAAAGAGGATACCTTACGTATTc                                                                                  <  1:3411527/98‑1 (MQ=255)

                  tGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGGTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTg                                                                               <  1:1712994/100‑1 (MQ=255)

                   gggTGGAAAATGCCGATGGTGTAAAGGATAGTCGTTTTATCGTGACGCCTCCTCTGTTTAAAAAAAAAAAAAAAAAAGAGAATACCTTACGTATTCTTg                                                                               <  1:2532912/99‑1 (MQ=255)

                      tGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGc                                                                           <  2:2576388/100‑1 (MQ=255)

                       ggAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGGTACCTTACGTATTCTTGATGca                                                                          >  2:540534/1‑100 (MQ=255)

                       ggAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGca                                                                          >  2:255333/1‑100 (MQ=255)

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

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

                         aaaaTGCCGATGCTGTAAAGGATGATCGTTTTATCGTGACGCCTCCTCTGTTTAAAAAAAAAAAAAAAAAAGAGGATACCTTACGTATTCTTGATGca                                                                          <  1:3603420/98‑1 (MQ=255)

                          aaaTGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGcaaca                                                                       >  2:2890700/1‑100 (MQ=255)

                             tGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAAt                                                                    <  1:4187852/100‑1 (MQ=255)

                               ccGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATaa                                                                  >  2:1113111/1‑100 (MQ=255)

                               ccGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATaa                                                                  >  2:2738505/1‑100 (MQ=255)

                                cGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTGCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAAc                                                                 <  2:474823/100‑1 (MQ=255)

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

                                       gtAAAGGATGGTCGTTTTATCGTGACGCTTCCTCTGTTTGTGATGCGGGGAAAAAAAGAGATTACCTTACGTATTCTTGATGCAACAAAGAACCAATTg                                                           >  1:2559967/1‑99 (MQ=255)

                                        tAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATCACCAAGTGcc                                                         >  1:387939/1‑100 (MQ=255)

                                         aaaGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCca                                                        >  2:3633863/1‑100 (MQ=255)

                                         aaaGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCca                                                        >  2:607776/1‑100 (MQ=255)

                                         aaaGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCca                                                        <  2:1658325/100‑1 (MQ=255)

                                            ggATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACTAATTGCCACAg                                                     >  2:4178299/1‑100 (MQ=255)

                                               tGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGAc                                                  >  2:3507413/1‑100 (MQ=255)

                                                 gTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCcacaagatcg                                                >  2:3994649/1‑94 (MQ=255)

                                                 gTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCg                                                <  2:1130588/100‑1 (MQ=255)

                                                    gTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGa                                             >  2:1732474/1‑100 (MQ=255)

                                                     ttttATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGaa                                            >  1:3262892/1‑100 (MQ=255)

                                                       ttATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAg                                          >  1:1729115/1‑100 (MQ=255)

                                                       tgaTCGTTACGCCCCCTCTGTTTGCGATGAAGGGAAAAATAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGa                                             <  1:794412/95‑1 (MQ=255)

                                                        tATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGt                                         <  2:3298271/100‑1 (MQ=255)

                                                        tATCGTGACGCCTCCTCTGGTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGt                                         <  2:1932503/100‑1 (MQ=255)

                                                           cGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGtttn                                      <  2:3876514/100‑2 (MQ=255)

                                                           cGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTa                                      <  2:2818035/100‑1 (MQ=255)

                                                            gTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAt                                     >  1:1119047/1‑100 (MQ=255)

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

                                                            gTGACGCCTCCTCTGTTTGCGAAAAAAAAAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAt                                     <  1:3444751/100‑1 (MQ=255)

                                                             tGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAtt                                    <  2:1610793/100‑1 (MQ=255)

                                                             tGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAtt                                    >  2:2611565/1‑100 (MQ=255)

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

                                                                cGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTg                                 <  2:3227066/100‑1 (MQ=255)

                                                                 gCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTgg                                <  1:3078008/100‑1 (MQ=255)

                                                                  cctcctCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGa                               <  2:433837/100‑1 (MQ=255)

                                                                  cctcctCTATTTGCGATGAGGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTg                                 <  1:1898026/98‑1 (MQ=255)

                                                                   ctcctcTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                              <  1:1391808/100‑1 (MQ=255)

                                                                   ctcctcTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                              >  2:1490116/1‑100 (MQ=255)

                                                                   ctcctcTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                              <  2:1182245/100‑1 (MQ=255)

                                                                   ctcctcTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                              <  1:890512/100‑1 (MQ=255)

                                                                   ctcctcTGTTTGAAAAAAAAAAAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                              <  1:4202871/100‑1 (MQ=255)

                                                                    tcctcTGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATg                             <  1:2944229/100‑1 (MQ=255)

                                                                      ctctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGaa                           <  2:1203037/100‑1 (MQ=255)

                                                                        ctGTTTGCGATGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACg                         <  2:2833858/100‑1 (MQ=255)

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

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

                                                                             tGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTaaa                    >  2:297722/1‑100 (MQ=255)

                                                                              gCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAg                   >  2:3767005/1‑100 (MQ=255)

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

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

                                                                                 ataaaaagAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAg                   <  1:3284626/90‑1 (MQ=255)

                                                                                 ataaaaagAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGa                <  1:2219707/93‑1 (MQ=255)

                                                                                 aaaaaaaaaaaaaaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGc                  <  1:2461435/90‑1 (MQ=255)

                                                                                 aaaaaaaaaaaaaaaGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTa                      <  1:577308/86‑1 (MQ=255)

                                                                                 aTGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCg                 >  2:1028670/1‑99 (MQ=255)

                                                                                  tGAAGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGa                >  1:3159410/1‑99 (MQ=255)

                                                                                    aaGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTc             >  2:216388/1‑100 (MQ=255)

                                                                                     aGGGAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTcc            >  2:1880543/1‑100 (MQ=255)

                                                                                     aGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTcc            >  2:3852892/1‑100 (MQ=255)

                                                                                     aGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTcc            >  1:3540435/1‑100 (MQ=255)

                                                                                      gggAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCg           >  1:3530036/1‑100 (MQ=255)

                                                                                      gggAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCg           <  2:305845/100‑1 (MQ=255)

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

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

                                                                                        gAAAAAAAGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTc         >  1:3730048/1‑100 (MQ=255)

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

                                                                                        gAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTATGGATGAACGTTAAAGCGATTCCGTc         >  1:1041609/1‑100 (MQ=255)

                                                                                          aaaaaaGAGGATACCTTACGTATTCTTGATGTAACAAATAACTAATAGAAAAAGGACCGGGAAAGTTTAAAATGAATTAACGTTAAGGCGATTcc            >  1:2043843/1‑95 (MQ=255)

                                                                                          aaaaaaGAGGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCaa       >  2:3060079/1‑100 (MQ=255)

                                                                                            aaaaGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATg     >  1:3928237/1‑100 (MQ=255)

                                                                                               agagGATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGAt  <  2:2259474/100‑1 (MQ=255)

                                                                                               agagAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGa   >  2:239703/1‑99 (MQ=255)

                                                                                               agagAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGAt  >  2:1194101/1‑100 (MQ=255)

ACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGAT  >  NC_000913/4544499‑4544693

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