BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,629,183	T→C	K346K (AAA→AAG)	yhiJ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,629,183	0	T	C	82.7%	778.8 / 154.6	352	K346K (AAA→AAG)	yhiJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (34/27); new base C (145/146); total (179/173)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 354 total aligned reads displayed.

TCCCATGATGATTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCATTGATACTCATTA  >  NC_000913/3629092‑3629274

tCCCATGATGATTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGtt                                                                                            <  1:4131649/93‑1 (MQ=255)

  ccATGATGATTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTc                                                                                          >  2:7504550/1‑93 (MQ=255)

   cATGATGATTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCa                                                                                         >  1:654849/1‑93 (MQ=255)

   cATGATGATTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCa                                                                                         >  2:5321048/1‑93 (MQ=255)

   cATGATGATTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCa                                                                                         <  1:3862341/93‑1 (MQ=255)

   cATGATGATTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCa                                                                                         <  1:3722927/93‑1 (MQ=255)

   cATGATGATTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCa                                                                                         >  2:177747/1‑93 (MQ=255)

    atgatgATTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCAt                                                                                        >  1:2898279/1‑93 (MQ=255)

     tgatgaTTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCATa                                                                                       <  2:1028880/93‑1 (MQ=255)

     tgatgaTTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATa                                                                                       >  2:4751753/1‑93 (MQ=255)

     tgatgaTTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATa                                                                                       >  1:7309077/1‑93 (MQ=255)

     tgatgaTTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATa                                                                                       <  1:744194/93‑1 (MQ=255)

      gatgatTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCATAc                                                                                      <  2:7097012/93‑1 (MQ=255)

      gatgatTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATAc                                                                                      >  2:4567696/1‑93 (MQ=255)

       atgatTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATAcg                                                                                     <  2:5657187/93‑1 (MQ=255)

        tgatTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATAcgc                                                                                    >  1:4185640/1‑93 (MQ=255)

         gatTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCATAcgcg                                                                                   <  2:6559967/93‑1 (MQ=255)

          atTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCATACGCGt                                                                                  <  2:3724608/93‑1 (MQ=255)

          atTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGt                                                                                  >  1:1561606/1‑93 (MQ=255)

          atTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGt                                                                                  <  1:7509870/93‑1 (MQ=255)

           tttGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTg                                                                                 <  1:4077233/93‑1 (MQ=255)

            ttGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGa                                                                                >  1:5745794/1‑93 (MQ=255)

              gCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCGGAGTAGGTATCCTCATCGAGTGGTTTCATACGCGTGAAg                                                                              >  2:1236751/1‑93 (MQ=255)

               ccATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGt                                                                             >  2:3096065/1‑93 (MQ=255)

                cATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCATACGCGTGAAGtt                                                                            >  2:2066531/1‑93 (MQ=255)

                cATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCATACGCGTGAAGtt                                                                            <  2:6050345/93‑1 (MQ=255)

                 attattTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTc                                                                           <  1:3018936/93‑1 (MQ=255)

                    attTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGc                                                                        <  1:6868976/93‑1 (MQ=255)

                     ttttGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCa                                                                       >  1:1094261/1‑93 (MQ=255)

                      tttGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCat                                                                      <  1:6639568/93‑1 (MQ=255)

                       ttGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCata                                                                     <  2:5442010/93‑1 (MQ=255)

                            tttGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCATACGCGTGAAGTTCTGCATATAGtt                                                                >  1:1687437/1‑93 (MQ=255)

                            tttGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGtt                                                                <  1:3931161/93‑1 (MQ=255)

                                gggTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCa                                                            <  2:4232653/93‑1 (MQ=255)

                                gggTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCa                                                            >  1:487639/1‑93 (MQ=255)

                                 ggTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCAt                                                           <  1:6939060/93‑1 (MQ=255)

                                  gTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCAtt                                                          <  1:2106256/93‑1 (MQ=255)

                                   tCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCATACGCGTGAAGTTCTGCATATAGTTACCATTg                                                         >  1:5625741/1‑93 (MQ=255)

                                   tCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTTAAGTTCTGTATATAGTTACCAttc                                                         >  2:7317195/1‑92 (MQ=255)

                                   tCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTg                                                         >  2:5044512/1‑93 (MQ=255)

                                      tAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCaa                                                      >  1:2285292/1‑93 (MQ=255)

                                      tAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCaa                                                      <  2:2752306/93‑1 (MQ=255)

                                      tAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCaa                                                      >  1:705947/1‑93 (MQ=255)

                                       aGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAAt                                                     >  1:6801452/1‑93 (MQ=255)

                                       aGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAAt                                                     <  2:5490944/93‑1 (MQ=255)

                                        gAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAAtt                                                    >  1:566263/1‑93 (MQ=255)

                                        gAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAAtt                                                    <  1:6744718/93‑1 (MQ=255)

                                         aaaCATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTc                                                   <  2:625143/93‑1 (MQ=255)

                                          aaCATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCg                                                  <  1:4982205/93‑1 (MQ=255)

                                          aaCATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCg                                                  <  1:5315277/93‑1 (MQ=255)

                                            cATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGtt                                                <  2:452492/93‑1 (MQ=255)

                                             aTACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGttt                                               <  1:3281771/93‑1 (MQ=255)

                                              tACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTc                                              <  1:6262853/93‑1 (MQ=255)

                                              tACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTc                                              >  1:4673791/1‑93 (MQ=255)

                                               aCAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCa                                             <  1:1247119/93‑1 (MQ=255)

                                               aCAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCa                                             <  2:5825986/93‑1 (MQ=255)

                                                 aGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATg                                           <  1:1927086/93‑1 (MQ=255)

                                                  ggATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGa                                          >  1:1369001/1‑93 (MQ=255)

                                                  ggATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGa                                          <  2:1682074/93‑1 (MQ=255)

                                                  ggATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGa                                          >  2:4843363/1‑93 (MQ=255)

                                                     tAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGAttt                                       >  2:926757/1‑93 (MQ=255)

                                                        tCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCa                                    >  2:3245165/1‑93 (MQ=255)

                                                         cTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAg                                   <  2:3970420/93‑1 (MQ=255)

                                                          ttttATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGc                                  >  1:1008581/1‑93 (MQ=255)

                                                          ttttATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGc                                  >  1:366646/1‑93 (MQ=255)

                                                           tttATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTCCCATTGCAATTCGTTTCATGATTTCCAGCt                                 <  1:2811503/93‑1 (MQ=255)

                                                           tttATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCt                                 <  2:2387894/93‑1 (MQ=255)

                                                            ttATCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTa                                >  1:4418488/1‑93 (MQ=255)

                                                              aTCCCAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAAt                              >  1:1492031/1‑93 (MQ=255)

                                                                 ccAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACa                           <  2:1564380/93‑1 (MQ=255)

                                                                 ccAGAGTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACa                           >  1:2001029/1‑93 (MQ=255)

                                                                   agagTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATg                         <  1:5942143/93‑1 (MQ=255)

                                                                    gagTAGGTATCCTCATCGAGTGGTTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGa                        >  2:3178712/1‑93 (MQ=255)

                                                                    gagTAGGTATCCTCATCGAGTGGTTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGa                        <  1:6354393/93‑1 (MQ=255)

                                                                     agTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGaa                       >  1:5261446/1‑93 (MQ=255)

                                                                     agTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGaa                       <  1:2705684/93‑1 (MQ=255)

                                                                      gTAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAAc                      <  2:7367852/93‑1 (MQ=255)

                                                                       tAGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACa                     >  1:1762249/1‑93 (MQ=255)

                                                                        aGGTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACAt                    <  2:2548019/93‑1 (MQ=255)

                                                                          gTATCCTCATCGAGTGGTTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACAttt                  <  2:1930084/93‑1 (MQ=255)

                                                                          gTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACAttt                  <  2:688808/93‑1 (MQ=255)

                                                                          gTATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACAttt                  >  1:914790/1‑93 (MQ=255)

                                                                           tATCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACAtttt                 >  1:5433793/1‑93 (MQ=255)

                                                                            aTCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTc                <  2:6983352/93‑1 (MQ=255)

                                                                            aTCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTc                >  2:5997835/1‑93 (MQ=255)

                                                                             tCCTCATCGAGTGGTTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCa               <  2:141363/93‑1 (MQ=255)

                                                                             tCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCa               <  1:883027/93‑1 (MQ=255)

                                                                             tCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCa               >  2:1743363/1‑93 (MQ=255)

                                                                             tCCTCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCa               >  2:3396396/1‑93 (MQ=255)

                                                                                tCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCATTg            >  1:2390953/1‑93 (MQ=255)

                                                                                tCATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCATTg            >  1:5149702/1‑93 (MQ=255)

                                                                                 cATCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCATTGa           <  2:3494674/93‑1 (MQ=255)

                                                                                  aTCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCATTGAt          <  2:2356611/93‑1 (MQ=255)

                                                                                  aTCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACACGAACATTTTCATTGAt          >  1:3345961/1‑93 (MQ=255)

                                                                                   tCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCATTGATa         >  2:2982736/1‑93 (MQ=255)

                                                                                   tCGAGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCATTGATa         <  1:5821757/93‑1 (MQ=255)

                                                                                      aGTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCATTGATACTc      >  2:671053/1‑93 (MQ=255)

                                                                                       gTGGCTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCATTGATACTCa     >  2:7240903/1‑93 (MQ=255)

                                                                                          gTTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCATTGATACTCATTa  <  2:2213577/93‑1 (MQ=255)

                                                                                          gTTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCATTGATACTCATTa  <  2:7310150/93‑1 (MQ=255)

TCCCATGATGATTTGCCATTATTTTGTATTTGGGGTCATAGAAACATACAGGATAATCTTTTATCCCAGAGTAGGTATCCTCATCGAGTGGTTTCATACGCGTGAAGTTCTGCATATAGTTACCATTGCAATTCGTTTCATGATTTCCAGCTAATACATGAACATTTTCATTGATACTCATTA  >  NC_000913/3629092‑3629274

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