BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,508,133	G→A	P161L (CCG→CTG)	yhfX ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,133	0	G	A	97.3%	648.4 / 0.6	185	P161L (CCG→CTG)	yhfX	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (4/1); new base A (77/103); total (81/104)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

Only 100 of 192 total aligned reads displayed.

GTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGAAACTTCCCGCGCT  >  NC_000913/3508043‑3508222

 tAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGAt                                                                                        >  1:509663/1‑93 (MQ=255)

 tAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGAt                                                                                        <  1:1502682/93‑1 (MQ=255)

  aaGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATa                                                                                       <  1:65676/93‑1 (MQ=255)

  aaGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATa                                                                                       >  1:3014312/1‑93 (MQ=255)

  aaGTCCGGCGAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATa                                                                                       <  2:1966318/93‑1 (MQ=255)

   aGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATaa                                                                                      <  2:2015367/93‑1 (MQ=255)

   aGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATaa                                                                                      <  1:1990692/93‑1 (MQ=255)

    gTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATaaa                                                                                     <  2:3665121/93‑1 (MQ=255)

    gTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATaaa                                                                                     >  1:437977/1‑93 (MQ=255)

    gTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATaaa                                                                                     <  2:2306556/93‑1 (MQ=255)

         gCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGaaaa                                                                                >  1:1046960/1‑93 (MQ=255)

         gCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGaaaa                                                                                <  1:4474530/93‑1 (MQ=255)

         gCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGaaaa                                                                                <  2:2162117/93‑1 (MQ=255)

         gCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGaaaa                                                                                <  2:4090634/93‑1 (MQ=255)

          cTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAat                                                                               >  2:448402/1‑93 (MQ=255)

           tAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAAGCGCTCTCCTGGCCCAGATAAAGAAAatc                                                                              >  2:2255056/1‑93 (MQ=255)

           tAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAatc                                                                              >  2:3204059/1‑93 (MQ=255)

           tAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAatc                                                                              <  1:906605/93‑1 (MQ=255)

             aaTGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAatcat                                                                            <  1:2449352/93‑1 (MQ=255)

             aaTGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAatcat                                                                            <  1:3327095/93‑1 (MQ=255)

             aaTGCAGCCCTGGCAGATGCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAatcat                                                                            <  2:3772640/93‑1 (MQ=255)

              aTGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAatcatc                                                                           <  1:197423/93‑1 (MQ=255)

                gCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCAtagt                                                                         <  2:3587027/93‑1 (MQ=255)

                   gCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGc                                                                      >  2:4260880/1‑93 (MQ=255)

                   gCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGc                                                                      <  2:2736934/93‑1 (MQ=255)

                   gCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGc                                                                      <  2:3513667/93‑1 (MQ=255)

                   gCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGc                                                                      <  1:3432284/93‑1 (MQ=255)

                   gCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGc                                                                      <  1:2798346/93‑1 (MQ=255)

                   gCCCTGGCAGATGCTGGATTTCGACGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGc                                                                      <  1:3886757/93‑1 (MQ=255)

                   gCCCGGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGc                                                                      <  2:1717295/93‑1 (MQ=255)

                     ccTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCta                                                                    <  1:4676727/93‑1 (MQ=255)

                     ccTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCta                                                                    <  2:398802/93‑1 (MQ=255)

                     ccTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCta                                                                    >  2:1083836/1‑93 (MQ=255)

                       tGGCAGATTCTGGATTTCCGCGACAATCTCCGGCAACCCCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCtata                                                                  <  2:3563351/93‑1 (MQ=255)

                        ggCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATaa                                                                 <  1:706402/93‑1 (MQ=255)

                         gCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATaaa                                                                <  2:3578938/93‑1 (MQ=255)

                           aGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACt                                                              <  2:4621919/93‑1 (MQ=255)

                            gATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACtt                                                             >  2:255884/1‑93 (MQ=255)

                            gATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACtt                                                             >  1:1858100/1‑93 (MQ=255)

                            gATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACtt                                                             >  1:2536620/1‑93 (MQ=255)

                            gATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACtt                                                             <  1:1294372/93‑1 (MQ=255)

                            gATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTAGCCCAGATAAAGAAAATCATCGTCGCTATAAACtt                                                             >  2:2183540/1‑93 (MQ=255)

                             aTTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATAAAGAAAATCATCGTCGCTATAAACttt                                                            >  2:1089969/1‑93 (MQ=255)

                               tCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTaa                                                          <  1:2139548/93‑1 (MQ=255)

                                cTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAag                                                         <  2:4571954/93‑1 (MQ=255)

                                 tGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAagc                                                        >  1:3297167/1‑93 (MQ=255)

                                 tGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAagc                                                        <  1:4565356/93‑1 (MQ=255)

                                 tGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAagc                                                        <  1:2676124/93‑1 (MQ=255)

                                  ggATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAagca                                                       <  2:3522985/93‑1 (MQ=255)

                                    aTTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAagcagc                                                     >  2:2591486/1‑93 (MQ=255)

                                    aTTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAagcagc                                                     <  2:1173078/93‑1 (MQ=255)

                                     tttCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAgcagca                                                    <  2:1749176/93‑1 (MQ=255)

                                     tttCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAgcagca                                                    >  1:4353139/1‑93 (MQ=255)

                                     tttCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAgcagca                                                    <  2:2803868/93‑1 (MQ=255)

                                         ggCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGa                                                <  2:2692603/93‑1 (MQ=255)

                                         ggCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGa                                                <  1:3835531/93‑1 (MQ=255)

                                          gCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGAc                                               <  1:3990861/93‑1 (MQ=255)

                                           cGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACt                                              <  1:233524/93‑1 (MQ=255)

                                            gACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTg                                             >  1:1188845/1‑93 (MQ=255)

                                            gACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTg                                             <  2:54402/93‑1 (MQ=255)

                                             aCAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTgg                                            >  2:74254/1‑93 (MQ=255)

                                               aaTCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGAt                                          >  2:2215977/1‑93 (MQ=255)

                                                aTCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGAtt                                         <  2:3414267/93‑1 (MQ=255)

                                                aTCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGAtt                                         <  1:1645476/93‑1 (MQ=255)

                                                aTCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGAtt                                         <  1:2062622/93‑1 (MQ=255)

                                                   tcCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGc                                      >  1:3241376/1‑93 (MQ=255)

                                                   tcCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGc                                      <  1:2885829/93‑1 (MQ=255)

                                                    ccGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGcc                                     >  1:4111326/1‑93 (MQ=255)

                                                     cGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGccc                                    <  1:2731298/93‑1 (MQ=255)

                                                     cGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGccc                                    >  2:4005440/1‑93 (MQ=255)

                                                      ggCAACACCTGGAGGACAAAACCGCTCTCGTGGCCAAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCGCAGACTGGATTCGCCCg                                   <  2:2890596/93‑1 (MQ=255)

                                                         aaCACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGcc                                <  2:1126738/93‑1 (MQ=255)

                                                         aaCACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGcc                                >  1:4149672/1‑93 (MQ=255)

                                                         aaCACCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGcc                                <  2:1004988/93‑1 (MQ=255)

                                                          acacCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCt                               <  1:325495/93‑1 (MQ=255)

                                                           cacCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCtt                              <  1:1556824/93‑1 (MQ=255)

                                                           cacCTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCtt                              <  2:3252078/93‑1 (MQ=255)

                                                             ccgtGACGGCTAAACCACTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCa                            <  2:962442/90‑1 (MQ=255)

                                                             ccTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCa                            <  1:203861/93‑1 (MQ=255)

                                                             ccTTGAGGCCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCa                            <  1:814169/93‑1 (MQ=255)

                                                              cTTGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCAc                           <  2:1522027/93‑1 (MQ=255)

                                                               ttGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCAcc                          <  1:1260033/93‑1 (MQ=255)

                                                               ttGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCAcc                          <  2:2190356/93‑1 (MQ=255)

                                                               ttGAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCAcc                          >  2:3819562/1‑93 (MQ=255)

                                                                 gAGGGCAAAACCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCAccgc                        >  1:4091833/1‑93 (MQ=255)

                                                                       aaaaCCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACcgccgccgc                  >  1:3761745/1‑93 (MQ=255)

                                                                        aaaCCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCa                 >  1:1102077/1‑93 (MQ=255)

                                                                          aCCGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGa               >  2:634211/1‑93 (MQ=255)

                                                                            cGCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGaaa             >  2:2662113/1‑93 (MQ=255)

                                                                             gCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGAAAc            <  2:1238387/93‑1 (MQ=255)

                                                                             gCTCTCCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGAAAc            <  2:1313928/93‑1 (MQ=255)

                                                                               tctcCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGAAACtt          >  1:3163301/1‑93 (MQ=255)

                                                                               tctcCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGAAACtt          >  1:2555417/1‑93 (MQ=255)

                                                                                ctcCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGAAACTTc         <  1:4128217/93‑1 (MQ=255)

                                                                                ctcCTGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGAAACTTc         >  2:2366438/1‑93 (MQ=255)

                                                                                  ccTGGCCCATATAAAAAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCCCCCGGCCTTCACCGCCGCCGCGGAAACTTccc       <  2:3428766/93‑1 (MQ=255)

                                                                                    tGGCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGAAACTTCCcgc     >  2:3831666/1‑93 (MQ=255)

                                                                                     ggCCCAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGAAACTTCCcgcg    <  2:3653637/93‑1 (MQ=255)

                                                                                       cccAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGAAACTTCCCGCGCt  <  1:3664829/93‑1 (MQ=255)

                                                                                       cccAGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGAAACTTCCCGCGCt  >  1:327017/1‑93 (MQ=255)

GTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATAAAGAAAATCATCGTCGCTATAAACTTTAAGCAGCACAGACTGGATTCGCCCGGCCTTCACCGCCGCCGCAGAAACTTCCCGCGCT  >  NC_000913/3508043‑3508222

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