BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,125,900	A→G	Q174R (CAG→CGG)	yceH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,125,900	0	A	G	82.8%	471.5 / 93.0	204	Q174R (CAG→CGG)	yceH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (18/17); new base G (83/86); total (101/103)
Fisher's exact test for biased strand distribution p-value = 8.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

Only 100 of 213 total aligned reads displayed.

CGAAGATGGTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAG  >  NC_000913/1125808‑1125991

  aaGATGGTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGc                                                                                           >  2:1503072/1‑93 (MQ=255)

   aGATGGTCCTTTTGTGGTGTGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGcc                                                                                          >  1:6490309/1‑93 (MQ=255)

   aGATGGTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGcc                                                                                          >  2:4013267/1‑93 (MQ=255)

   aGATGGTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGcc                                                                                          <  2:1438256/93‑1 (MQ=255)

   aGATGGTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGcc                                                                                          <  2:4033581/93‑1 (MQ=255)

      tGGTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATcggccggc                                                                                       >  2:1779254/1‑93 (MQ=255)

      tGGTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATcggccggc                                                                                       >  2:4353737/1‑93 (MQ=255)

      tGGTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATcggccggc                                                                                       >  1:4264855/1‑93 (MQ=255)

      tGGTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATcggccggc                                                                                       <  2:516707/93‑1 (MQ=255)

      tGGTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCcggc                                                                                       <  2:4954318/93‑1 (MQ=255)

       ggTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCcggcg                                                                                      >  2:5528319/1‑93 (MQ=255)

       ggTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCcggcg                                                                                      <  1:571424/93‑1 (MQ=255)

        gTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCcggcgg                                                                                     <  1:5806092/93‑1 (MQ=255)

         tCCTTTTGTGGTGTGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGt                                                                                    <  1:4697716/93‑1 (MQ=255)

         tCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGt                                                                                    <  2:2600044/93‑1 (MQ=255)

           cTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGa                                                                                  <  1:6303596/93‑1 (MQ=255)

           cTTTTGTGGTGCGTCTGGCCCGCGAACCGGGGAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGa                                                                                  >  1:5937094/1‑93 (MQ=255)

            ttttGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGAc                                                                                 <  2:3500160/93‑1 (MQ=255)

              ttGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACgg                                                                               <  1:3684084/93‑1 (MQ=255)

                 tggtgCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGata                                                                            <  2:3147415/93‑1 (MQ=255)

                  ggtgCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGatat                                                                           <  1:3039550/93‑1 (MQ=255)

                    tgCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGt                                                                         <  2:2072101/93‑1 (MQ=255)

                       gTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGa                                                                      <  2:3818013/93‑1 (MQ=255)

                       gTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGa                                                                      >  1:1055076/1‑93 (MQ=255)

                        tCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGaa                                                                     <  2:1090592/93‑1 (MQ=255)

                         cTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAAc                                                                    <  2:6013905/93‑1 (MQ=255)

                          tGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAAcg                                                                   >  2:6183023/1‑93 (MQ=255)

                          tGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAAcg                                                                   >  1:3373/1‑93 (MQ=255)

                           ggCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAAcgc                                                                  >  2:199550/1‑93 (MQ=255)

                           ggCCCACGAACCGAGTAAACGCGACATCCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAAcgc                                                                  <  1:6516009/93‑1 (MQ=255)

                            gCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAAcgcg                                                                 >  2:2689806/1‑93 (MQ=255)

                             cccGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGTCGGATATGGCGACCGCgg                                                                >  1:1903588/1‑93 (MQ=255)

                              ccGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGt                                                               <  1:455548/93‑1 (MQ=255)

                              ccGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGt                                                               <  2:6599377/93‑1 (MQ=255)

                              ccGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGt                                                               <  1:1654868/93‑1 (MQ=255)

                               cgcgAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGtt                                                              <  2:5563991/93‑1 (MQ=255)

                               cgcgAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGtt                                                              <  1:2239355/93‑1 (MQ=255)

                                gcgAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTg                                                             >  1:5857840/1‑93 (MQ=255)

                                  gAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGGTGAc                                                           >  1:3884602/1‑93 (MQ=255)

                                   aaCCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAAAGCGGTTGACg                                                          >  1:2466757/1‑93 (MQ=255)

                                     ccGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGt                                                        <  1:5772226/93‑1 (MQ=255)

                                      cGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTg                                                       >  1:6146821/1‑93 (MQ=255)

                                       gggTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGa                                                      <  1:5495934/93‑1 (MQ=255)

                                          tAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGAttt                                                   >  1:4477047/1‑93 (MQ=255)

                                             aCGGGAAAACCGCTACATGCATCTTTTCAGTGGTGATGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACa                                                <  1:1383026/93‑1 (MQ=255)

                                              cgcgAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGCTTTTCAg                                               >  2:16715/1‑93 (MQ=255)

                                              cgcgAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAg                                               <  1:2758653/93‑1 (MQ=255)

                                              agcgAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAg                                               <  2:1736599/92‑1 (MQ=255)

                                               gcgAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAgg                                              >  1:4966631/1‑93 (MQ=255)

                                                cgAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGc                                             <  2:4624821/93‑1 (MQ=255)

                                                cgAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGGTGACGGTGATTTACAGGc                                             >  2:5508883/1‑93 (MQ=255)

                                                 gAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGcc                                            <  1:1517531/93‑1 (MQ=255)

                                                 gAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGcc                                            <  2:1418489/93‑1 (MQ=255)

                                                 gAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGcc                                            >  2:3514062/1‑93 (MQ=255)

                                                  aaaaCCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGccc                                           <  2:221180/93‑1 (MQ=255)

                                                  aaaaCCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGccc                                           <  1:3657105/93‑1 (MQ=255)

                                                   aagccGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCcg                                          <  1:68822/90‑1 (MQ=255)

                                                    acccGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCcgc                                         <  2:2473992/91‑1 (MQ=255)

                                                    aaCCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCcgc                                         <  1:41134/93‑1 (MQ=255)

                                                    aaCCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCcgc                                         <  1:6032428/93‑1 (MQ=255)

                                                       cGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTc                                      >  2:4499556/1‑93 (MQ=255)

                                                       cGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTc                                      <  2:331085/93‑1 (MQ=255)

                                                         cTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGa                                    <  2:4755111/93‑1 (MQ=255)

                                                         cTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGa                                    >  2:4527995/1‑93 (MQ=255)

                                                         cTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGa                                    >  2:5858437/1‑93 (MQ=255)

                                                         cTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGa                                    <  2:3086974/93‑1 (MQ=255)

                                                         cTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGa                                    <  2:1014067/93‑1 (MQ=255)

                                                          tACATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGaa                                   <  2:1425395/93‑1 (MQ=255)

                                                          tACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGaa                                   >  1:4762793/1‑93 (MQ=255)

                                                           aCATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAg                                  >  2:2057348/1‑93 (MQ=255)

                                                           aCATGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAg                                  >  2:5339954/1‑93 (MQ=255)

                                                             aTGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGcc                                >  2:735067/1‑93 (MQ=255)

                                                             aTGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGcc                                >  2:1586881/1‑93 (MQ=255)

                                                              tGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGccc                               <  1:5594584/93‑1 (MQ=255)

                                                              tGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGccc                               <  1:5385556/93‑1 (MQ=255)

                                                              tGCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGGGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGccc                               >  1:2660507/1‑93 (MQ=255)

                                                               gCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCt                              <  2:2850221/93‑1 (MQ=255)

                                                               gCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCt                              <  1:2837570/93‑1 (MQ=255)

                                                               gCATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCt                              <  1:4795752/93‑1 (MQ=255)

                                                               gCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCt                              <  2:5224587/93‑1 (MQ=255)

                                                                cATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTg                             >  1:3537812/1‑93 (MQ=255)

                                                                cATCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTg                             <  2:6555328/93‑1 (MQ=255)

                                                                 aTCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTgg                            <  2:3009272/93‑1 (MQ=255)

                                                                  tCTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGa                           <  2:3067826/93‑1 (MQ=255)

                                                                   cTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGaa                          >  1:4756115/1‑93 (MQ=255)

                                                                   cTTTTCAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGaa                          <  1:2337979/93‑1 (MQ=255)

                                                                    ttttCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGaaa                         >  2:391196/1‑93 (MQ=255)

                                                                        cAGTGGTGAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGa                     >  1:1570181/1‑93 (MQ=255)

                                                                          gtggtgAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAg                   >  2:6370216/1‑93 (MQ=255)

                                                                           tggtgAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGt                  <  1:448310/93‑1 (MQ=255)

                                                                           tggtgAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGt                  <  2:1937083/93‑1 (MQ=255)

                                                                             gtgAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTgg                <  1:1636692/93‑1 (MQ=255)

                                                                              tgAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGc               >  2:3689783/1‑93 (MQ=255)

                                                                               gAGGTTGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCa              <  2:5638231/93‑1 (MQ=255)

                                                                               gAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCa              >  1:1761897/1‑93 (MQ=255)

                                                                                   ttGAAGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAAc          <  2:3580624/93‑1 (MQ=255)

                                                                                      aaGATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGa       <  2:5679120/93‑1 (MQ=255)

                                                                                       aGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGaa      >  1:6356449/1‑93 (MQ=255)

                                                                                        gATCGGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGaaa     <  1:2714154/93‑1 (MQ=255)

                                                                                           cAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAg  >  1:2352614/1‑93 (MQ=255)

CGAAGATGGTCCTTTTGTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAG  >  NC_000913/1125808‑1125991

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