BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,935,434	A→G	E172G (GAA→GGA)	rbsK →	ribokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,935,434	0	A	G	85.8%	644.6 / 76.5	304	E172G (GAA→GGA)	rbsK	ribokinase
Reads supporting (aligned to +/- strand): ref base A (18/25); new base G (133/127); total (152/152)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 304 total aligned reads displayed.

TGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGAT  >  NC_000913/3935344‑3935525

tGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGAc                                                                                            <  2:5998244/93‑1 (MQ=255)

 ggAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACt                                                                                           <  2:2685050/93‑1 (MQ=255)

  gAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACtt                                                                                          <  1:1851393/93‑1 (MQ=255)

    aTCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGGCGACCTTc                                                                                         >  2:1735803/1‑93 (MQ=255)

    aTCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGGCTTcc                                                                                        >  2:2459078/1‑93 (MQ=255)

    aTCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTcc                                                                                        <  1:6607045/93‑1 (MQ=255)

    aTCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTcc                                                                                        <  2:326365/93‑1 (MQ=255)

     tCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCt                                                                                       <  2:341046/93‑1 (MQ=255)

       accacTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGa                                                                                     >  2:5466734/1‑93 (MQ=255)

         cacTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGGCTTCCTGACg                                                                                   >  1:2248274/1‑93 (MQ=255)

          acTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGa                                                                                  >  1:5534446/1‑93 (MQ=255)

          acTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑CCCGACTTCCCGACgg                                                                                  >  1:7251160/1‑92 (MQ=255)

           cTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGaa                                                                                 >  2:1516226/1‑93 (MQ=255)

           cTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGAACTTCCTGACGaa                                                                                 >  2:2946538/1‑93 (MQ=255)

            tCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAAc                                                                                >  2:112495/1‑93 (MQ=255)

            tCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑CCCGGCTTCCTGACGAAc                                                                                >  1:4725840/1‑93 (MQ=255)

              gAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAActg                                                                              >  2:5810349/1‑93 (MQ=255)

               aaaGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAActgc                                                                             <  2:3350393/93‑1 (MQ=255)

               aaaGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAActgc                                                                             <  2:4093735/93‑1 (MQ=255)

               aaaGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAActgc                                                                             <  1:5375431/93‑1 (MQ=255)

               aaaGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAActgc                                                                             <  1:3920518/93‑1 (MQ=255)

               aaaGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAActgc                                                                             <  2:6406870/93‑1 (MQ=255)

               aaaGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑CCGGGCTTCCTTACCAActgc                                                                             >  1:3757785/1‑93 (MQ=255)

                aaGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAActgct                                                                            >  2:667050/1‑93 (MQ=255)

                aaGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGAACTTCCTGACGAActgct                                                                            <  2:747160/93‑1 (MQ=255)

                  gtgtGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTgg                                                                          >  2:495859/1‑93 (MQ=255)

                  gtgtGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGAACTTCCTGACGAACTGCTgg                                                                          <  2:5681420/93‑1 (MQ=255)

                      gATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCt                                                                      <  2:3814662/93‑1 (MQ=255)

                       aTGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCtg                                                                     <  2:4810299/93‑1 (MQ=255)

                          gCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGAACTTCCTGACGAACTGCTGGCGCtggtg                                                                  <  1:2625019/93‑1 (MQ=255)

                          gCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑CCGGGCTTCCTGACGAACTTCCGGCGCtgggg                                                                  >  1:293723/1‑91 (MQ=255)

                           cAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGAACTTCCTGACGAACTGCTGGCGCtggtgg                                                                 <  1:277258/93‑1 (MQ=255)

                            aGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGa                                                                <  2:2143760/93‑1 (MQ=255)

                            aGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCGTAACCCGGCTCCGGCTC‑GCGAACTTCCTGACGAACTGCTGGCGCTGGTGGa                                                                >  2:4488654/1‑93 (MQ=255)

                                gcgAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACatt                                                            >  2:2043116/1‑93 (MQ=255)

                                gcgAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACatt                                                            >  2:3484060/1‑93 (MQ=255)

                                  gAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACattat                                                          <  1:5341187/93‑1 (MQ=255)

                                  gAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGGGGGCattat                                                          >  1:6028886/1‑93 (MQ=255)

                                  gAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑CCGAACATCCTCCCGACCCGCTTGCGCTGGGGGACattat                                                          >  1:7034260/1‑93 (MQ=255)

                                   aaaaTCGCCCATCAAAATAAGACTATCGTTGCGCTTACCCCGGCTCCGGCTC‑GCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACattatt                                                         <  2:2578875/93‑1 (MQ=255)

                                   aaaaTCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACattatt                                                         <  1:5705773/93‑1 (MQ=255)

                                    aaaTCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACAttatta                                                        <  2:3018730/93‑1 (MQ=255)

                                    aaaTCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACAttatta                                                        <  1:7069963/93‑1 (MQ=255)

                                    aaaTCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑CCCGACTTCCTGACGAACTGCTGGCGCTGGGGGACAttatta                                                        >  1:1113210/1‑93 (MQ=255)

                                     aaTCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTAc                                                       <  2:6825551/93‑1 (MQ=255)

                                     aaTCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTAc                                                       <  2:182183/93‑1 (MQ=255)

                                       tCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑c                                                    >  2:5320429/1‑93 (MQ=255)

                                         gCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCa                                                  <  1:101742/93‑1 (MQ=255)

                                          cccATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑CCGGGCCTCCTTACGAACACCTGGCGCGGGGGGACATTTTTAAG‑CCaa                                                 >  1:2433934/1‑93 (MQ=255)

                                            cATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGGCTTCCTGACGAACTGCTGGCGGTGGTGGGCATTATTTTG‑CCaaac                                               >  1:1918102/1‑93 (MQ=255)

                                            cATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCaaac                                               <  1:6946430/93‑1 (MQ=255)

                                             aTCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCaaacg                                              >  2:5573158/1‑93 (MQ=255)

                                             aTCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCaaacg                                              >  2:759918/1‑93 (MQ=255)

                                              tCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCaaacga                                             >  2:1218461/1‑93 (MQ=255)

                                               cAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCaaacgaa                                            <  1:6366295/93‑1 (MQ=255)

                                                aaaaTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCaaacgaaa                                           >  2:5121076/1‑93 (MQ=255)

                                                aaaaTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCaaacgaaa                                           <  1:2860160/93‑1 (MQ=255)

                                                aaaaTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCCAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCaaacgaaa                                           >  1:4877214/1‑93 (MQ=255)

                                                 aaaTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCaaacgaaac                                          <  1:4685173/93‑1 (MQ=255)

                                                  aaTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCaaacgaaacg                                         >  2:2987457/1‑93 (MQ=255)

                                                  aaTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCaaacgaaacg                                         <  2:4256246/93‑1 (MQ=255)

                                                     aaGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGaa                                      <  2:1386949/93‑1 (MQ=255)

                                                     aaGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGaa                                      <  2:3666057/93‑1 (MQ=255)

                                                     aaGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑CCGAACTTCCTGACGAACTGCTGGGGCTGGGGGGCATTATTACG‑CCAAAACAAAAGGaa                                      >  1:5824675/1‑93 (MQ=255)

                                                      aGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAg                                     >  2:6213588/1‑93 (MQ=255)

                                                       gACTGTCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCCGGTGGACATTATTACG‑CCAAACGAAACGGAAGc                                    >  2:7371624/1‑93 (MQ=255)

                                                       gACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑CCGGGCTTCCTGACCAACTGCTGGGGCTGGGGGGCATTATTAAG‑CCCAACGAAACGGGAGc                                    >  1:2111391/1‑93 (MQ=255)

                                                       gACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑CCCGGCTTCCTGACGAACTGCTGGCGCCGGGGGGCATTTTTACGCCCAAACAAAAGGGagc                                     >  1:6737781/1‑91 (MQ=255)

                                                        aCTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCa                                   <  1:7527647/93‑1 (MQ=255)

                                                        aCTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCa                                   >  2:5631004/1‑93 (MQ=255)

                                                        aCTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCa                                   <  2:6863144/93‑1 (MQ=255)

                                                          tATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGa                                 >  2:2824614/1‑93 (MQ=255)

                                                            tCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGaaa                               >  2:5942709/1‑93 (MQ=255)

                                                            tCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGaaa                               <  1:5929670/93‑1 (MQ=255)

                                                             cGTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGaaaa                              <  1:5307006/93‑1 (MQ=255)

                                                              gTTGCGCTTAACCCGGCTCCGGCTC‑GCGGGCTTCCTGACGAACTGCTGGCGCTGGGGGGCATTATTACG‑CCAAACGAAACGGAAGCAGaaaaa                             >  1:3101744/1‑92 (MQ=255)

                                                              gTTGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGGGGGCCTTATTACG‑CCAAACAAAACGGAAGCAGaaaaa                             >  1:330549/1‑92 (MQ=255)

                                                               ttGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGc                            >  2:1294548/1‑93 (MQ=255)

                                                                tGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCt                           <  1:5043098/93‑1 (MQ=255)

                                                                tGCGCTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCt                           <  2:7120003/93‑1 (MQ=255)

                                                                tGCGCTTAACCCGGCTCCGGCTC‑CCGGGCTTCCCTACGAACTGCTGGCGCTGGGGGGCATTATTTCC‑CCAAACGAAACGGAAGCAGAAAAACt                           >  1:4476083/1‑93 (MQ=255)

                                                                 gcgcTTAACCCGGCTCCGGCTC‑CCCAAATTCCTGACGAACTGCTGGCGCCGGGGGGCATTATTTCG‑CCAAACGAAACGGGAACAGAAAAGCTc                          >  1:715314/1‑93 (MQ=255)

                                                                  cgcTTAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCa                         <  1:4216701/93‑1 (MQ=255)

                                                                    cTTAACCCGGCTCCGGCTC‑GCGGGCTTCCTGACGAACTGCTGGCGCTGGGGGGCATTATTACG‑CCAAACGAAAAGGGAGCAGAAAAAATCAcc                       >  1:1644825/1‑93 (MQ=255)

                                                                     ttAACCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCg                      >  2:2195170/1‑93 (MQ=255)

                                                                     ttAACCCGGCTCCGGCTC‑GCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCg                      <  1:126141/93‑1 (MQ=255)

                                                                        aCCCGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTa                   >  2:6958776/1‑93 (MQ=255)

                                                                         ccccGCTCCGGCTC‑GGGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTAt                  <  2:7158672/93‑1 (MQ=255)

                                                                         cccGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTAt                  <  1:2806396/93‑1 (MQ=255)

                                                                         cccGGCTCCGGCTC‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTAt                  >  2:3592567/1‑93 (MQ=255)

                                                                           cggctccggctc‑GCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTc                >  2:2301462/1‑93 (MQ=255)

                                                                             gctccggctc‑GCGGACTTCCTGACGAACTGCTGGGGCTGGGGGGCATTATTACG‑CCAAACGAAACGGAAGCAGAAAAAATCACCGGGATTCgt              >  1:878897/1‑93 (MQ=255)

                                                                               tccggctc‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCgtgt            <  1:6421972/93‑1 (MQ=255)

                                                                               tccggctc‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCgtgt            <  1:300260/93‑1 (MQ=255)

                                                                                ccggctc‑GCGGGCTTCCTGACGAAACGCTGGCGCTGGGGGGCATTATTACG‑CCAAACGAAACGGAAGCAAAAAAACTCACCCGGATTCGTGtt           >  1:5186125/1‑93 (MQ=255)

                                                                                ccggctc‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGtt           <  1:988242/93‑1 (MQ=255)

                                                                                  ggctc‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGa         <  1:3553425/93‑1 (MQ=255)

                                                                                    ctc‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGaaa       >  2:6565426/1‑93 (MQ=255)

                                                                                     tc‑GCGGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGaaaa      >  1:2469410/1‑93 (MQ=255)

                                                                                        gcgGACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAatg    >  1:6151702/1‑93 (MQ=255)

TGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC‑GCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACG‑CCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGAT  >  NC_000913/3935344‑3935525

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