BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,991,845	T→C	N62N (AAT→AAC)	b2853 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,991,845	0	T	C	99.8%	2027.2 / ‑5.2	613	N62N (AAT→AAC)	b2853	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (326/285); total (327/286)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

Only 100 of 621 total aligned reads displayed.

CACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCAT  >  NC_000913/2991755‑2991937

   cAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCTGAAGATAACCCTAt                                                                                         >  2:1566920/1‑93 (MQ=255)

   cAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAt                                                                                         <  2:3701910/93‑1 (MQ=255)

    aGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAtc                                                                                        >  2:9654817/1‑93 (MQ=255)

    aGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAtc                                                                                        >  1:2778005/1‑93 (MQ=255)

    aGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAtc                                                                                        <  2:7827539/93‑1 (MQ=255)

     gTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAtct                                                                                       >  2:1440062/1‑93 (MQ=255)

      tGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTAtctc                                                                                      <  2:3991868/93‑1 (MQ=255)

          aCTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCt                                                                                  >  2:5641646/1‑93 (MQ=255)

          aCTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCt                                                                                  <  2:6002601/93‑1 (MQ=255)

           cTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCtt                                                                                 >  2:1144524/1‑93 (MQ=255)

            tGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTa                                                                                >  1:4120555/1‑93 (MQ=255)

            tGATACAGATCAAATACAAGCATTTACAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTAGTCAGAAGATAACCCTATCTCACCCTa                                                                                >  2:8932523/1‑93 (MQ=255)

             gATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTAc                                                                               <  2:8735006/93‑1 (MQ=255)

               tACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACta                                                                             <  1:5430970/93‑1 (MQ=255)

               tACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACta                                                                             >  1:3300621/1‑93 (MQ=255)

               tACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACta                                                                             >  1:1520541/1‑93 (MQ=255)

                 cAGATCAACTACAAGCATTTATAGTTGCAACATTGATGGCGCCTTTTCAAAATTATATGTATTCAGAAGATAACCCTATCTCACCTTACtata                                                                           <  1:5542082/93‑1 (MQ=255)

                   gATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATaaa                                                                         >  2:7761527/1‑93 (MQ=255)

                     tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAAt                                                                       >  1:8472378/1‑93 (MQ=255)

                     tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAAt                                                                       <  2:1850539/93‑1 (MQ=255)

                     tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAAt                                                                       <  2:3299074/93‑1 (MQ=255)

                       aaaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTg                                                                     <  1:3192074/93‑1 (MQ=255)

                         aTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAg                                                                   <  1:9655108/93‑1 (MQ=255)

                         aTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAg                                                                   >  2:4728038/1‑93 (MQ=255)

                          tACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGt                                                                  <  1:1308062/93‑1 (MQ=255)

                           aCAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTg                                                                 >  1:3935897/1‑93 (MQ=255)

                            cAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTgg                                                                >  1:7520996/1‑93 (MQ=255)

                            cAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTgg                                                                >  1:6832977/1‑93 (MQ=255)

                             aaGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGt                                                               >  2:6950782/1‑93 (MQ=255)

                              aGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTa                                                              >  1:9821628/1‑93 (MQ=255)

                              aGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTa                                                              <  2:6534270/93‑1 (MQ=255)

                               gCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAg                                                             <  2:9503989/93‑1 (MQ=255)

                               gCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAg                                                             <  1:6563514/93‑1 (MQ=255)

                                cATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGa                                                            >  2:5593125/1‑93 (MQ=255)

                                cATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGa                                                            <  1:2170696/93‑1 (MQ=255)

                                  tttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACg                                                          >  2:7371987/1‑93 (MQ=255)

                                  tttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACg                                                          <  1:896790/93‑1 (MQ=255)

                                  tttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACg                                                          >  2:9246513/1‑93 (MQ=255)

                                   ttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGa                                                         >  1:2969444/1‑93 (MQ=255)

                                   ttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGa                                                         >  2:10024801/1‑93 (MQ=255)

                                    tataGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGAt                                                        >  1:2022086/1‑93 (MQ=255)

                                    tataGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGAt                                                        <  1:6273512/93‑1 (MQ=255)

                                    taaaGTTTCACCATGGATGGCGCCTTTTCAAAATGATACGTATTCAGAAGATAACCCTATCTCCCCTTACTAT‑AAATTGAGTGGTAGACGAtt                                                       <  1:8711607/90‑1 (MQ=255)

                                     ataGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGAtt                                                       <  2:6329863/93‑1 (MQ=255)

                                     ataGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGAtt                                                       >  2:7173990/1‑93 (MQ=255)

                                      taGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTg                                                      >  1:7284121/1‑93 (MQ=255)

                                       aGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGa                                                     >  2:2176160/1‑93 (MQ=255)

                                           tCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaa                                                 >  1:2543077/1‑93 (MQ=255)

                                           tCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaa                                                 <  1:446593/93‑1 (MQ=255)

                                             aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaat                                               >  1:6710641/1‑93 (MQ=255)

                                              aCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaata                                              <  2:4558678/93‑1 (MQ=255)

                                                atggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaataaa                                            <  1:9489003/93‑1 (MQ=255)

                                                 tggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaataaat                                           <  1:1140033/93‑1 (MQ=255)

                                                 tggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCaaataaat                                           >  2:4485797/1‑93 (MQ=255)

                                                   gatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCa                                         >  1:5750477/1‑93 (MQ=255)

                                                   gatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCa                                         <  1:7929096/93‑1 (MQ=255)

                                                   gatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCa                                         <  2:4198696/93‑1 (MQ=255)

                                                    atggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCAt                                        <  1:360120/93‑1 (MQ=255)

                                                    atggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCAt                                        <  1:4948513/93‑1 (MQ=255)

                                                      ggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATaa                                      <  1:1829128/93‑1 (MQ=255)

                                                      ggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATaa                                      >  1:9066942/1‑93 (MQ=255)

                                                      ggCGCCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATaa                                      <  2:2918467/93‑1 (MQ=255)

                                                       gcgcCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAAc                                     <  2:8202297/93‑1 (MQ=255)

                                                         gcCTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACga                                   >  2:4417072/1‑93 (MQ=255)

                                                          ccTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACgag                                  <  2:8404530/93‑1 (MQ=255)

                                                           cTTTTCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACgaga                                 <  1:5359837/93‑1 (MQ=255)

                                                            ttttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGAt                                <  2:511750/93‑1 (MQ=255)

                                                            ttttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGAt                                >  1:6020897/1‑93 (MQ=255)

                                                            ttttCAAAATGATATGTATTCAGAACATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATAGACCAAATAAATCATAACGAGAt                                >  2:5499274/1‑93 (MQ=255)

                                                             tttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATa                               <  2:8076866/93‑1 (MQ=255)

                                                             tttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATa                               >  1:3450171/1‑93 (MQ=255)

                                                              ttCAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATAc                              <  1:8163156/93‑1 (MQ=255)

                                                                cAAAATGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtt                            <  1:3708973/93‑1 (MQ=255)

                                                                 aaaaTGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACttt                           <  2:2618105/93‑1 (MQ=255)

                                                                  aaaTGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtttt                          <  2:4815546/93‑1 (MQ=255)

                                                                    aTGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtat                        >  1:9350083/1‑93 (MQ=255)

                                                                     tGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtata                       >  2:2452993/1‑93 (MQ=255)

                                                                     tGATATGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtata                       <  2:3429478/93‑1 (MQ=255)

                                                                        tatGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataata                    >  1:6369813/1‑93 (MQ=255)

                                                                         atGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                   <  1:8299852/93‑1 (MQ=255)

                                                                         atGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                   <  1:5231638/93‑1 (MQ=255)

                                                                         atGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                   >  1:5110799/1‑93 (MQ=255)

                                                                         atGTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                   >  1:428438/1‑93 (MQ=255)

                                                                           gTATTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa                 <  1:7154656/93‑1 (MQ=255)

                                                                             aTTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc               <  1:5805942/93‑1 (MQ=255)

                                                                             aTTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc               <  2:9365508/93‑1 (MQ=255)

                                                                             aTTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc               <  1:7748244/93‑1 (MQ=255)

                                                                             aTTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc               <  2:8544277/93‑1 (MQ=255)

                                                                             aTTCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc               >  2:5305528/1‑93 (MQ=255)

                                                                              ttCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtct              >  1:6966523/1‑93 (MQ=255)

                                                                              ttCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtct              >  1:1643094/1‑93 (MQ=255)

                                                                               tCAGAAGATAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtctc             <  2:2674737/93‑1 (MQ=255)

                                                                                 agaagaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt           <  1:10009094/93‑1 (MQ=255)

                                                                                  gaagaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATTAATCATAACGAGATACTTTTATAATAGTATCTCGtt          >  1:9956230/1‑93 (MQ=255)

                                                                                  gaagaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGtt          >  2:7476418/1‑93 (MQ=255)

                                                                                   aagaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTc         >  2:6730341/1‑93 (MQ=255)

                                                                                    agaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTcc        <  1:5019043/93‑1 (MQ=255)

                                                                                    agaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTcc        <  1:6141714/93‑1 (MQ=255)

                                                                                     gaTAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCt       >  1:4781705/1‑93 (MQ=255)

                                                                                       tAACCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt     >  1:4654118/1‑93 (MQ=255)

CACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCAT  >  NC_000913/2991755‑2991937

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