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	97.9%	444.6 / ‑0.7	143	Q243R (CAG→CGG)	proP	proline permease II
Reads supporting (aligned to +/- strand): ref base A (0/3); new base G (79/61); total (79/64)
Fisher's exact test for biased strand distribution p-value = 8.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 147 total aligned reads displayed.

TCGGGCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCAGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGCCACCAACG  >  NC_000913/4328717‑4328896

tCGGGCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCggg                                                                                         >  1:3837719/1‑93 (MQ=255)

 cGGGCTTTACCTGCGCCATGCGCTGGAAGAGCCTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGa                                                                                        <  1:4281670/93‑1 (MQ=255)

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

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

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

   ggCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATg                                                                                      <  2:2311963/93‑1 (MQ=255)

    gCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATgc                                                                                     >  1:923674/1‑92 (MQ=255)

      tttACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGcc                                                                                   >  2:1097125/1‑93 (MQ=255)

      tttACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGcc                                                                                   >  2:1222922/1‑93 (MQ=255)

      tttACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGCGAAGGTTTGCGGGATGGcc                                                                                   >  1:4618129/1‑93 (MQ=255)

       ttACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGccc                                                                                  <  1:762457/93‑1 (MQ=255)

       ttACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGccc                                                                                  <  1:2155012/93‑1 (MQ=255)

       ttACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGCTTGCGGGATGGccc                                                                                  >  2:277228/1‑93 (MQ=255)

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

         aCCTGCGCCACGCGCTGGAAGAGACTCGGGCCTTCCAGCAGTATGTCGATAAACTGGAACAGGGGGGCCGTGAAGGTTTGCGGGATGGCCAGa                                                                                <  2:3134115/93‑1 (MQ=255)

           cTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAAAAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGaag                                                                              >  1:1325221/1‑92 (MQ=255)

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

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

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

                ccATGGGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtctc                                                                         <  2:2782078/93‑1 (MQ=255)

                ccATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtctc                                                                         <  1:1676857/93‑1 (MQ=255)

                ccATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtctc                                                                         >  2:256245/1‑93 (MQ=255)

                ccATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtctc                                                                         >  2:2634051/1‑93 (MQ=255)

                ccATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtctc                                                                         <  2:4306821/93‑1 (MQ=255)

                ccATGCACTGGAGGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtctc                                                                         <  2:4209681/93‑1 (MQ=255)

                  aTGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGt                                                                       >  1:4293170/1‑93 (MQ=255)

                  aTGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGt                                                                       >  1:4052093/1‑93 (MQ=255)

                   tGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCAGGATGGCCCGAAAGTCTCGtt                                                                      <  1:1551936/93‑1 (MQ=255)

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

                        tGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAag                                                                 <  1:1066575/93‑1 (MQ=255)

                         ggAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAaga                                                                <  1:3629144/93‑1 (MQ=255)

                          gAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTGGTTTAAagag                                                               >  1:837686/1‑93 (MQ=255)

                          gAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAagag                                                               <  1:4654029/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                   ccGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACt                                                      <  2:1816115/93‑1 (MQ=255)

                                    cGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTa                                                     >  2:88583/1‑93 (MQ=255)

                                    cGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTa                                                     >  1:3116228/1‑93 (MQ=255)

                                    cGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAAATTGCCACTa                                                     >  2:762856/1‑93 (MQ=255)

                                     ggCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGCCTCGTTTAAAGAGATTGCCACTaa                                                    <  1:1920050/93‑1 (MQ=255)

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

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

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

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

                                         ttCAAGCAGCATGTC‑AAAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACt                                               <  1:2676995/93‑1 (MQ=255)

                                               cagcaTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGgcgc                                          <  2:2767027/93‑1 (MQ=255)

                                               cagcaTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGgcgc                                          >  2:3246363/1‑93 (MQ=255)

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

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

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

                                                   aTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGcc                                      >  1:4461161/1‑93 (MQ=255)

                                                     gTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCAGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTg                                    <  1:4267127/93‑1 (MQ=255)

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

                                                       cGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGtt                                  <  1:4127578/93‑1 (MQ=255)

                                                        gATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTg                                 >  2:1531700/1‑93 (MQ=255)

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

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

                                                           aaaCTGGAACAGGGCGACCGTGAAGGTTTGCAGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACa                              <  1:1593794/93‑1 (MQ=255)

                                                            aaCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACAt                             >  2:321512/1‑93 (MQ=255)

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

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

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

                                                             aCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCAGAAAGTCTCGTTTAAAGAGATTGCCCCTAAATACTGGCGCAGCCTGTTGACATg                            <  2:1619634/93‑1 (MQ=255)

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

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

                                                               tGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTATTGACATGTa                          >  2:396208/1‑93 (MQ=255)

                                                               tGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAACCTGTTGACATGTa                          >  1:1811777/1‑93 (MQ=255)

                                                                ggAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGACACTAAATACTGGCGCAGCCTGTTGACATGTAt                         >  1:4574777/1‑93 (MQ=255)

                                                                 gAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAATCTCGGTTAAAGAGATTGCCACTAAATACTGGCTCAGCCTGTTGACATGTAtt                        >  2:3979909/1‑93 (MQ=255)

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

                                                                  aaCAGGGCGACAGTGAAGGTTTGCGGTATGGCCCGAAAGTCTAGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTg                       <  2:633237/93‑1 (MQ=255)

                                                                   aCAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTgg                      >  1:3308005/1‑93 (MQ=255)

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

                                                                      gggCGACCGTGAAGGTTTGCGTGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTTGCGCAGCCTGTTGACATGTATTGGTCt                   >  2:3664096/1‑93 (MQ=255)

                                                                      gggCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCt                   >  1:4509975/1‑93 (MQ=255)

                                                                      gggCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCt                   >  2:3974361/1‑93 (MQ=255)

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

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

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

                                                                            ccGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAAt             >  2:3707649/1‑93 (MQ=255)

                                                                            ccGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAAt             >  1:2917578/1‑93 (MQ=255)

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

                                                                              gTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTg           >  2:543466/1‑93 (MQ=255)

                                                                              gTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTg           >  2:2089402/1‑93 (MQ=255)

                                                                              gTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTAACATGTATTGGTCTGGTAATTg           >  1:1902608/1‑93 (MQ=255)

                                                                              gTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGAAGCCTGTTGACATGTATTGGTCTGGTAATTg           >  1:3939513/1‑93 (MQ=255)

                                                                               tGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGc          >  1:1977060/1‑93 (MQ=255)

                                                                               tGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGc          >  1:2461438/1‑93 (MQ=255)

                                                                                  aGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGGCTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTgtgac       >  2:3054584/1‑89 (MQ=255)

                                                                                   agTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGccacc      <  2:2510137/92‑1 (MQ=255)

                                                                                    gTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGccacca     <  2:4412713/93‑1 (MQ=255)

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

                                                                                      ttGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGCCACCAAc   <  1:3292383/93‑1 (MQ=255)

                                                                                       tGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGCCACCAACg  <  1:4029243/93‑1 (MQ=255)

                                                                                       tGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGCCACCAACg  <  1:4735577/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` `< 28 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`