BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	52,554	A→G	L288P (CTG→CCG)	pdxA ←	PdxA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	52,554	0	A	G	99.5%	644.2 / ‑4.8	203	L288P (CTG→CCG)	pdxA	PdxA
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (103/99); total (103/100)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 220 total aligned reads displayed.

AGCGCCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCAGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCACGGCGTCG  >  NC_000913/52464‑52643

aGCGCCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCggg                                                                                         >  1:111618/1‑93 (MQ=255)

aGCGCCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCggg                                                                                         <  1:2087337/93‑1 (MQ=255)

aGCGCCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCggg                                                                                         <  2:6221789/93‑1 (MQ=255)

 gcgcCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAgggcgggc                                                                                        >  1:3923444/1‑93 (MQ=255)

  cgcCGTAATAAAACTGCCGACATCGGCTTTGCCATGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGcc                                                                                       >  1:5692139/1‑93 (MQ=255)

  cgcCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGcc                                                                                       <  2:2555569/93‑1 (MQ=255)

       taataaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCg                                                                                  <  2:358391/93‑1 (MQ=255)

          taaACCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTaa                                                                               <  1:194118/93‑1 (MQ=255)

          taaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTaa                                                                               <  1:6505513/93‑1 (MQ=255)

          taaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTaa                                                                               >  2:2625707/1‑93 (MQ=255)

          taaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGCGCGGGCCCAGCGTaa                                                                               >  2:3137822/1‑93 (MQ=255)

           aaaaCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAAt                                                                              >  2:1797015/1‑93 (MQ=255)

            aaaCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATg                                                                             >  2:3480368/1‑93 (MQ=255)

            aaaCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATg                                                                             <  1:4723014/93‑1 (MQ=255)

            aaaCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATg                                                                             <  1:3756615/93‑1 (MQ=255)

            aaaCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGGCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATg                                                                             >  2:4021043/1‑93 (MQ=255)

              aCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGtt                                                                           <  2:1362678/93‑1 (MQ=255)

               cTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTc                                                                          >  1:2645107/1‑93 (MQ=255)

               cTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTc                                                                          >  2:4593125/1‑93 (MQ=255)

               cTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTc                                                                          <  2:2491704/93‑1 (MQ=255)

                tGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTca                                                                         <  1:2923600/93‑1 (MQ=255)

                 gCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTcac                                                                        >  2:4428367/1‑93 (MQ=255)

                     aCATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACcg                                                                    >  2:4045918/1‑93 (MQ=255)

                       aTCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACcgcg                                                                  >  1:4412651/1‑93 (MQ=255)

                       aTCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACcgcg                                                                  >  1:6392719/1‑93 (MQ=255)

                        tCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCgcgc                                                                 >  2:1859679/1‑93 (MQ=255)

                        tCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCgcgc                                                                 >  2:396816/1‑93 (MQ=255)

                        tCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCgcgc                                                                 >  1:3049156/1‑93 (MQ=255)

                           gCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCATCGTAATGTTCACACCGCGCCCg                                                              >  1:1658705/1‑93 (MQ=255)

                           gCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCg                                                              <  2:5210992/93‑1 (MQ=255)

                              ttGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAg                                                           >  2:6663160/1‑93 (MQ=255)

                               tGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGc                                                          <  2:2835906/93‑1 (MQ=255)

                                gCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGcc                                                         >  2:6058088/1‑93 (MQ=255)

                                 ccACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGccc                                                        <  2:2745010/93‑1 (MQ=255)

                                 ccACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGccc                                                        <  2:5001815/93‑1 (MQ=255)

                                  cACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCt                                                       >  2:6524781/1‑93 (MQ=255)

                                  cACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCt                                                       <  2:2269771/93‑1 (MQ=255)

                                   aCGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTg                                                      >  2:2865189/1‑93 (MQ=255)

                                   aCGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTg                                                      >  1:287192/1‑93 (MQ=255)

                                    cGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTgg                                                     <  1:598532/93‑1 (MQ=255)

                                      gcccGCCAGTTCAAGCGCGGTGCCGTCGTCCACTGATGTGCGAATAAAGGCCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTa                                                   <  1:4118540/92‑1 (MQ=255)

                                        ccGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTAtt                                                 <  2:5457932/93‑1 (MQ=255)

                                        ccGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTAtt                                                 >  2:420013/1‑93 (MQ=255)

                                        ccGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTAtt                                                 <  2:4635342/93‑1 (MQ=255)

                                          gCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTAtttt                                               <  1:6137914/93‑1 (MQ=255)

                                           ccAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTa                                              <  2:5084318/93‑1 (MQ=255)

                                             aGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGc                                            <  2:1156526/93‑1 (MQ=255)

                                              gTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCa                                           <  1:3433482/93‑1 (MQ=255)

                                              gTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCa                                           >  2:5295535/1‑93 (MQ=255)

                                                 cAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGTGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACgg                                        <  2:3876242/93‑1 (MQ=255)

                                                 cAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACgg                                        >  2:6456348/1‑93 (MQ=255)

                                                   aGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGa                                      >  1:3366959/1‑93 (MQ=255)

                                                     cgcgGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAg                                    <  2:6314144/93‑1 (MQ=255)

                                                      gcgGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGa                                   <  2:4467461/93‑1 (MQ=255)

                                                        ggTGCCGTGGTCCACTGATGTGCGAATAAAGGGCCGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGAcc                                 <  1:5244863/93‑1 (MQ=255)

                                                         gTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGAccc                                >  2:5774490/1‑93 (MQ=255)

                                                          tGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCt                               <  1:5132120/93‑1 (MQ=255)

                                                          tGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCt                               <  1:3255935/93‑1 (MQ=255)

                                                          tGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCt                               <  2:3462200/93‑1 (MQ=255)

                                                          tGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCt                               <  1:2174845/93‑1 (MQ=255)

                                                            ccGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGa                             <  1:2259113/93‑1 (MQ=255)

                                                             cGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGAt                            <  1:3966889/93‑1 (MQ=255)

                                                              gTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATc                           >  2:4809891/1‑93 (MQ=255)

                                                              gTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATc                           >  2:2980853/1‑93 (MQ=255)

                                                              gTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATc                           >  1:4608604/1‑93 (MQ=255)

                                                              gTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGGATTTTAGCACGGGAAGACCCTGATc                           >  2:6481572/1‑93 (MQ=255)

                                                                ggTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGt                         <  1:6664892/93‑1 (MQ=255)

                                                                 gTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTg                        >  1:5722826/1‑93 (MQ=255)

                                                                 gTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTg                        >  1:1393392/1‑93 (MQ=255)

                                                                   ccACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGt                      <  2:6236487/93‑1 (MQ=255)

                                                                   ccACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGt                      <  1:5524469/93‑1 (MQ=255)

                                                                    cACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTa                     >  1:925726/1‑93 (MQ=255)

                                                                    cACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTa                     >  1:1790869/1‑93 (MQ=255)

                                                                      cTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACa                   <  2:5939397/93‑1 (MQ=255)

                                                                       tGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACAt                  <  2:6578388/93‑1 (MQ=255)

                                                                       tGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACAt                  <  2:6035645/93‑1 (MQ=255)

                                                                        gATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATc                 <  1:2966610/93‑1 (MQ=255)

                                                                        gATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATc                 <  1:4597590/93‑1 (MQ=255)

                                                                        gATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATc                 <  2:3393485/93‑1 (MQ=255)

                                                                        gATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATc                 >  1:3705637/1‑93 (MQ=255)

                                                                         aTGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCg                >  2:5515761/1‑93 (MQ=255)

                                                                           gtgCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGcc              <  1:4942266/93‑1 (MQ=255)

                                                                           gtgCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGcc              >  1:3324390/1‑93 (MQ=255)

                                                                             gCGAATAAAGGGCGGGTCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCCCGGGAAGACCCTGATCGTGGTACATCGCCAg            <  1:692406/93‑1 (MQ=255)

                                                                             gCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAg            <  1:5943441/93‑1 (MQ=255)

                                                                             gCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAg            >  1:4459801/1‑93 (MQ=255)

                                                                               gAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCa          >  2:1515391/1‑93 (MQ=255)

                                                                               gAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCa          >  1:5358274/1‑93 (MQ=255)

                                                                               gAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCa          <  1:331298/93‑1 (MQ=255)

                                                                               gAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCa          <  1:239658/93‑1 (MQ=255)

                                                                                 aTAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcg        >  2:1949844/1‑93 (MQ=255)

                                                                                 aTAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcg        >  1:5628914/1‑93 (MQ=255)

                                                                                  tAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcgg       <  2:6355167/93‑1 (MQ=255)

                                                                                   aaaGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggc      <  1:1311089/93‑1 (MQ=255)

                                                                                    aaGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcg     >  2:3440213/1‑93 (MQ=255)

                                                                                    aaGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcg     >  2:1905844/1‑93 (MQ=255)

                                                                                     aGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgt    >  2:5072213/1‑93 (MQ=255)

                                                                                     aGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgt    <  2:2417013/93‑1 (MQ=255)

                                                                                     aGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgt    <  1:5035260/93‑1 (MQ=255)

                                                                                       ggCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgtcg  <  2:126466/93‑1 (MQ=255)

AGCGCCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCAGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCACGGCGTCG  >  NC_000913/52464‑52643

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