BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,720,435	A→G	L821L (TTA→TTG)	yfiQ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,720,435	0	A	G	99.5%	608.5 / ‑4.8	196	L821L (TTA→TTG)	yfiQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (99/96); total (100/96)
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.54e-01

Only 100 of 206 total aligned reads displayed.

AGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAA  >  NC_000913/2720345‑2720524

agagATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTggg                                                                                         >  1:4140643/1‑93 (MQ=255)

 gagaTCCTCGGCGTCACGCGTTCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGc                                                                                        <  2:4122733/93‑1 (MQ=255)

 gagaTCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGc                                                                                        >  2:84485/1‑93 (MQ=255)

 gagaTCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGc                                                                                        <  1:3330576/93‑1 (MQ=255)

    aTCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTa                                                                                     <  2:505202/93‑1 (MQ=255)

     tCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAg                                                                                    <  2:2113065/93‑1 (MQ=255)

     tCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGGTGGGCTTAg                                                                                    >  2:88114/1‑93 (MQ=255)

       cTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGt                                                                                  >  2:1219951/1‑93 (MQ=255)

        tCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTc                                                                                 <  1:1719059/93‑1 (MQ=255)

        tCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTc                                                                                 <  2:3020233/93‑1 (MQ=255)

        tCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTc                                                                                 <  2:393956/93‑1 (MQ=255)

        tCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTc                                                                                 >  1:192786/1‑93 (MQ=255)

        tCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGGTGGGCTTAGGTc                                                                                 >  2:2814012/1‑93 (MQ=255)

           gCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGAc                                                                              >  1:1104991/1‑93 (MQ=255)

           gCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGAc                                                                              >  1:2570538/1‑93 (MQ=255)

           gCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGAc                                                                              <  1:1039434/93‑1 (MQ=255)

           gCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGAc                                                                              >  2:1946233/1‑93 (MQ=255)

              tCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCt                                                                           >  2:3344107/1‑93 (MQ=255)

              tCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCt                                                                           >  1:765569/1‑93 (MQ=255)

                aCGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCttt                                                                         >  2:4015438/1‑92 (MQ=255)

                aCGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTa                                                                         >  1:3484589/1‑93 (MQ=255)

                 cgcgTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTaa                                                                        <  2:443083/93‑1 (MQ=255)

                  gcgTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAAt                                                                       >  1:2721068/1‑93 (MQ=255)

                   cgTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATg                                                                      >  2:638678/1‑93 (MQ=255)

                    gTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTCGGCTTAGGTCGACGCTTAATgg                                                                     >  2:4235216/1‑93 (MQ=255)

                     tGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGa                                                                    <  1:898308/93‑1 (MQ=255)

                       cGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGaaa                                                                  >  2:2262250/1‑93 (MQ=255)

                        gATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGaaaa                                                                 <  1:4478140/93‑1 (MQ=255)

                          tttCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGt                                                               <  2:2748716/93‑1 (MQ=255)

                          tttCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGt                                                               >  2:1125948/1‑93 (MQ=255)

                          tttCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTACGCTCGGATCACAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGt                                                               <  2:1157803/93‑1 (MQ=255)

                           ttCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGtt                                                              <  2:1835312/93‑1 (MQ=255)

                           ttCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGtt                                                              <  1:65806/93‑1 (MQ=255)

                            tCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTg                                                             >  1:4417238/1‑93 (MQ=255)

                             ccGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGa                                                            <  1:3373026/93‑1 (MQ=255)

                              cGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGAt                                                           <  2:587372/93‑1 (MQ=255)

                              cGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGAt                                                           >  2:3566072/1‑93 (MQ=255)

                              cGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCACTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGAt                                                           <  1:1782218/93‑1 (MQ=255)

                              cGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCACTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGAt                                                           >  1:1330647/1‑93 (MQ=255)

                               gATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGAtt                                                          >  2:3061568/1‑93 (MQ=255)

                               gATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGAtt                                                          <  2:1019129/93‑1 (MQ=255)

                                   cTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCt                                                      >  2:172255/1‑93 (MQ=255)

                                   cTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCt                                                      >  2:1165805/1‑93 (MQ=255)

                                     gATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCtat                                                    <  1:981048/93‑1 (MQ=255)

                                     gATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCtaa                                                    >  2:627096/1‑92 (MQ=255)

                                      aTAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCtata                                                   <  2:2654048/93‑1 (MQ=255)

                                       tAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATAc                                                  <  2:1807278/93‑1 (MQ=255)

                                         aCATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATAcgc                                                >  2:3485206/1‑93 (MQ=255)

                                         aCATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATAcgc                                                >  2:3121070/1‑93 (MQ=255)

                                          cATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATAcgcg                                               >  1:1841553/1‑93 (MQ=255)

                                           aTCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGTCTTAGGTCGACGTTTAATGGAAAAGTTGTTTACCTATACGCga                                              >  1:4683553/1‑93 (MQ=255)

                                           aTCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCga                                              >  2:2544685/1‑93 (MQ=255)

                                            tCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCgag                                             >  2:2043619/1‑93 (MQ=255)

                                            tCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCgaa                                             >  1:1664636/1‑92 (MQ=255)

                                                tGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCa                                         >  1:4037221/1‑93 (MQ=255)

                                                tGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCa                                         >  1:1453859/1‑93 (MQ=255)

                                                tGCCGAATTGGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGTCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCa                                         <  1:1998890/93‑1 (MQ=255)

                                                  ccGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACg                                       <  1:3617028/93‑1 (MQ=255)

                                                  ccGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACg                                       <  1:3906420/93‑1 (MQ=255)

                                                  ccGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACg                                       <  1:423856/93‑1 (MQ=255)

                                                     aaTTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGAc                                    >  1:2152708/1‑93 (MQ=255)

                                                        ttGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCCTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTAc                                 >  2:3285166/1‑93 (MQ=255)

                                                         tGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACa                                >  2:343009/1‑93 (MQ=255)

                                                         tGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACa                                >  1:1611030/1‑93 (MQ=255)

                                                          tcTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACaa                               <  1:1390500/92‑1 (MQ=255)

                                                          gCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACaa                               <  1:3821927/93‑1 (MQ=255)

                                                          gCTGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACaa                               <  1:4548472/93‑1 (MQ=255)

                                                            tGTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACg                             >  1:2665237/1‑93 (MQ=255)

                                                             gTACTGGTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGt                            <  2:1275944/93‑1 (MQ=255)

                                                               tcTGGTTCGCTCGGATCTCAAAGGGTTGGCCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCt                          <  2:4603615/92‑1 (MQ=255)

                                                                  ggTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGaa                       <  2:1631897/93‑1 (MQ=255)

                                                                  ggTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGaa                       <  1:1072101/93‑1 (MQ=255)

                                                                  ggTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGaa                       <  2:430394/93‑1 (MQ=255)

                                                                  ggTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGaa                       <  2:1401544/93‑1 (MQ=255)

                                                                   gTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAAt                      <  1:3106212/93‑1 (MQ=255)

                                                                   gTTCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAAt                      <  1:2864299/93‑1 (MQ=255)

                                                                    ttCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATg                     >  1:3266344/1‑93 (MQ=255)

                                                                    ttCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATg                     >  1:2646812/1‑93 (MQ=255)

                                                                    ttCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATg                     >  1:3392674/1‑93 (MQ=255)

                                                                     tCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATgg                    <  1:3501815/93‑1 (MQ=255)

                                                                     tCGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATgg                    <  1:853367/93‑1 (MQ=255)

                                                                      cGCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGt                   >  1:2082255/1‑93 (MQ=255)

                                                                       gCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTa                  >  1:499637/1‑93 (MQ=255)

                                                                       gCTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTa                  >  2:349677/1‑93 (MQ=255)

                                                                        cTCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTAt                 >  2:3181232/1‑93 (MQ=255)

                                                                         tCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTAtt                >  1:1588763/1‑93 (MQ=255)

                                                                         tCGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTAtt                <  2:133737/93‑1 (MQ=255)

                                                                          cGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTa               <  2:2915865/93‑1 (MQ=255)

                                                                          cGGATCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTa               <  2:3783998/93‑1 (MQ=255)

                                                                             aTCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGa            >  2:1517329/1‑93 (MQ=255)

                                                                             aTCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGa            >  1:2786201/1‑93 (MQ=255)

                                                                             aTCTCAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGa            <  1:3054302/93‑1 (MQ=255)

                                                                              tctcAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGAt           >  2:606605/1‑93 (MQ=255)

                                                                              tctcAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGAt           <  1:726798/93‑1 (MQ=255)

                                                                                tcAAAGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGc         >  2:2502903/1‑93 (MQ=255)

                                                                                   aaGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCaa      >  2:2704154/1‑93 (MQ=255)

                                                                                   aaGGGTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCaa      >  2:2366858/1‑93 (MQ=255)

                                                                                       gTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACaa  <  1:2443039/93‑1 (MQ=255)

                                                                                       gTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACaa  <  1:2295233/93‑1 (MQ=255)

                                                                                       gTTGGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACaa  <  2:3697874/93‑1 (MQ=255)

AGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAA  >  NC_000913/2720345‑2720524

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