BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,013,656	T→C	M159T (ATG→ACG)	ygfJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,013,656	0	T	C	85.0%	699.7 / 115.6	299	M159T (ATG→ACG)	ygfJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (17/28); new base C (109/145); total (126/173)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

Only 100 of 299 total aligned reads displayed.

TCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGT  >  NC_000913/3013566‑3013745

tCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACgg                                                                                         <  2:3300894/93‑1 (MQ=255)

tCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACgg                                                                                         >  2:3656454/1‑93 (MQ=255)

  cccGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGga                                                                                       >  1:2024154/1‑93 (MQ=255)

   ccGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGgag                                                                                      <  2:3528323/93‑1 (MQ=255)

   ccGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGgag                                                                                      >  1:5778787/1‑93 (MQ=255)

      gCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATc                                                                                   <  1:5205367/93‑1 (MQ=255)

      gCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATc                                                                                   <  1:5433423/93‑1 (MQ=255)

      gCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATc                                                                                   <  1:6003368/93‑1 (MQ=255)

      gCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATc                                                                                   <  2:3766231/93‑1 (MQ=255)

       ccATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCa                                                                                  >  1:2733290/1‑93 (MQ=255)

       ccATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCa                                                                                  >  1:1618252/1‑93 (MQ=255)

          tCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCatta                                                                               >  2:2720509/1‑93 (MQ=255)

          tCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCatta                                                                               >  2:6231238/1‑93 (MQ=255)

          tCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCatta                                                                               >  1:5917845/1‑93 (MQ=255)

           ccGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCattat                                                                              >  1:55632/1‑93 (MQ=255)

           ccGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCattat                                                                              <  2:6144040/93‑1 (MQ=255)

           ccGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCattat                                                                              <  1:4971139/93‑1 (MQ=255)

            cGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCattatt                                                                             <  1:6276486/93‑1 (MQ=255)

               ttttAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCg                                                                          <  2:7362621/93‑1 (MQ=255)

               ttttAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCg                                                                          <  1:929217/93‑1 (MQ=255)

               ttttAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCg                                                                          >  1:39779/1‑93 (MQ=255)

               ttttAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCg                                                                          >  2:2048615/1‑93 (MQ=255)

               ttttAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCg                                                                          >  2:3711433/1‑93 (MQ=255)

                tttAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGt                                                                         <  2:4346689/93‑1 (MQ=255)

                tttAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGt                                                                         <  1:2229969/93‑1 (MQ=255)

                   aGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGa                                                                      <  2:2111404/93‑1 (MQ=255)

                   aGTGTCAAAACCATACCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGa                                                                      <  1:4830166/93‑1 (MQ=255)

                    gtgtCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGaa                                                                     >  1:1418421/1‑93 (MQ=255)

                       tCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATa                                                                  >  1:1199983/1‑93 (MQ=255)

                       tCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATa                                                                  <  1:4730844/93‑1 (MQ=255)

                       tCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATa                                                                  <  1:4537107/93‑1 (MQ=255)

                         aaaaCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGa                                                                <  1:6190053/93‑1 (MQ=255)

                         aaaaCCATCCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGa                                                                <  1:6592081/93‑1 (MQ=255)

                          aaaCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGaa                                                               <  1:2526491/93‑1 (MQ=255)

                           aaCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGaaa                                                              >  2:5531306/1‑93 (MQ=255)

                           aaCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGaaa                                                              <  1:5512139/93‑1 (MQ=255)

                           aaCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGaaa                                                              >  2:2429130/1‑93 (MQ=255)

                           aaCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGaaa                                                              <  2:6178076/93‑1 (MQ=255)

                           aaCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGaaa                                                              <  2:6649607/93‑1 (MQ=255)

                             ccATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAAt                                                            >  2:6269200/1‑93 (MQ=255)

                             ccATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAAt                                                            >  2:1380227/1‑93 (MQ=255)

                               aTGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGc                                                          <  2:5972946/93‑1 (MQ=255)

                               aTGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGc                                                          <  1:6929742/93‑1 (MQ=255)

                                tGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCa                                                         <  1:1798705/93‑1 (MQ=255)

                                tGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCa                                                         >  1:6408604/1‑93 (MQ=255)

                                tGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCa                                                         >  2:342864/1‑93 (MQ=255)

                                 gCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAg                                                        >  2:845004/1‑93 (MQ=255)

                                   cTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGaa                                                      >  1:4238165/1‑93 (MQ=255)

                                    tGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGaaa                                                     >  1:2987056/1‑93 (MQ=255)

                                    tGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGaaa                                                     >  2:6834104/1‑93 (MQ=255)

                                     gATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAaat                                                    >  1:2645014/1‑93 (MQ=255)

                                     gATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAaat                                                    >  2:1546829/1‑93 (MQ=255)

                                       tGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAaataa                                                  >  2:509415/1‑93 (MQ=255)

                                        gCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCGCTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAaataat                                                 >  2:1886469/1‑93 (MQ=255)

                                        gCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAaataat                                                 <  1:3498729/93‑1 (MQ=255)

                                          aGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAAttt                                               >  2:392822/1‑93 (MQ=255)

                                            gCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTa                                             <  2:3218386/93‑1 (MQ=255)

                                            gCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTa                                             <  1:344758/93‑1 (MQ=255)

                                              aaTCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGa                                           <  1:3607381/93‑1 (MQ=255)

                                               aTCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGat                                          >  1:7479400/1‑93 (MQ=255)

                                               aTCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGat                                          <  2:2966217/93‑1 (MQ=255)

                                               aTCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGat                                          <  1:5653239/93‑1 (MQ=255)

                                                tCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGata                                         <  2:5375078/93‑1 (MQ=255)

                                                tCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGata                                         >  1:1244246/1‑93 (MQ=255)

                                                 ccAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGatat                                        >  2:4852464/1‑93 (MQ=255)

                                                  cAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATAtt                                       >  2:5445874/1‑93 (MQ=255)

                                                      gACCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATa                                   <  1:3055208/93‑1 (MQ=255)

                                                       aCCCAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATAc                                  <  1:3687073/93‑1 (MQ=255)

                                                        cccAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACt                                 <  1:3892920/93‑1 (MQ=255)

                                                          cAATGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTcc                               <  2:4171417/93‑1 (MQ=255)

                                                            aTGTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCgg                             >  1:6769451/1‑93 (MQ=255)

                                                              gTCACCAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGgac                           <  1:2402223/93‑1 (MQ=255)

                                                                  ccAATATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACt                       >  2:1100341/1‑93 (MQ=255)

                                                                    aaTATGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACttt                     >  1:722001/1‑93 (MQ=255)

                                                                     atatGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTtta                    <  2:1915892/93‑1 (MQ=255)

                                                                       atGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTttatt                  <  1:3116694/93‑1 (MQ=255)

                                                                        tGCGTCAGGCACTCCTTATGGGAGGTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGGCTttatta                 >  1:6171899/1‑93 (MQ=255)

                                                                        tGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTttatta                 <  1:192925/93‑1 (MQ=255)

                                                                        tGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTttatta                 <  2:76100/93‑1 (MQ=255)

                                                                        tGCGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTttatta                 >  2:1590160/1‑93 (MQ=255)

                                                                         gCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTAc                <  1:7067436/93‑1 (MQ=255)

                                                                          cGTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACa               <  2:2265018/93‑1 (MQ=255)

                                                                           gTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAg              <  2:5835479/93‑1 (MQ=255)

                                                                           gTCAGGCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAg              >  2:7312983/1‑93 (MQ=255)

                                                                               ggCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGaa          <  2:477210/93‑1 (MQ=255)

                                                                               ggCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGaa          <  2:4051787/93‑1 (MQ=255)

                                                                               ggCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGaa          <  1:5255981/93‑1 (MQ=255)

                                                                               ggCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGaa          >  2:2792210/1‑93 (MQ=255)

                                                                               ggCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGaa          >  1:4359367/1‑93 (MQ=255)

                                                                               ggCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGaa          <  2:7193016/93‑1 (MQ=255)

                                                                                gCACTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAg         <  2:1448420/93‑1 (MQ=255)

                                                                                  aCTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGa       >  2:3000429/1‑93 (MQ=255)

                                                                                  aCTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGa       <  1:5038203/93‑1 (MQ=255)

                                                                                   cTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGaa      >  2:7087931/1‑93 (MQ=255)

                                                                                   cTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGaa      >  2:3833057/1‑93 (MQ=255)

                                                                                   cTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGaa      <  2:6061884/93‑1 (MQ=255)

                                                                                   cTCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGaa      <  2:5253253/93‑1 (MQ=255)

                                                                                    tCCTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGTAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGaaa     >  2:5194102/1‑93 (MQ=255)

                                                                                      cTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAgg   <  2:6565160/93‑1 (MQ=255)

                                                                                      cTTACGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAgg   <  2:232713/93‑1 (MQ=255)

TCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGT  >  NC_000913/3013566‑3013745

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