BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	857,539	T→C	D249G (GAC→GGC)	ybiU ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	857,539	0	T	C	97.6%	385.2 / ‑0.3	126	D249G (GAC→GGC)	ybiU	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (52/70); total (53/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.30e-01

Only 100 of 141 total aligned reads displayed.

ACGACGTGCAGCAGCCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAAC  >  NC_000913/857449‑857627

acgacgTGCAGCAGCCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGccc                                                                                        >  2:2714175/1‑93 (MQ=255)

  gacgTGCAGCAGCCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCAc                                                                                      <  2:4107791/93‑1 (MQ=255)

  gacgTGCAGCAGCCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCAc                                                                                      >  1:282258/1‑93 (MQ=255)

   acgTGCAGCAGCCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCAcc                                                                                     >  2:2612725/1‑93 (MQ=255)

        cagcagCCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTa                                                                                <  1:1421516/93‑1 (MQ=255)

        cagcagCCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTa                                                                                <  1:2266708/93‑1 (MQ=255)

         agcagcCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTAc                                                                               <  2:389588/93‑1 (MQ=255)

         agcagcCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTAc                                                                               >  2:2905698/1‑93 (MQ=255)

          gcagcCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTAct                                                                              >  1:1000327/1‑93 (MQ=255)

          gcagcCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTAct                                                                              >  2:4672388/1‑93 (MQ=255)

          gcagcCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTAct                                                                              >  2:4339348/1‑93 (MQ=255)

              ccccTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTc                                                                          <  2:1481130/93‑1 (MQ=255)

               cccTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCa                                                                         >  2:4126098/1‑93 (MQ=255)

               cccTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCa                                                                         >  1:4340506/1‑93 (MQ=255)

                    aCCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTc                                                                    <  1:1977424/93‑1 (MQ=255)

                    aCCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTc                                                                    <  1:30700/93‑1 (MQ=255)

                     ccAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTcc                                                                   >  1:1797591/1‑93 (MQ=255)

                      cAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACCTTTGGCGGTGTTGCCCACCGTGTCCTCTTCACCTTCCg                                                                  >  1:522119/1‑93 (MQ=255)

                        ggCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTa                                                                <  2:1041119/93‑1 (MQ=255)

                        ggCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTa                                                                <  1:1090618/93‑1 (MQ=255)

                         gcagcaTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTAc                                                               <  1:3542685/93‑1 (MQ=255)

                         gcagcaTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTAc                                                               <  1:3065658/93‑1 (MQ=255)

                           agcaTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGa                                                             >  2:1468392/1‑93 (MQ=255)

                           agcaTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGa                                                             <  2:2793860/93‑1 (MQ=255)

                           agcaTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGa                                                             <  1:3307652/93‑1 (MQ=255)

                            gcaTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAt                                                            <  1:1432923/93‑1 (MQ=255)

                             caTATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtg                                                           <  2:874331/93‑1 (MQ=255)

                              atatCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtgt                                                          <  2:2744129/93‑1 (MQ=255)

                              atatCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtgt                                                          >  2:3588533/1‑93 (MQ=255)

                              atatCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtgt                                                          <  2:3987043/93‑1 (MQ=255)

                               tatCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGAtgtg                                                         >  1:2727429/1‑93 (MQ=255)

                                 tCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGcc                                                       <  1:2487500/93‑1 (MQ=255)

                                 tCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGcc                                                       >  1:2755858/1‑93 (MQ=255)

                                   agagagCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGc                                                     >  1:2472544/1‑93 (MQ=255)

                                     agagCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCAt                                                   <  1:3412746/93‑1 (MQ=255)

                                     agagCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCAt                                                   <  2:832435/93‑1 (MQ=255)

                                      gagCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATg                                                  <  1:1495212/93‑1 (MQ=255)

                                        gcgcTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCTTGTACTCTTCAACTTCCGTACGATGTGCCGCATGcc                                                <  1:1151193/93‑1 (MQ=255)

                                        gcgcTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGCACTCTTCAACTTCCGTACGATGTGCCGCATGcc                                                >  2:4570855/1‑93 (MQ=255)

                                         cgcTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCa                                               <  1:898421/93‑1 (MQ=255)

                                          gcTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAg                                              <  1:931379/93‑1 (MQ=255)

                                            tGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAggg                                            <  2:2669837/93‑1 (MQ=255)

                                            tGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAggg                                            <  1:3888956/93‑1 (MQ=255)

                                             gTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGa                                           <  1:1101594/93‑1 (MQ=255)

                                               ccAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATc                                         >  1:371921/1‑93 (MQ=255)

                                                  gCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCata                                      >  2:337/1‑93 (MQ=255)

                                                  gCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCata                                      <  2:2262396/93‑1 (MQ=255)

                                                  gCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCata                                      <  2:2262395/93‑1 (MQ=255)

                                                    ccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCTTACGATGTGCCGCATGCCAGGGATCATAtt                                    >  1:1406804/1‑93 (MQ=255)

                                                     cTGGAATGTCCGAAACACGGAACATTTTGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTg                                   >  2:4668842/1‑93 (MQ=255)

                                                     cTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTg                                   >  1:2504179/1‑93 (MQ=255)

                                                     cTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTg                                   >  2:4630375/1‑93 (MQ=255)

                                                       ggAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTgcg                                 >  1:1756772/1‑93 (MQ=255)

                                                        gAATGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTgcgc                                >  2:1120422/1‑93 (MQ=255)

                                                          aTGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCa                              >  2:2141576/1‑93 (MQ=255)

                                                          aTGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCa                              <  2:177292/93‑1 (MQ=255)

                                                          aTGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCa                              >  1:3657124/1‑93 (MQ=255)

                                                          aTGTCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCa                              <  2:2852209/93‑1 (MQ=255)

                                                           tGTCCGAAACAGTGAACATGTTGTGGGTGTGCCCACCGTGTTCTCTTCAACTTCCGTACGATGTGCCGCATGCCAAGGATCATATTGCACCAg                             <  1:4176560/93‑1 (MQ=255)

                                                             tCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGAt                           <  2:2069103/93‑1 (MQ=255)

                                                             tCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGAt                           >  1:406390/1‑93 (MQ=255)

                                                             tCCGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGAt                           >  1:406389/1‑93 (MQ=255)

                                                              ccGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGAtt                          >  2:189837/1‑93 (MQ=255)

                                                              ccGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGAtt                          >  2:537358/1‑93 (MQ=255)

                                                              ccGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGAtt                          <  2:631947/93‑1 (MQ=255)

                                                                gAAACACGGAACATTTGGAGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGc                        >  2:1652680/1‑93 (MQ=255)

                                                                 aaaCACGGAACATTTGGTGGTGTTGCCCACCGTGTAATCTTCAACTTCCGTACGATGTGCCGAATGCCAGGGATCATATTGCGCCAGATTGcc                       >  1:1208353/1‑93 (MQ=255)

                                                                  aaCACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCa                      <  1:2164002/93‑1 (MQ=255)

                                                                  aaCACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCa                      >  1:4670482/1‑93 (MQ=255)

                                                                  aaCACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTACCa                      <  2:4644428/93‑1 (MQ=255)

                                                                   acacGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCAt                     <  1:1990957/93‑1 (MQ=255)

                                                                   acacGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCAt                     >  1:774364/1‑93 (MQ=255)

                                                                    cacGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCAtt                    <  1:1709617/93‑1 (MQ=255)

                                                                      cGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTaa                  >  2:2060420/1‑93 (MQ=255)

                                                                         aaCATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGa               >  2:2565392/1‑93 (MQ=255)

                                                                           cATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACg             <  2:4695695/93‑1 (MQ=255)

                                                                           cATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACg             >  2:3377746/1‑93 (MQ=255)

                                                                            aTTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGt            <  1:3175392/93‑1 (MQ=255)

                                                                            aTTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGt            >  2:2901792/1‑93 (MQ=255)

                                                                             tttGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGtt           <  2:1543665/93‑1 (MQ=255)

                                                                             tttGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGtt           >  1:2092264/1‑93 (MQ=255)

                                                                             tttGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGtt           <  1:2458358/93‑1 (MQ=255)

                                                                             tttGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGtt           <  1:1466790/93‑1 (MQ=255)

                                                                             tttGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGtt           <  2:1207562/93‑1 (MQ=255)

                                                                             tttGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGtt           <  1:2359745/93‑1 (MQ=255)

                                                                              ttGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTg          <  2:4008267/93‑1 (MQ=255)

                                                                              ttGGTGGTGTTGCCCACCGTGTACTCTTCAACTTACGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTg          <  2:3566600/93‑1 (MQ=255)

                                                                               tggtggTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTgg         <  1:3438584/93‑1 (MQ=255)

                                                                               tggtggTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATAATATTGCGCCAGATTGCCATTAAAGACGTTgg         <  2:1867126/93‑1 (MQ=255)

                                                                                ggtggtGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGc        <  1:3542905/93‑1 (MQ=255)

                                                                                ggtggtGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGc        >  2:2787555/1‑93 (MQ=255)

                                                                                 gtggtgTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCg       >  1:2010714/1‑93 (MQ=255)

                                                                                  tggtgTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGa      <  1:4537245/93‑1 (MQ=255)

                                                                                     tgttgcCCGCCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGaaaa   <  1:4024985/93‑1 (MQ=255)

                                                                                     tgttgcCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGNCATTAAAGACGTTGGCGaaaa   >  1:3567614/1‑93 (MQ=255)

                                                                                     tgttgcCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGaaaa   >  2:44800/1‑93 (MQ=255)

                                                                                     tgttgcCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGaaaa   <  1:3705566/93‑1 (MQ=255)

                                                                                     tgttgcCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGaaaa   <  1:3551141/93‑1 (MQ=255)

                                                                                      gttgcCCACCGTGTACTCTTCCACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGCTTGCCATTAAAGTCGTTGGCGAacac  >  1:671475/1‑93 (MQ=255)

                                                                                      gttgcCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAAc  >  1:671476/1‑93 (MQ=255)

ACGACGTGCAGCAGCCCCTGACCAGGCAGCATATCAGAGAGCGCTGTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAAC  >  NC_000913/857449‑857627

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