BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,328,807	A→G	Q243R (CAG→CGG)	proP →	proline permease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,328,807	0	A	G	99.7%	1016.9 / ‑5.3	324	Q243R (CAG→CGG)	proP	proline permease II
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (158/165); total (158/166)
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.42e-01

Only 100 of 336 total aligned reads displayed.

TCGGGCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCAGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGCCACCAACG  >  NC_000913/4328717‑4328896

  gggCTTTACCTGCGCCATGCGCTGGAAGAGCCTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGAt                                                                                       <  2:3267152/93‑1 (MQ=255)

  gggCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGag                                                                                       >  1:1024554/1‑92 (MQ=255)

  gggCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGag                                                                                       >  1:6786539/1‑92 (MQ=255)

  gggCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGAt                                                                                       >  1:7863355/1‑93 (MQ=255)

  gggCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGAt                                                                                       <  2:4628319/93‑1 (MQ=255)

  gggCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGAt                                                                                       <  2:1956574/93‑1 (MQ=255)

   ggCTTTCCCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATg                                                                                      <  2:72397/93‑1 (MQ=255)

   ggCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATg                                                                                      >  2:1585459/1‑93 (MQ=255)

   ggCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATg                                                                                      >  1:1858794/1‑93 (MQ=255)

    gCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATgg                                                                                     <  2:1799488/93‑1 (MQ=255)

       ttACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGccc                                                                                  >  2:4375045/1‑93 (MQ=255)

       ttACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGccc                                                                                  >  2:3417912/1‑93 (MQ=255)

        tACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCg                                                                                 >  1:2863566/1‑93 (MQ=255)

          ccTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGaa                                                                               >  2:4303011/1‑93 (MQ=255)

          ccTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGaa                                                                               <  2:5296902/93‑1 (MQ=255)

          ccTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGaa                                                                               >  1:1361568/1‑93 (MQ=255)

           cTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGaaa                                                                              <  2:4901476/93‑1 (MQ=255)

           cTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGaaa                                                                              <  2:5840735/93‑1 (MQ=255)

           cTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGaaa                                                                              >  1:5738019/1‑93 (MQ=255)

            tGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAg                                                                             <  2:2389859/93‑1 (MQ=255)

               gcCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtct                                                                          >  2:114731/1‑93 (MQ=255)

                   tGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGtt                                                                      >  1:2850405/1‑93 (MQ=255)

                     cgcTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTa                                                                    <  1:4356487/93‑1 (MQ=255)

                     cgcTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTa                                                                    >  2:3773501/1‑93 (MQ=255)

                      gcTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGTGATGGCCCGAAAGTCTCGTTTaa                                                                   >  2:5652887/1‑93 (MQ=255)

                       cTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTaaa                                                                  <  1:1977130/93‑1 (MQ=255)

                          gAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAagag                                                               >  2:903434/1‑93 (MQ=255)

                           aaGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAgaga                                                              <  2:7078227/93‑1 (MQ=255)

                            agagACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGAt                                                             <  1:2679686/93‑1 (MQ=255)

                            agagACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGAt                                                             <  2:2825929/93‑1 (MQ=255)

                             gagaCTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGAtt                                                            <  2:3181133/93‑1 (MQ=255)

                             gagaCTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGAtt                                                            >  2:3672600/1‑93 (MQ=255)

                              agaCTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTg                                                           <  1:2290649/93‑1 (MQ=255)

                               gaCTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGc                                                          >  2:1982030/1‑93 (MQ=255)

                               gaCTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGc                                                          >  1:6275584/1‑93 (MQ=255)

                                aCTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGcc                                                         >  1:6457430/1‑93 (MQ=255)

                                  tCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCAc                                                       <  1:3449654/93‑1 (MQ=255)

                                     ggCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTaa                                                    >  1:2157436/1‑93 (MQ=255)

                                     ggCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTaa                                                    >  1:370104/1‑93 (MQ=255)

                                      gCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGGGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTaaa                                                   <  1:1625405/93‑1 (MQ=255)

                                      gCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTaaa                                                   <  2:4079329/93‑1 (MQ=255)

                                      gCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTaaa                                                   <  2:6298597/93‑1 (MQ=255)

                                       cGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAAt                                                  <  1:3244367/93‑1 (MQ=255)

                                        gTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATa                                                 <  2:3974217/93‑1 (MQ=255)

                                         ttCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATAc                                                >  2:1787538/1‑93 (MQ=255)

                                          tCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACt                                               <  2:6900175/93‑1 (MQ=255)

                                          tCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACt                                               >  1:6937030/1‑93 (MQ=255)

                                          tCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACt                                               <  1:978936/93‑1 (MQ=255)

                                            cagcagCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTgg                                             >  1:5173616/1‑93 (MQ=255)

                                            cagcagCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTgg                                             >  1:4697214/1‑93 (MQ=255)

                                            cagcagCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTgg                                             <  1:3137072/93‑1 (MQ=255)

                                             agcagcATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGgc                                            >  1:5521742/1‑93 (MQ=255)

                                             agaagcaTGTCGATAATCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCGGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGgc                                            <  2:613419/90‑1 (MQ=255)

                                              gcagcaTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGgcg                                           >  1:5318988/1‑93 (MQ=255)

                                              gcagcaTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGgcg                                           >  2:2002497/1‑93 (MQ=255)

                                                agcaTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCa                                         <  1:5001917/93‑1 (MQ=255)

                                                agcaTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCa                                         <  2:7167995/93‑1 (MQ=255)

                                                 gcaTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAg                                        >  2:255732/1‑93 (MQ=255)

                                                 gcaTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAg                                        >  1:1422763/1‑93 (MQ=255)

                                                   aTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGcc                                      <  2:2654580/93‑1 (MQ=255)

                                                   aTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGcc                                      <  2:502596/93‑1 (MQ=255)

                                                      tCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGt                                   >  2:4851439/1‑93 (MQ=255)

                                                      tCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGt                                   >  1:7675869/1‑93 (MQ=255)

                                                      tCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGt                                   <  1:7511139/93‑1 (MQ=255)

                                                         aTAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGa                                >  2:2873104/1‑93 (MQ=255)

                                                         aTAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGa                                >  1:6603726/1‑93 (MQ=255)

                                                          tAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGAc                               >  1:2924280/1‑93 (MQ=255)

                                                          tAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGAc                               <  1:6390043/93‑1 (MQ=255)

                                                          tAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTCGCGCAGCCTGTTGAc                               >  1:1067535/1‑93 (MQ=255)

                                                           aaaCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACa                              <  2:1749/93‑1 (MQ=255)

                                                           aaaCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACa                              >  1:5403629/1‑93 (MQ=255)

                                                           aaaCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACa                              <  2:448790/93‑1 (MQ=255)

                                                            aaCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACAt                             <  2:6643605/93‑1 (MQ=255)

                                                            aaCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACAt                             <  1:2419213/93‑1 (MQ=255)

                                                            aaCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACAt                             <  2:2732494/93‑1 (MQ=255)

                                                             aCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACt                                               >  2:5955890/1‑74 (MQ=255)

                                                             aCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATg                            >  1:4986261/1‑93 (MQ=255)

                                                             aCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATg                            <  2:741900/93‑1 (MQ=255)

                                                              cTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGt                           <  1:1002455/93‑1 (MQ=255)

                                                              cTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGt                           <  1:3936319/93‑1 (MQ=255)

                                                               tGGAACAGGGCGACCGTGAAGGTTTGCGGGGTGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGGCATgaa                          >  1:1733464/1‑91 (MQ=255)

                                                               tGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTa                          <  1:7769992/93‑1 (MQ=255)

                                                               tGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTa                          <  2:4593518/93‑1 (MQ=255)

                                                               tGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTa                          >  1:3841026/1‑93 (MQ=255)

                                                                ggAACAGGGCGACCGTGAAGGTTTGCAGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTAt                         <  2:5446154/93‑1 (MQ=255)

                                                                  aaCAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTg                       <  1:1585459/93‑1 (MQ=255)

                                                                    cAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGt                     >  2:7813230/1‑93 (MQ=255)

                                                                     aGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTc                    >  2:1202662/1‑93 (MQ=255)

                                                                      gggCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCt                   <  1:4766927/93‑1 (MQ=255)

                                                                      gggCGACCGTGAAGGTTTGAGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCt                   >  2:6602049/1‑93 (MQ=255)

                                                                      gggCGACCGTGAAGGTTGGCGGGATGGCCCGAAAGTCCCGTTTAAAGAGATCGCCACTAAATACTGGCGCAGCCTGCTGACATGGATAGGGCt                   >  2:6994930/1‑93 (MQ=255)

                                                                       ggCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTg                  >  1:3778277/1‑93 (MQ=255)

                                                                         cGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTTAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGt                >  1:4211575/1‑93 (MQ=255)

                                                                          gACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTa               >  1:6212936/1‑93 (MQ=255)

                                                                            ccGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAAt             <  2:6765925/93‑1 (MQ=255)

                                                                             cGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAAtt            >  1:6014176/1‑93 (MQ=255)

                                                                                gAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGcc         >  1:251377/1‑93 (MQ=255)

                                                                                  aGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGccac       >  2:3581719/1‑93 (MQ=255)

                                                                                   ggTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGccacc      >  2:6201456/1‑93 (MQ=255)

                                                                                     tttGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGCCACCaa    <  1:7157110/93‑1 (MQ=255)

TCGGGCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCAGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGCCACCAACG  >  NC_000913/4328717‑4328896

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